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Abstract 



In this thesis properties and the origin of black hole entropy are investigated 
from various points of view. First, laws of black hole thermodynamics are reviewed. 
In particular, the first and generalized second laws are investigated in detail. It is in 
these laws that the black hole entropy plays key roles. Next, three candidates for the 
origin of the black hole entropy are analyzed: the D-brane statistical-mechanics, the 
brick wall model, and the entanglement thermodynamics. Finally, discussions are 
given on semiclassical consistencies of the brick wall model and the entanglement 
thermodynamics and on the information loss problem. 
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Chapter 1 

Introduction 



Thermodynamics describes behavior of coarse-grained or averaged quantities of a 
system with a large number of physical degrees of freedom. The behavior is traced 
by a small number of parameters. Mathematically, the microscopic description of 
thermodynamics, statistical mechanics, is grounded by the ergodic hypothesis. On 
the other hand, in the theory of black holes, the no hair theorem [jj] allows us 
to describe a stationary black hole by a small number of parameters. Since the 
cosmic censorship conjecture Q combined with the singularity theorem jl| predicts 
inevitable occurrence of black holes, it seems that the no hair theorem plays the 
same role as the ergodic hypothesis plays in thermodynamics. 

In fact, it is well known that black holes have many properties analogous to 
those of thermodynamics. Those are as a whole called black hole thermodynamics. 
In particular, four laws of black holes combined with the generalized second law 
make up a main framework of the black hole thermodynamics. In these laws, black 
hole entropy defined as follows plays an important role. 

(1.1, 

where A is area of black hole horizon. Moreover, it was suggested that the black 
hole entropy, or the horizon area, is an adiabatic invariant B and that it can be 
used as a potential function in catastrophe theory to judge stability of black hole 
solutions ||. The formula (1.1) for black hole entropy is often called Bekenstein- 



Hawking formula since the concept of black hole entropy was first introduced by 
Bekenstein |j| as a quantity proportional to the horizon area and the proportional- 
ity coefficient was fixed by Hawking's discovery of thermal radiation from a black 
hole (see arguments below). He showed that a black hole radiates thermal radi- 
ation with temperature given by 

, „ %K , 

k B T BH = ^— , 1.2 
2nc 

where k is the surface gravity of a background black hole. This thermal radia- 
tion and its temperature are called Hawking radiation and Hawking temperature, 
respectively. 

Let us recall basic properties of the black-hole thermodynamics by taking a 
simple example. We consider a one-parameter family of Schwarzschild black holes 
parameterized by the mass Mbh- We assume that a relation analogous to the first 
law of thermodynamics holds for a black-hole system. In the present example, there 
is only one parameter M B h characterizing a black hole. Therefore, this relation 
should be of the simplest form 

SEbh = TbhSSbh, (1-3) 
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where Ebh, Sbh and Tbh are quantities that are identified with the energy, the 



entropy and the temperature of a black hole, respectively. The relation Eq.(1.3) is 
called the first law of the black- hole thermodynamics Thus, if two of the quan- 
tities Ebh, Sbh and Tbh are given, ^f^g. ( ]1 .3| ) determines the remaining quantity. 

In the present example, Mbh is the only parameter characterizing the family 
of black holes. Therefore the simplest combination which yields the dimension of 
energy is 

E BH = M BH c 2 . (1.4) 

This is the energy of the black hole. 

There is also a natural choice for Tbh @- Hawking showed that a black hole 
with surface gravity n emits thermal radiation of a quantum matter field (which 
plays the role of a thermometer) at temperature given by (|l.2|). Moreover, as shown 



in section 2.3, if a matter field in a thermal-equilibrium state at some temperature is 
scattered by a black hole, then it always becomes closer to the thermal-equilibrium 
state at the Hawking temperature fll.2|) M, nn. Thus it is natural to define the 



temperature of a Schwarzschild black hole with mass Mbh by 



Tic 3 

kBTBH = 8^' (L5) 



since n = c 4 /4GM B h 




From Eqs.(|l.3|)-(|l.5|), we get an expression for Sbh as 



S ™ = k -§& A+C > ^ 

where A = 16ttG 2 M B h/ c4 is the area of the event horizon and C is some constant. 
Since a value of C is not essential in our discussions, we shall set hereafter 



C* = 0. (1.7) 



Note that Eq. (L6) with (1/7) is a special case of the Bekenstein-Hawking formula 



It is well-known that classically the area of the event horizon does not decrease 
in time (the area law |p^| or the second law of black hole) just as the ordinary ther- 



modynamical entropy. The Bekenstein-Hawking formula (1.1) looks reasonable in 
this sense. Indeed this observation was the original motivation for the introduction 
of the black-hole entropy [^J. Moreover, when quantum effects are taken into ac- 
count, it is believed that a sum of the black hole entropy and matter entropy does 
not decrease (the generalized second law). 

The zeroth law of black hole thermodynamics states that surface gravity of a 
Killing horizon is constant throughout the horizon. This supports our choice of the 
black hole temperature. (For a proof of the zeroth law, see Refs. || |ll|, Q.) 

At this stage we would like to point out that for a black hole the third law 
does not hold in the sense of Planck: Sbh — ► oo as Tbh — ¥ for the family of 
Schwarzschild black holes, irrespective of the choice of the value for C as is seen 



from Eqs.(1.5) and (1.6). Rather, the third law does hold in the sense of Nernst: it 
is impossible by any process, no matter how idealized, to reduce n to zero in a finite 
sequence of operations [||. (See Ref. |L4| for a precise expression and a proof.) 

Thermodynamics has a well-established microscopic description: the quantum 
statistical mechanics. In the thermodynamical description, information on each mi- 
croscopic degree of freedom is lost, and only macroscopic variables are concerned. 
However, the number of all microscopic degrees of freedom is implemented in a 
macroscopic variable: entropy S is related to the number of all consistent micro- 
scopic states N as 

S = k B \nN. (1.8) 
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In analogy, it is expected that there might be a microscopic description of the black 
hole thermodynamics, too. In particular, it is widely believed that the black hole 
entropy might be related to a number of microscopic states. Since the microscopic 
description seems to require a quantum theory of gravity, detailed investigations of 
the black hole entropy should contribute a lot toward construction of the theory of 
quantum gravity. This is one among the several reasons why the origin of the black 
hole entropy needs to be understood at the fundamental level. 

Another strong motivation to investigate the black hole entropy is the so-called 
information loss problem. Hawking argued that, if a black hole is formed by gravi- 
tational collapse, then evolution of quantum fields becomes non-unitary because of 
evaporation of the black hole due to the Hawking radiation p5| . This means that 
some information is lost in the process of the black hole evaporation. Moreover, 
this suggests that the conventional field-theoretical approach may be useless for 
the purpose of construction of the theory of quantum gravity since the field theory 
is based on the unitarity. Hence, the evaporation of a black hole makes people, 
who wish to construct a unitary theory of quantum gravity, be in difficulties: if 
the evaporation of a black hole would actually occur and lead to the information 
loss, then they would be obliged to give up the unitarity. Thus, we have to clarify 
whether information is lost or not due to the Hawking radiation in order to take a 
step forward. This problem is called the information loss problem. On the other 
hand, since entropy is strongly connected with information in the theory of infor- 
mation, it is natural to expect that the black hole entropy might be related to some 
information. Therefore, investigations of the origin of the black hole entropy seem 
to provide important insight toward the information loss problem. 

In these senses, the origin of the black hole entropy is one of the most important 
issues at the present stage of black hole physics. 

Recently a microscopic derivation of the black hole entropy was given in su- 
perstring theory pi FT?]] by using the so-called D-brane technology |lq|. In this 



approach, as will be shown in section 3.1, the black hole entropy is identified with 



the logarithm of the number of states of massless strings attached to D-branes, 
with D-brane configuration and total momentum of the strings along a compacti- 
fied direction fixed to be consistent with the corresponding black hole jl9|, ^o) . The 
analysis along this line was extended to the so-called M-theory In particular, 
by using a conjectured correspondence (the Matrix theory) between the M-theory 
in the infinite momentum frame and a 10-dimcnsional U(N) supersymmetric Yang- 
Mills theory dimensionally reduced to (0 + l)-dimension with N — > oo f2^| , the black 
hole entropy was calculated by means of the Yang-Mills theory. The result gives 
the correct Bekenstein-Hawking entropy for BPS black holes and their low lying 
excitations Moreover, in Ref. the black hole entropy of a Schwarzschild 

black hole was derived in the Matrix theory up to a constant of order unity. On the 
other hand, in loop quantum gravity |2j|, black hole entropy was identified with the 
logarithm of the number of different spin-network states for a fixed eigenvalue of 
the area operator pq] . The result coincides with the Bekenstein-Hawking entropy 
up to a constant of order unity. 

The derivations in these candidate theories of quantum gravity depend strongly 
on details of the theories. In this sense, the success of the derivations can be 
considered as non-trivial consistency checks of the theories. However, it is believed 
that proportionality of the black hole entropy to horizon area is more universal 
and does not depend on details of the theory. Hence, one should be able to give 
a statistical or thermodynamical derivation of the black hole entropy, which does 
not depend on details of theory, while we are proceeding with theory-dependent 
derivations of it by using the well-established candidate theories of quantum gravity. 

There were many attempts to explain the origin of the black hole entropy be- 
sides the above theory-dependent approaches. (See Ref. 27 for an up-to-date re- 



CHAPTER 1. INTRODUCTION 



6 



view.) For example, in Euclidean gravity the black hole entropy is associated with 
the topology of an instanton which corresponds to a black hole ]28|, |3(], |l| q 
Wald [|53| defined the black hole entropy as a Noether charge associated with a 
bifurcating Killing horizon^] (See section ^j].); 'tHooft |35| identified the black hole 
entropy with the statistical entropy of a thermal gas of quantum particles with a 
mirror-like boundary just outside the horizon (This model is called the brick wall 
model and is analyzed in detail in section pT^ ) ; Pretorius et al. |3(| identified the 
black hole entropy with the thermodynamical entropy of a shell in thermal equilib- 
rium with acceleration radiation due to the shell's gravity in the limit that the shell 
forms a black hole. 

There remains another strong candidate for the statistical origin of the black 
hole entropy, called entanglement entropy |?7|, p9| . It is a statistical entropy 
measuring the information loss due to a spatial division of a system ||]J. The 
entanglement entropy is based only on the spatial division, and can be defined 
independently of the theory, although explicit calculations in the literature are 
dependent on the model employed. Moreover, as will be explained in section |3.3| , it 
is expected independently of the details of the theory that the entanglement entropy 
is proportional to the area of the boundary of the spatial division. In this sense, 
the entanglement entropy is considered to be a strong candidate for the statistical 
origin of the black hole entropy. 

In chapter p|, laws of black hole thermodynamics are reviewed. In particular 
the first laws of black hole statics and dynamics, and the generalized second law 
are studied in detail. The first law derived in the above simple arguments on 
Schwarzschild black holes relates changes of physical quantities of stationary black 
holes corresponding to a variation in a space of stationary black hole solutions. In 
this sense we can call it the first law of black hole statics. Historically, the first law of 
black hole statics is derived in Ref. |8| in general relativity, and is extended by Wald 
to a general covariant theory of gravity p3| . In section 2.1 we re-analyze the first 

we consider a generalization 



law of Wald in detail following Ref. [j40|. In section 2.2 
of the first law to a purely dynamical situation |4l|. 



42 



We call the dynamical 
a proof of the generalized 



version the first law of black hole dynamics. In sectio n |2.3 
second law is given for a quasi-stationary black hole |l(f| . 

In chapter |3J, three ca ndid ates for the origin of the black hole entropy are ana- 
lyzed in detail. In section a microscopic derivation of black hole entropy by the 
D-brane technology is shown. We consider a 5-dimensional black hole solution in 
the low energy effective theory of Type IIB superstring. This black hole is, in fact, a 
black brane in 10-dimensional sense and can be interpreted as a configuration of D- 
branes wrapped on T 5 = T 4 x S 1 . We calculate statistical- mechanical entropy and 
temperature of open strings on the D-branes and compare them with the Bekenstein- 
Hawking entropy and the Hawking temperature of the original 5-dimensional black 
hole (43). In section 3.2 we re-examine the brick wall model in detail and solve 
problems concerning this model f44|. In section 3.2 we construct a thermodynamics 
(entanglement thermodynamics [ [45| , f4(| ) which includes the entanglement entropy 
as the entropy variable, for a massless scalar field in Minkowski, Schwarzschild and 
Reissner-Nordstrom spacetimes to understand the statistical origin of black-hole 
thermodynamics. In section 3.4 a new interpretation of entanglement entropy is 
proposed 

Chapter W is devoted to a summary of this thesis, discussions on semiclassical 
consistencies and the information loss problem, and speculations. 



1 The 1-loop correction to the black hole entropy was also calculated and compared with the 
brick wall model and the conical singularity method [j^] . 

2 Relations to the approach by Euclidean gravity was investigated in Ref. ]34|. 



Chapter 2 



Laws of black hole 
thermodynamics 

2.1 The first law of black hole statics 

In Ref. [Q, the first law of black hole mechanics was derived not only in general 
relativity but also in a general covariant theory of gravity for stationary variations 
around a stationary black hole. It is formulated as a relation among variations of 
those quantities such as energy, angular momentum and entropy, each of which is 
defined in terms of a Noether charge. The first law was extended to non-stationary 
variations around a stationary black hole in Ref. p8| ]. 

The Noether charge form of the first law has many advantages over the original 
first law of Ref. ||. For example, it gives a general method to calculate stationary 
black hole entropy in general covariant theories of gravity [pf8|| ; it connects various 
Euclidean methods for computing black hole entropy ||34f| ; it suggests a possibility 
of defining entropy of non-stationary black holes jl8j |49|; etc. 

However, in their derivation there are several issues to be discussed in more 
detail. 

(a) In Ref. |p3 , unperturbed and perturbed stationary black holes are identified 
so that horizon generator Killing fields with unit surface gravity coincide in 
a neighborhood of the horizons and that stationary Killing fields and axial 
Killing fields coincide in a neighborhood of infinity. This corresponds to tak- 
ing a certain gauge condition in linear perturbation theory. For a complete 
understanding of the first law, we have to clarify whether such a gauge con- 
dition can be imposed or not. If it can, then we wish to know whether such a 
gauge condition is necessary. Note that, on the contrary, the original deriva- 
tion in general relativity by Bardeen, Carter and Hawking Q is based on a 
gauge condition such that the stationary Killing fields and the axial Killing 
fields coincide everywhere on a spacelike hypersurface whose boundary is a 
union of a horizon cross section and spatial infinity. 

(b) In Ref. |l8|, the first law is extended to non-stationary perturbations around 
a stationary black hole. In the derivation, change of black hole entropy is 
calculated on a (n — 2)-surface, which is a bifurcation surface for an unper- 
turbed black hole, but which is not a cross section of an event (nor apparent) 
horizon for a perturbed non-stationary black hole in general. Does this mean 
that black hole entropy would be assigned to a surface which is not a horizon 
cross section for a non-stationary black hole? It seems more natural to assign 
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black hole entropy to a horizon cross section also for a non-stationary black 
hole. 

In this section these two issues are discussed and it is concluded that there are no 
difficulties in the derivation of the Noether charge form of the first law for both 
stationary and non-stationary perturbations about a stationary black hole. In its 
course, we give an alternative derivation of the first law based on a variation in 
which a horizon generator Killing field with unit surface gravity is fixed. 



In subsection 2.1.1 gauge conditions are analyzed. In subsection 2.1.2 the first 



law of black holes is derived for stationary variations around a stationary black hole. 



In subsection 2.1.3 the derivation is extended to non-stationary variations around 



a stationary black hole. 

2.1.1 Gauge conditions 

Consider a stationary black hole in n-dimensions, which has a bifurcating Killing 
horizon. Let £° be a generator Killing field of the Killing horizon, which is nor- 
malized as £° = t a + il^cpfy and £ be the bifurcation surface. Here, t a is the 
stationary Killing field with unit norm at infinity, {f®^} (p — 1, 2, ■ • ■) is a family 

of axial Killing fields, and {D,^} is a family of constants (angular velocities). Let 
k be the surface gravity corresponding to £ a : 

£ fe v 6 r = (2.i) 

on the horizon. 

Now let us show that it is not possible in general to impose a gauge condition 
such that <5£ a = in a neighborhood of the bifurcation surface. For this purpose 
we shall temporarily assume that <5£ a — and show a contradiction. 

On E, the covariant derivative of £° is given by 

V„r = (2.2) 
where e ab is binormal to S. However, the variation of the l.h.s. is zero: 

s(v b a = srie = o (2.3) 

since £ a = 0, where 5T% C is given by 

5?bc = \g ad {V c 5g db + V b 5g dc - V d 5g bc ). (2.4) 



Hence, 



5e b a = --e b a . (2.5) 

K 

b 



Substituting this into the identity S(e b a e a ) = 0, we obtain 

= 5(e b "e\) = -^. (2.6) 

K 

Thus, the assumption S^ a = leads to 5k = 0, which implies, for example, that 
5M = for the vacuum general relativity in a static case, where M is mass of 
Schwarzschild black holes. This peculiar behavior can be understood as appearance 
of a coordinate singularity at the bifurcation surface of a coordinate fixed by the 
gauge condition 5^ a = since in the above argument finiteness of <5r^ c has been 
assumed implicitly. Therefore, it is impossible to impose the condition <5£ a = in 
a neighborhood of the bifurcation surface whenever 5k 0. 
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As mentioned above, the original derivation of the first law in Ref. g adopt the 
gauge condition 5t a = 5ip a — 0. This leads to <5£ a = when <5f2 = (for example, 
when we consider static black holes). Of course, in Ref. ||, a general horizon cross 
section (not necessary a bifurcation surface) is considered as a surface on which 
black hole entropy is calculated. Hence, the above argument arises no difficulties 
unless the cross section is taken to be the bifurcation surface. The derivation in 
Ref. ||] suffers from the above argument if and only if black hole entropy is estimated 
on the bifurcation surface. 

On the other hand, arguments like the above do not lead to any contradiction if 
we adopt a gauge condition such that £ a is fixed in a neighborhood of the bifurcation 
surface under variations, where £ a = £ a /k is a horizon generator Killing field with 
unit surface gravity. Moreover, it is concluded that, if we intend to fix a horizon 
generator Killing field, then it must have the same value of surface gravity for 
unperturbed and perturbed black holes. Hence, the gauge condition S^ a — in a 
neighborhood of the bifurcation surface adopted in Ref. |33|, fl8) is very natural one. 

In fact, it is always possible to identify unperturbed and perturbed stationary 
black holes so that the Killing horizons and the generator Killing fields with unit 
surface gravity coincide. As stated in Ref. |33| , such an identification can be done 
at least in a neighborhood of the horizon by using the general formula for Kruskal- 
type coordinates (U, V) given in Ref. Jl3| . (The identified Killing horizon is given 
by U — and V = 0. The identified Killing field with unit surface gravity is given 
by |° = U(d/dU) a - V{d/dV) a .) 

The purpose of the next subsection is to discuss the remaining gauge condition 
5t a = Sip a — at infinity. It is evident that this gauge condition at infinity can 
be imposed by identifying the perturbed and unperturbed spacetimes suitably. So, 
our question now is whether this gauge condition is necessary or not. For this 
purpose we temporarily adopt a gauge condition such that £ a is fixed everywhere 
on a hypersurface connecting the bifurcation surface and spatial infinity. In deriving 
the first law in this gauge condition, the gauge condition St a = S<p a = at infinity 
is found to be necessary for a prop er interpretation of the first law. On the other 
hand, as shown in subsection (2.1.3 , it is not necessary to fix £ a in a neighborhood 



of the bifurcation surface, strictly speaking. Hence, it can be concluded that the 
minimal set of gauge conditions necessary for the derivation of the first law is that 
t a and cp a are fixed at spatial infinity. 

2.1.2 The first law for stationary black holes 

Before deriving the first law, we review basic ingredients of the formalism. 

We consider a classical theory of gravity in n-dimensions with arbitrary matter 
fields, which is described by a diffcomorphism invariant Lagrangian n-form L(</>), 
where (f> denotes dynamical fields [Q . 

The Noether current (n — l)-form j[V] for a vector field V a is defined by 

}[V] = ®(<l>,£v<P)-V-L(cl>), (2.7) 
where the (n — l)-form &(4>, <5</>) is defined by 

SL{<f>) = E(cj))6<p + d®(<p,5<p). (2.8) 

It is easily shown that the Noether current is closed as 

dj[V] = --E(4>)C V 4> = 0, (2.9) 

where we have used the equations of motion E(</>) = 0. Hence, by using the machin- 
ery developed in Ref. |5C[ ], we obtain the Noether charge (n — 2)-form Q[V] such 
that 

j[V] = dQ[V]. (2.10) 
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Hereafter, we assume that in an asymptotically flat spacetime there exists an 
(n — l)-form B such that 

r V-SBtf) = [ V-®{cj>,84>), (2.11) 

J oo 

where the integral is taken over an (n— 2)-dimcnsional sphere at infinity. By using B, 
we can write a Hamiltonian H[V] corresponding to evolution by V a as follows |33|. 



H[V] = { (Q[V]-V-B). (2.12) 

J OO 



The symplectic current density uj((j), <5i0, 624*) is defined by 

u>((l>,5i(l>,52<t>) = 6i[&((t>,6 2 <j>)] - S 2 [&(4>,Si4>)] (2.13) 

and is linear both in 8\<j> and its derivatives, and 8 2 § and its derivatives |pl[ . 

Now we define a space of solutions in which we take a variation to derive the 
first law. 

Let £° be a fixed vector field, which vanishes on a (n— 2)-surface S. (Note that £ a 
and £ can be defined without referring to any dynamical fields, eg. the metric g a b-) 
In the following arguments, we consider a space of stationary, asymptotically flat 
solutions of the equations of motion E(0) = 0, each of which satisfies the following 
three conditions, (a) There exists a bifurcating Killing horizon with the bifurcation 
surface E. (b) £ a is a generator Killing field of the Killing horizon, (c) Surface 
gravity corresponding to £ a is 1: 

i b w b i a = r (2.i4) 

on the Killing horizon. 

For each element in this space, there exist constants n and (fi = 1, 2, • • •) 
such that 

Kl^f + n^), (2.15) 

where t a is the stationary Killing field with unit norm at infinity, {f^} (a* — 

1, 2, • • •) is a family of axial Killing fields. Hence, k is surface gravity and are 
angular velocities of the horizon. 

Note that, by definition, the vector field £ a is fixed under a variation of dynamical 
fields. We express this explicitly by denoting the variation by 5: 

Si a = 0. (2.16) 

We now derive the first law of black hole mec han ics. 

First, by taking a variation of the definition (2_7) for and using ( 2.16 ) and 
(|2.8|), we obtain 

= uj(<p,UXi4>) + d{i-®(4>,H)). (2.17) 

Here we have used the equations of motion E(</>) = and the following identity for 
an arbitrary vector V a and an arbitrary differential form A to obtain the last line. 

C V A = V ■ dA + d{V ■ A). (2.18) 

Since u>((j),5(f>,jC,?<f>) is linear in C^4> and its derivatives, we obtain 

d{5Q[£}) = d(£-©{<i>M)) (2-19) 
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by using Cg(f> — and E q. ([2.1C| ). 

Next we ingrate Eq. ( 2.1 9| ) over an asymptotically flat spacelike hypersurface C, 
which is parallel to at infinity and the interior boundary of which is E. Since 

£ ct = on E, we obtain 

' ' Q[£]=6H[t]. (2.20) 



Finally we rewrite the r.h.s. and the l.h.s. of (|2.20D in a form useful to be 
estimated at infinity and the horizon, respectively. 

A relation among variations of k, t a and is obtained by substituting 

( pup in (gTTeJ). 



By using this relatio n and the fact that H[V] is linear in the vector field V, we can 
rewrite the r.h.s. of ( 2.20 ) as follows. 

_ _ i _ . o^) . 

6H[$ = -(6H[t]-H[6t}) + -e-(5H[<p. M }-H[5<p M ]) 

hi hi 

i fi M 

= -SnHft + ^-SnHbto], (2.22) 
where the variation Sqq is denned for linear functionals F[t] and G(^) [^(^)] so that 

S^Flt] = SF[t]-F[k], 
SooG^ip^] = <5G( M )[v? (M )] - G (fl) [Sip {tl) }. (2.23) 

This newly introduced variation corresponds to a variation at infinity such that t a 
and tp a are fixed: 

<W a = <W( M ) = 0. (2.24) 
In Ref. [Est a useful expression of the Noether charge was given as follows. 

Q[v] = w c (0)y c + x crf (^)v [c y d] + Y(<f>, c v <t>) + dz(& v), (2.25) 

where W c , X cd , Y and Z are locally constructed covariant quantities. In particular, 
Y(0, Cv4>) is linear in Cv<j> and its derivatives, and X cd is given by 

(X cd (0)) C3 ... Cji =-^e afcC3 ... c „. (2.26) 

Here E^ cd is the would-be equations of motion form (48) for the Riemann tensor 
Rabcd and e Q &c 3 ■■■£:„ is the volume ro-form. 



By using the form of Q we can rewrite the integral in the l.h.s. of (2.2C) as 

/ Q[£] - f X rf (0)V [c efl, (2.27) 

where we have used the Killing equation C^4> — and the fact that £ a = on E. 
Using the relation 

V c | rf = e cd (2.28) 



on E, for any stationary solutions we can eliminate explicit dependence of Eq. ( 2.27 ) 
on £, where e c d is the binormal to E. Hence, at least within the space of stationary 
solutions, we can take the variation 5 of the integral without any difficulties. 
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Thus, we obtain the first law for stationary black holes by rewriting Eq. (2.20) 



as 



where entropy S is defined by 

S = 2^X cd (0)e cd , 
and energy £ and angular momenta J^) are defined by 



(2.29) 



(2.30) 



£ = H[t] = f (Q[t] — t ■ B), 

J OO 

J oo 



(2.31) 



Note that, in the r.h.s. of Eq. (2.29), variations of £ and are taken with t a 



and fixed. This condition is explicitly implemented by the definition (2.22) of 
<5oo and is necessary for a proper interpretation of the first law. 

We conclude this subsection by giving another expression of the entropy. 

Since £ a is a generator Killing field of the Killing horizon, we have Cgcf) = and 



the pull-back of £ • L(0) to the horizon vanishes. Hence, the definition (2.7) says 
that the pull-back of to the horizon is zero [Q. Thus, the integral of Q[£] is 
independent of the choice of the horizon cross section. 



Moreover, it can be shown that the integral in ( 2.30 ) does not change even if 
we replace the integration surface E by an arbitrary horizon cross section E' f49| . 
Therefore we obtain 

S = 2n f X cd (0K„ (2.32) 
where e' cd denotes the binormal to E'. 



2.1.3 Non-stationary perturbation 

In this subsection, we shall derive the first law for a non-stationary perturbation 
about a stationary black hole with a bifurcating Killing horizon. Unfortunately, for 
non-stationary perturbations, 5k and 5fl^ do not have meaning of perturbations of 
surface gravity and angular velocity of the Killing horizon, even if they are defined. 
However, since the first law ( |2 . 29| ) does not refer to 5k and Sfl^ but only to the 

unperturbed values of k and , we expect that the first law holds also for non- 
stationary perturbations. In the following, we shall show that it does hold. 

First, we specify a space of solutions in which we take a variation. 

Let £q be a fixed vector field, which vanishes on an fixed (n — 2)-surface E. 
In this subsection, we consider a space of asymptotically flat solutions of the field 
equation E(0) = 0, for each of which £q is an asymptotic Killing field. 

For each solution in this space, there exist constants k and (/x = 1, 2, • • •) 
such that at spatial infinity 

Kl^f+fi&Vfo, ( 2 - 33 ) 

where t a is a timelike asymptotic Killing field with unit norm at infinity, {y?^)} 
(fi = 1,2, •••) is a family of axial asymptotic Killing fields orthogonal to t a at 
infinity and {fl^} is a family of constants. Note that the constants k and 
do not have meaning of surface gravity and angular velocities unless we consider 
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a stationary solution. Moreover, in general, £g and £ have no meaning but an 
asymptotic Killing field and a fixed (n — 2)-surface, respectively. 

Note that, by definition, the vector field £q is fixed under the variation. We 
denote the variation by <5: 

= 0. (2.34) 

On the contrary, t a , ^f^y k and are not fixed under the variation since def- 
initions of them refer to dynamical fields, which are varied. Their variations are 



related by ( |2.2l| ) 



Suppose that an element (fro of the space of solutions satisfies the following three 
conditions, (a') <f>Q is a stationary solution with a bifurcating Killing horizon with 
the bifurcation surface E. (b') £q is a generator Killing field of the Killing horizon 
of <fro- (c') Surface gravity of <fr corresponding to £q is 1: 

l§V 6 |g = | a (2.35) 

on the Killing horizon. 

Now we derive the first law for the non- stationary perturbation Scfr about the 
stationary solution (fro. 

First, we mention that the validity of Eq. ( |2.20 ) in the previous section depends 



on the following three facts, (i) The equations of motion E(</>) = hold for both 
unperturbed and perturbed fields. (Unless they hold also for perturbed fields, 5j 
can not be rewritten as d(SQ).) (ii) £ a (corresponding to £q) is a Killing field of the 
unperturbed solution, (iii) £ a = (corresponding to £q — 0) on E for unperturbed 
solution. 

These three are satisfied for the unperturbed solution (fr and the non-stationary 
variation 5(fr around <po, too. Thus, Eq. ( [2.20 ) is valid, provided that £ a is replaced 

by is- 



Since 6t a , ^V 3 ?,,)! $ K an d 5fl^ are related by Eq. ( 2.21 ), we can transform the 



r.h.s. of ( |2.20D to obtain 



kS f QKo] = SvoS - ntfSnJw, (2.36) 

whe re, a s in the previous section, energy £ and angular momenta 3^) are defined 
by ( 2.3l| ), and the variation 8^ is defined at infinity so that t a and ip?^ are fixed. 



Here note that k and are surface gravity and angular velocities, respectively, 
for (fro . 

Up to this point we have not yet used explicitly the fact that ^ = on E for 
the perturbed solution, although we have used it implicitly. By using it explicitly, 



we can rewrite the l.h.s. of (2.36) in a useful form. The result is 



8 I Q[lo] = ^SS, (2.37) 



where S is defined by (2.30). (For explicit manipulations, see the proof of Theorem 
6.1 of Iyer-Wald @.) 



Finally, we obtain the first law (2.29) for non-stationary perturbations 5(fr about 
a stationary black hole solution (fro- 

Now we comment on entropy for the perturbed, non-stationary black hole. 

As stated above, the (n— 2)-surface S has no meaning for the perturbed solution. 
(It is nothing but a surface on which £q vanishes.) In general, it does not lie on the 
event (or apparent) horizon for the perturbed solution. Hence, entropy evaluated on 
£ may not coincide with that on a cross section of the perturbed horizon, provided 
that the entropy is defined as 2tt times an integral of Q [£ ] for both (n — 2)-surfaces. 
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Note that it is in general impossible to make gauge transformation so that E lie on 
a horizon cross section, if entropy (eg. a quarter of area in general relativity) on E 
is different from entropy on a horizon cross section. The difference is given by 2ir 
times an integral of the Noether current j [£o] over a hypersurface whose boundary 
is a union of E and a cross section of the perturbed horizon. Since it is natural 
to assign black hole entropy to the horizon cross section j4^] , it might be expected 
that there appears an extra term corresponding to the integral of j[£o] m the first 
law. 

However, as shown in the next paragraph, the integral of j[£o] vanishes to first 
order in 5(f> []. Thus, 5S evaluated on E gives the correct variation of entropy 
defined on the horizon to first order in 5<j). This means that the extra term does not 
appear and that the first law of Ref. derived in this section for non-stationary 
perturbation about a stationary black hole is the correct formula. 

Let us show the above statement. Since £q = on E and £e Q 4>Q = 0> the Noether 
current j[£o] vanishes on E for the unperturbed solution by the definition (2.7). 
Hence, for the perturbed solution, the Noether current is at least first order in 5<f) 
on E. On the other hand, deviation of a horizon cross section from E is at least 
first order. Therefore, the integral of j[£o] over a hypersurface connecting E and the 
perturbed horizon cross section is at least second order in 8<f>. 

Finally, let us apply the first law of this subsection to a stationary perturbation. 
The result is the same as that derived in the previous subsection, ft is evident that 
the gauge condition used in this subsection is weaker than that used in the previous 
subsection. In fact, £ a (^ £q for a perturbed solution) is not fixed in the former 
condition. Hence, it can be concluded that the minimal set of gauge conditions 
necessary for the derivation of the first law is that t a and (p9*. are fixed at spatial 
infinity. 



2.2 The first law of black hole dynamics 

As will be seen in section |2.3| , the first law of black hole statics is used in a (quasi- 
stationary but) dynamical situation to prove the generalized second law |5^, [Tc|| , 
which is a natural generalization of both the second law (or area law [jjjj ) of black 
hole and the second law of usual thermodynamics. In the proof, by assuming quasi- 
stationarity, the use of the first law of statics can be justified to relate a small 
change from an initial stationary black hole to final stationary one. However, if we 
intend to extend the proof of the generalized second law to finite changes between 
two stationary black holes or a purely dynamical situation, the first law of black 
hole statics can not be used. 

So, we want to extend the first law to a dynamical situation and call it a first 



law of black hole dynamics (BHD). It will be discussed in subsection 2.3.3 that 
the generalized second law might be extended to not quasi-stationary situations by 
using the first law of black hole dynamics. 

In this section, for simplicity, we consider general relativity only. In subsec- 
tion 2.2. 1| we consider two non-statistical definitions of entropy for dynamic (non- 
stationary) black holes in spherical symmetry. The first is analogous to the original 
Clausius definition of thermodynamic entropy: there is a first law containing an 
energy-supply term which equals surface gravity times a total differential. The sec- 
ond is Wald's Noether-charge method, adapted to dynamic black holes by using 
the Kodama flow. Both definitions give the same answer for Einstein gravity: one- 



quarter the area of the trapping horizon fl41|. In subsection 2.2.2, the first law of 



BHD is derived without assuming any symmetry and any asymptotic conditions 

The author thanks Professor R. M. Wald for helpful comments on this point. 
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In the derivation, a definition of dynamical surface gravity is proposed. 

2.2.1 Dynamic black hole entropy in general relativity with 
spherical symmetry 

It is generally thought that black-hole entropy should have a statistical origin, pre- 
sumably in a quantum theory of gravity. This is, of course, due to the definition 
of entropy in statistical mechanics. However, it should be remembered that the 
original concept of entropy was not statistical p3| . The original argument of Clau- 
sius was that, in a cyclic reversible process, the total heat supply SQ divided by 
temperature z9 should vanish. Thus in any reversible process, 8Q/$ should be the 
total differential dS of a state function S, the entropy. Moreover, in irreversible 
processes, there should be a second law dS > SQ/i). The heat supply also occurs 
in a first law dU = 5Q + SW, where U is the internal energy and SW the work 
being done. These are basic laws of thermodynamics as stated in typical textbooks 
and originally formulated by Clausius before the invention of statistical mechanics. 
In this subsection, we argue that there is a similar concept of entropy for dynamic 
black holes, suggested by the mathematical structure of the first law: it contains an 
energy-supply term which equals surface gravity times a total differential. 

Kodama vector and the first law of black hole dynamics in spherical 
symmetry 

The relevant quantities and equations in spherical symmetry may be summarized 
as follows. The area A or areal radius r = y/A/Aii of the spheres of symmetry 
determines the 1-form 

k = *dr (2.38) 

where d is the exterior derivative and * is the Hodge operator of the two-dimensional 
space normal to the spheres of symmetry. Henceforth, 1-forms and their vector 
duals with respect to the space-time metric will not be distinguished. Then k is the 
divergence- free vector introduced by Kodama |54| , which generates a preferred flow 
of time and is a dynamic analogue of a stationary Killing vector | p5| . The active 
gravitational energy or mass is 

E = (1 - dr ■ dr)r/2 (2.39) 

where the dot denotes contraction. Misner and Sharp originally defined E and, 
thus, is called the Misner-Sharp energy. (Ref. [j57j described its physical properties.) 
The dynamic surface gravity 

n = *dk/2 (2.40) 

was defined in Ref. |55| ] by analogy with the standard definition of stationary surface 
gravity. This reference also introduced two invariants of the energy tensor T: the 
energy density (work density) 

w = -trT/2 (2.41) 
and the energy flux (localized Bondi flux) 

i/j = T-dr + wdr (2.42) 

where tr denotes the the two-dimensional normal trace. One may say that (A, k) 
are the basic kinematic quantities, (E, k) the gravitational quantities and (w, ip) 
the relevant matter quantities. Instead of ip one may also use the divergence-free 



energy-momentum vector 55, 54, p7\ 



j = + wk. 



(2.43) 
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Finally, the relevant components of the Einstein equation are [ p5[ 

E = r 2 n + ATrr 3 w (2.44) 

and |7) 

Aj = *dE. (2.45) 

The latter may be rewritten as 

dE = A1P + wdV (2.46) 

where V = ^7rr 3 is the areal volume. This is the unified first law of Ref.j55|]. One 
may regard wdV as a type of work and Aip as an energy supply, analogous to heat 
supply 5Q — § q, where q is the heat flux. The energy supply can be written as 

ndA , f E\ 

Aip = h rd — . (2.47) 

8n \ r ) 

The second term vanishes when projected along a black-hole horizon, defined as in 
Refs.|35|, [58| , |57f by a trapping horizon: a hypersurface where dr is null, so that 
E = r/2. This also occurs for any hypersurface on which E/r is constant, thereby 
covering any smooth space-time. The key point is that the first term is the product 
of surface gravity n and a total differential. Identifying k/2tt as a temperature, this 
total differential therefore determines a Clausius entropy A/ A. Note that this stems 
from a purely mathematical property of the energy supply occurring in the first 



law. The restriction to spherical symmetry will be removed in subsection 2.2.2 



Wald-Kodama entropy 

Wald |3^] also gave a definition of entropy 

kS = 2tt <j> Q[£] (2.48) 

where Q henceforth denotes a Noether charge 2-form obtained from a certain type 
of Lagrangian. For Einstein gravity, Q a {, = — e a fc C( iV c £ d /167r p8[ , where e is the 
space-time volume form, £ is a generating vector for the diffeomorphisms which is 
taken to be the horizon generating Killing vector of a stationary black hole, and 
k is the surface gravity corresponding to £. In spherical symmetry, we propose 
using the Kodama vector k for £ to give an alternative definition of the entropy of 
dynamic black holes. We call it Wald-Kodama entropy. This prescription effectively 
corresponds to also replacing £ by k in S2 of Jacobsen et al. [Q. Then from the 
above expression for Q by Wald's method, 



j>Q[k] = ^, (2.49) 



where the integral is over a sphere of symmetry. Thus the Wald-Kodama entropy 
is 

S = A/A. (2.50) 



Agreement and speculations 

So, the Wald-Kodama entropy agrees with the Clausius entropy. The motivations 
also seem similar, since Wald's construction involved a first law of black-hole statics 
based on perturbations S of a stationary solution. This agreement suggests that 
some combination of the two methods may be useful in general, assuming neither 
stationarity nor Einstein gravity. 
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It is also interesting that both methods formally hold not just on a black-hole 
horizon, but anywhere in the space-time, as do all the equations displayed above. 
Whether any surface in any space-time should have an entropy related to its area 
is arguable, but this does concur with the entanglement entropy approach, which is 



discussed in section 3.3 



2.2.2 Qusi-local first law of black hole dynamics in general 
relativity 

A general definition of a black hole 

Here we would like to treat a dynamical, not necessarily asymptotically flat space- 
time. Even for such a general situation, there is a definition of a black hole. Namely, 
a black hole is defined as a future outer trapping horizon |58| . The future outer trap- 
ping horizon is the closure of a three-surface foliated by marginal surfaces on which 
(9_ < and C-Q + < 0, where the marginal surface is a spatial two-surface on which 
one of two null expansions (which we have denoted by 9 + ) vanishes. Here (9_ is 
another null expansion and £_ is a Lie derivative w.r.t. a null vector defined below. 
For the purpose of this subsection, we only need the fact that 8 + 6- = along the 
horizon. Hence, the first law we shall obtain remains to hold for a general trapping 
horizon, i.e. a hypersurface foliated by marginal surfaces. We mention here that 
a trapping horizon can be regarded as a black hole, a white hole and a wormholc 
when it is future outer, past outer and temporal outer, respectively |5{l] . 



The double-null formalism 

To investigate behavior of the trapping horizon, the so-called double-null formalism, 
or the (2 + 2) decomposition, of general relativity is useful. Among several (2 + 2)- 
formalisms |30[ [H], we adopt one based on Lie derivatives w.r.t null vectors devel- 
oped by Hayward ]6l| ]. Let us review basic ingredients of the formalism. Suppose 
that a four-dimensional spacetime manifold (M, g) is foliated (at least locally) by 
two families of null hypersurfaces S ± , each of which is parameterized by a scalar 
respectively. The null character is described by (7 _1 (n ± , n^) — 0, where = —d^ 
are normal 1-forms to T, ± . The relative normalization of the null normals defines a 
function / as g~ 1 (n + ,n~) = — e* . The intersections of £ + (£ + ) and £~(£~) define 
a two-parameter family of two-dimensional spacelike surfaces S(£, + ,l;~). Hence, by 
introducing an intrinsic coordinate system (Q 1 ,0 2 ) of the 2-surfaces, the foliation is 
described by the imbedding x — 6 1 , 6 2 ). 

For the imbedding, the intrinsic metric on the 2-surfaces is found to be h = 
g + e~f(n + ® n~ + n~ ® n + ). Correspondingly, the vectors u± — d/d£ have 'shift 
vectors' s± =_L it±, where _L indicates projection by h. The 4-dimensional metric 
is written in terms of (h,f,s±) as 

/ h(s+,s+) h(s + ,s^)-e~ f h(s+) \ 
g=\ h(s^,s+)~e~ f /i(s_,s_) ft(s_) . (2.51) 
\ h(s + ) /i(s_) h J 

Geometrical quantities such as expansions 9± , shears o~± and the twist lo are defined 
by 6± — *C± * 1, cr± =-L C±h — 9±h and u = e' [£_, l+]/2, where * denotes the 
Hodge-dual operator of h, l± = u± — s± = e _ ^g _1 (n =F ) are null normal vectors to 
S ± , and C± denotes the Lie derivative along l±, respectively. It is possible to write 
down the Einstein tensor in terms of these geometrical quantities. The component 
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useful for our purpose is = G(Z+, 1-), which is given by 



2e / G+_ = (2) i? + e / (^ 



29+9 ) - 2 



V 2 f. 



(2.52) 

Here h$ = g^ x hg^ x is h raised by g^ 1 , V 2 and ^R are the two-dimensional Lapla- 
cian and the Ricci scalar both associated with the metric h. 



Hawking energy 

Before deriving the first law we have to define energy and surface gravity in a 
quasi-local way. In spherical symmetry there is a widely accepted energy: the 
Misner-Sharp (MS) energy (2.39) [Q. In the previous subsection the MS energy is 
used to derive the first law of BHD in spherical symmetry |S . In this subsection 
we adopt the Hawking energy p2] |, which reduces to the MS energy in spherical 
symmetry It is defined by 



16tt 



E{£, + ,C) = -^r d 2 9Vh\ (2) R + e f 9 + 9 



(2.53) 



where h is the determinant of the two-dimensional metric h a b and the area radius 
r is defined by 



r = y/A/Air, A = 

Js 



d 2 9Vh. 



(2.54) 



A proposal of dynamical surface gravity 

In Ref. Ipql , a definition of dynamical surface gravity was proposed in spherical 
symmetry as (2.40). A natural generalization to a non spherically-symmetric case 



-i 



167rr 



d 2 9Vh e f {£+9-+£-9 + + i 



(2.55) 



This is the most simple generalization in the sense that it includes neither the shear 
a Aab nor the twist uj a . 



Note that this definition of surface gravity and the definition (|2.53| ) of the Hawk- 
ing energy are both quasi-local in the sense originally introduced by Penrose |63|| . 
Hence we call the corresponding first law, which we shall derive below, the quasi- 
local first law of BHD. 



The dynamical first law 



We now derive the quasi-local first law of BHD for the Hawking energy ( p.53| ) and 
the surface gravity defined by ( 2.55 ). It is easy to show that 



dE 



where w is defined by 



-dA = wAdr 



1 E 



E 



w = — nr 

A\r 



(2.56) 



(2.57) 



Here note that 'cf in Eq. (2.56) is not a variation in a space of stationary solu- 
tions of the Einstein equation as in the first law of black hole statistics, but is the 
differentiation w.r.t. the parameters £ of the spacetime foliation. (For example, 
dE = d£ + d + E + d£~d-E.) We mention that Eq. (2.56) holds independently of the 
definitions of E and n while the following arguments depend on the definitions. 
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Since the Gauss-Bonnet theorem says that 

d 2 9Vh (2) i? = 87r(l- 7 ), 



where 7 is the genus or number of handles of S(£ + , £ - )> the energy divided by area 
radius is given by E/r = (1 — 7) /2 on a marginal surface and is a constant. Thus, 

E 1 = — A' + wAr', (2.58) 

87T 

where the prime denotes the derivative along the trapping horizon. This is the quasi- 
local first law of BHD. Note that this also holds along any hypersurface foliated by 
2-surfaces on which E/r is constant. 

The work term 



By using Eq. ( 2.52 ) it is easy to show that w is written as follows. 

w = p m + pj, (2.59) 

where the averaged matter energy density p m and the effective angular energy 
density pj are defined by 



'•>■■ -^r / d 2 eVhe f G+., 



8ttA 

Pj = 



8ttA 



d 2 0Vh 



h(cj,cj) + ^h i (df,df) 



(2.60) 



The Einstein equation G — 8ttT says that p m is e* T(l+, Z_) averaged over the 
2-surface. It seems that pj represents effective energy density due to angular mo- 
mentum. 

The term wAr' should be a work term done along the horizon. For example, 
for an electromagnetic field, the term p m Ar' reduces to the electromagnetic work 
done along the horizon j55j. It seems that the term PjAr' is a work associated with 
angular momentum of the trapping horizon. 

Comments 

In this subsection the quasi-local first law of black hole dynamics has been derived 
without assuming any symmetry and any asymptotic condition. In the derivation 
we have given a new definition of dynamical surface gravity. In spherical symmetry 
it reduces to that defined in Ref. [5^ ]. 

By using the quasi-local first law derived in this subsection, it might be possible 
to extend a proof of the generalized second law to not quasi-stationary situations. 



(See subsection |2. 3. 3| .) 

Besides the first law derived in this subsection, there exist the second law 
and the third law |I4[ ] for the trapping horizon (or apparent horizon) . It seems that 
by using these laws we can formulate black hole thermodynamics consistently as 
trapping horizon dynamics. However, for this purpose, there is an important open 
question: we have to associate temperature of quantum fields with the trapping 
horizon. All we can say here is that the temperature may be given by %k/2it, where 
k is the surface gravity introduced here. 

The final comment is in order. The surface gravity k(£ + ,£~) is an invariant of 
a double-null foliation at the surface. Since a non-null trapping horizon locally de- 
termines a unique double-null foliation, the surface gravity is definitely an invariant 
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of the trapping horizon if the horizon is not null. On the other hand, the null case 
is ambiguous because of the freedom to rescale the other null direction. Fixing this 
would require some kind of limiting argument that might be effectively a zcroth 
law. Therefore, we have to impose an auxiliary condition for the surface gravity 
/■t(£ + ,£ - ) to work well when the trapping horizon is null. Since surface gravity 
seems to be related to temperature of quantum fields as stated above, it will be 
valuable to investigate the auxiliary condition in detail p2fl . 



2.3 The generalized second law 

The generalized second law of black hole thermodynamics is a statement that a 
sum of black hole entropy and thermodynamic entropy of matter fields outside the 
horizon does not decrease 0, |(35|, [32[ where the black hole entropy is defined 
as a quarter of the area of the horizon. Namely it says that an entropy of the whole 
system does not decrease. It interests us in a quite physical sense since it links a 
world inside a black hole and our thermodynamic world. In particular it gives a 
physical meaning to black hole entropy indirectly since it concerns the sum of black 
hole entropy and ordinary thermodynamic entropy, and since physical meaning of 
the latter is well-known by statistical mechanics. 

Frolov and Page [ p2[ proved the generalized second law for a quasi-stationary 
eternal black hole by assuming that a state of matter fields on the past horizon is 
thermal one and that a set of radiation modes on the past horizon and a set of 
radiation modes on the past null infinity are quantum mechanically uncorrelated. 
The assumption is reasonable for the eternal case since a black hole emit a thermal 
radiation (the Hawking radiation). When we attempt to apply their proof to a 
non-eternal black hole which is formed by gravitational collapse, we might expect 
that things would go well by simply replacing the past horizon with a null surface 
at a moment of a formation of a horizon (v — vq surface in Figure [2.1| ). However, it 
is not the case since the above assumption does not hold in this case. The reason 
is that on a background describing gravitational collapse the thermal radiation is 
observed not at the moment of the horizon formation but at the future null infinity 
and that any modes on the future null infinity have correlation with modes on the 
past null infinity located after the horizon formation. The correlation can be seen 



in the equation ( 2.65 ) of this section explicitly. Thus, their proof does not hold for 
the case in which a black hole is formed by gravitational collapse. Since a black 
hole is thought to be formed by gravitational collapse in astrophysical situation, we 
want to prove the generalized second law in this case. 



The rest of this section is organized as follows. In subsection 2.3.1 we consider 
a real massless scalar field in a background of a gravitational collapse to show that 
a thermal state with special values of temperature and chemical potential evolves 
to a thermal state with the same temperature and the same chemical potential. 
These special values are determined by the background geometry. In subsection 



2.3.2, first, the generalized second law is rewritten as an inequality which states 
that there is a non-decreasing functional of a density matrix of matter fields. After 
that, we give a theorem which shows an inequality between functionals of density 
matrices. Finally, we apply it to the scalar field investigated in subsection 2.3.1 to 



prove the generalized second law for the quasi-stationary background. In subsection 



2.3.3 we summarize this section. 



2.3.1 A massless scalar field in black hole background 

In this subsection we consider a real massless scalar field in a curved background 
which describes formation of a quasi-stationary black hole. Let us denote the past 
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null infinity by T~ , the future null infinity by T + and the future event horizon by 
H + . Introduce the usual null coordinates u and v, and suppose that the formation 
of the event horizon H + is at v = vq (see Figure 2.1). On I~ and Z + , by virtue of 
the asymptotic flatness, there is a natural definition of Hilbert spaces Tlx- and Tlx+ 
of mode functions with positive frequencies |37| . The Hilbert spaces J-(Tlx± ) of all 
asymptotic states are defined as follows with a suitable completion (symmetric Fock 
spaces): 

T{Hx±) = C®Hx±@ (H x ± ® Hx±) sym © • • • , 

where (• • -) S y m denotes the symmetrization ((£ ® r/) sym = |(£ r\ + -q ® £), etc.). 
Physically, C denotes the vacuum state, Tlx± one particle states, (Tlx± <8> Tlx±) sym 
two particle states, etc.. We suppose that all our observables are operators on 
T(Tlx± ) since we observe a radiation of the scalar field radiated by the black hole 
at places far away from it. In this sense J r (Tlx+) are quite physical. Next let us 
consider how to set an initial state of the scalar field. We want to see a response of 
the scalar field on the quasi-stationary black hole background which is formed by 
gravitational collapse of other materials (a dust, a fluid, etc.). Hence as the initial 
state at X~ we consider a state such that it includes no excitations of modes located 
before the formation of the horizon (no excitation at v < vq). A space of all such 
states is a subspace of J 7 (Tlx-), and we denote it by Tx-( v ->v Q )- We like to derive a 
thermal property of a scattering process of the scalar field by the quasi-stationary 
black hole. Hence we consider density matrices on Tx-(v>v Q ) and T(Tlx+). Denote 
a space of all density matrices on Tx- ( v >v ) by V and a space of all density matrices 
on T(Hx+) by V. 

Let us discuss the evolution of a state at 1~ to future. Since T + is not a Cauchy 
surface because of the existence of H + , J- (Tlx-) is mapped not to T(7ix+) but 
to J-(Hx+) <8> J-(TIh+) by a unitary evolution, where TIh+ is a Hilbert space of 
mode functions on the horizon with a positive frequency, and J-(TIh+) is a Hilbert 
space of all states on H + defined as a symmetric Fock space (see the definition of 
3~(Hz±)). Although there is no natural principle to determine positive frequency 
modes (equivalently, there is no natural definition of the particle concept) on H + , 
how to define Hh+ docs not affect the result since we shall trace out the degrees of 
freedom of !F(TIh+ ) (see (2.61)). To describe the evolution of a quantum state of the 
scalar field from J- (Tlx- ) to T(TIx+)®J : (TIh+ ) an S-matrix is introduced |Q. For a 
given initial state \tp) in J 7 (Tlx-), the corresponding final state in F(TIx+)®J :: (TIh+) 
is S\ip) . Then the corresponding evolution from Tx-{v>v Q ) to T(Tlx+) is obtained 
by restricting S to J~x-( v >v n )i an d we denote it by S, too. The S-matrix elements 
was given by Wald |67j . 



Superscattering matrix T 

Suppose that the initial state of the scalar field is \<j)} (G F(Tlx-)) and that the 
corresponding final state is observed at T + (see the argument after the definition 
of T(Tlx±)). Formally the observation corresponds to a calculation of a matrix 
element (4>\S^ OS\<i>) , where S is the S-matrix which describes the evolution of the 
scalar field from Fx-(v>vo) to ? r (Tlx+) ® T(Tl H +) and O is a self-adjoint operator 
on J r (Tlx+) corresponding to a quantity we want to observe. The matrix element 
can be rewritten in the following convenient fashion: 

(<f>\S^OS\<f>) = Tr I+ [Op red ], 

where 

Pred = Tr H+ [S\ct>){4>\&] , 

TYj+ and Tr H + denote partial trace over T(Hx+) and T(Tl H +), respectively. In 
viewing this expression we are lead to an interpretation that the corresponding final 



CHAPTER 2. LAWS OF BLACK HOLE THERMODYNAMICS 



22 



state at X + is represented by the reduced density matrix p re d- Next we generalize 
this argument to a wider class of initial states, which includes all mixed states. For 
this case an initial state is represented not by an element of ^Fx~(v>v ) but by an 
element of V (a density matrix on Tx- („>„ )). Its evolution to T + is represented by 
the so-called superscattering matrix T defined as follows: let p (e V) be an initial 
density matrix then the corresponding final density matrix T(p) (£ V) is 

T(p) = Tr H+ [SpS^] . (2.61) 
Note that T is a linear map from V into V . 



Thermodynamic property of T 

Let us calculate a conditional probability defined as follows: 

P({n iP }|{n i7 })^({n sP }|T(|{n i7 })({n s7 }|)|{n sP }), 

where 



(2.62) 



IR?» 
IK.}) 



n^H^H 1°)' 

n i (atfc,))V j.n. 



(2.63) 



|{n i7 }) is a state in Fx-(v>vo) characterized by a set of integers n i7 (i = 1, 2, • ■ •) 
and |{%,}) is a state in J-(Ttx+) characterized by a set of integers n iP (i = 1, 2, • • •). 
Therefore ^{"splK"^}) i s a conditional probability for the final state to be |{n iP }) 
when the initial state is specified to be |{ji <7 }). In the expressions, A is a part 
of a representation of a Bogoliubov transformation ^7j, which represents a map 
from Hj+ © T~Ih+ to Ti-x- > an d »7 is a unit vector in Hj+ © 7i_f/+ such that A ,7 
corresponds to a wave packet whose peak is located at a point on T~ later than the 
formation of the horizon (y > Vo) |6^| . On the other hand, jp i s a unit vector in 
?i x+ and corresponds to a wave packet on T + |6^] (see Figure 2A ) . The probability 
( 2 .62) ) is a generalization of P{k\j) investigated by Panangaden and Wald ||. (Our 
P({n iP }|{n i7 }) reduces to P(k\j) of Panangaden- Wald when n il =j,n iP — k 
and n i7 = n iP = for all i other than iq. Here iq is an arbitrary fixed value of 
i.) Evidently, our conditional probability P({w i p}|{w i 7}) includes more abundant 
information 2 about a response of the scalar field than P(k\j). In fact, any initial 
states on T~ , which include no excitation before the formation of the horizon (v < 
vq.), can be represented by using the basis {|{n. i7 })} and any final states on T + can 



be expressed by the basis {|{n iP })}, i.e. a set of all |{n i7 }) generates T^Hx- 



(v>V 



)) 



and a set of all |{n iP }) generates J 7 (Hx+)- This is the very reason why we have 
generalized P(k\j) to P({n iP }\{n n }). 

By using the S-matrix elements given in |67|, the conditional probability is 



rewritten as follows (see appendix A.l for its derivation): 



P(M\{n a }) 

i 

min(n. 7 ,n. p ) min(n. 7 ,n. p ) 

x E ' E ' 



[-\Ri\ 2 /(l-\Ri\ 2 )] 



"1 n. .. t 



'it iP' 



U=0 



h\(n n - h)l{n iP - li)\mi\(n n - mi)!(ri iP - m»)! 



{ n -7}{ n/ T } 

All the information about the response of the scalar field is included in T, , . ? . defined 
in Lemma bl 
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Figure 2.1: A conformal diagram of a background geometry which describes a 
gravitational collapse. T~ and I + are the past null infinity and the future null 
infinity, respectively and H + is the future event horizon. Shaded region represents 
collapsing materials which forms the black hole. Besides the collapsing matter, we 
consider a real massless scalar field and investigate a scattering problem by the 
black hole after its formation (v > v n ). Thus we specify possible initial states at I~ 
to those states which are excited from the vacuum by only modes whose support 
is within v > v (elements of J 7 I -( V>V0 )), and possible mixed states constructed 
from them (elements of V). In the diagram, A ,7 (i = 1, 2, • • •) is a mode function 
corresponding to a wave packet whose peak is at v > vq on Z~ , ip (i = 1, 2, • • •) is a 
mode function corresponding to a wave packet on I + . 
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E 



rii—n. p —m\n{li ^m^) 



(m - n iP + li)\{n t - n iP + rm) 



(x 2 i\Ri\ 



2\rii —n. , 



, (2-64) 



where Ri is a reflection coefficient for the mode specified by the integer i on the 
Schwarzschild metric (see appendix A..1) and constant defined by Xi = 

exp(— 7r(wi — Osjjmi)/fc). fn the expression, Ui and m, are the frequency and the 
azimuthal angular momentum quantum number of the mode specified by the integer 
i, SIbh and k are the angular velocity and the surface gravity of the black hole. 

Now, the expression in the squared bracket in ( [2.64 ) appears also in the calcu- 
lation of P(k\j). Using the result of M, it is easily shown that 



P({n iP }\{n a }) = H 



min(n.p,n. T ) 



(»,/ 



Si\{n iP - Sj)!(n 47 - S;)! 



(2.65) 



where 



Ki 



(1 



\R. 



\2\ n ii+ n iP 



(f- \R l \ 2 x 2 J^ +n ' p+1 

(|i? 4 | 2 -^) (f-|^l 2 *?) 



This is a generalization of the result of 
derived by using this expression. 



(l-l^l 2 ) 2 *? 

and the following lemma is easily 



Lemma 1 For the conditional probability defined by (2.6i) the following equality 
holds: 



P{{n iP = h}\{n il =k})e- pBH ^ k{ ' 
- P{{n iP =Ji}|{n i7 - ki })e- t}n "^i k ^- n "" mi \ 



(2.66) 



where u>i and mi are the frequency and the azimuthal angular momentum quantum 
number of the mode specified by i, £Ibh is the angular velocity of the horizon and 

Pbh = 2tt/k. 
Here k is the surface gravity of the black hole. 

Note that /3g# is the Hawking temperature of the black hole. This lemma 
states that a detailed balance condition holds Summing up about all fc's, we 
expect that a thermal density matrix pth(PBHi ^bh) in V with a temperature f3g\j 
and a chemical potential Qbh for azimuthal angular momentum quantum number 
will be mapped by T to a thermal density matrix p~th{PBH i Qbh) in V with the 
same temperature and the same chemical potential. To show that this expectation 
is true, we have to prove that all off-diagonal elements of T {pth( §BH > ^bh)) are 



zero. For this purpose the following lemma is proved in Appendix A. 2 
Lemma 2 Denote a matrix element of T as 

T™$$ ee ({n iP }|T(|{n <7 })({< 7 }|) 



(2.67) 



s It guarantees that a thermal distribution of any temperature is mapped to a thermal dis- 
tribution of some other temperature closer to the Hawking temperature, as far as the diagonal 
elements are concerned. 
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Then 

T Uk} - °> (2 - 68) 

unless 

n n — n' a = n iP — n' p (2.69) 

for Vi . 

Lemma || shows that all off-diagonal elements of T(p) in the basis {|{?\p}}} 
vanish if all off-diagonal elements of p in the basis is zero. Thus, combining 

it with Lemma |l| and the well-known fact that |0) is mapped to the thermal state, 
the following theorem is easily proved. Note that a set of all |{7Vy})({^ 7 }| generates 



V and a set of all |{n- iP })({n i ' p }| generates V (see the argument below ( |2.63| )) 



Theorem 3 Consider the linear map T defined by (2.61) for a real, massless scalar 
field on a background geometry which describes a formation of a quasi- stationary 
black hole. Then 

T (pth(/3BH,^BH)) = Pth(pBH, n BH ), (2.70) 

where 

PthifiBH^BH) = Z- 1 £ e-ft^E,^"*- nm) \{n n })({n n }\, 

i n i-i} 

PthWsH^BH) = z-' J2 ^ BH ^ nA ^ BHm ' ] \M){Ml 

i n ip} 

% = e~^ BH Ei i.( u .- n sH m i) (2 71) 

bi} 

Pthif^BH, &bh) and pthif^BH, ^bh) can be regarded as 'grand canonical ensem- 
ble' in V and V respectively, which have a common temperature P BH and a common 
chemical potential Qbh for azimuthal angular momentum cjuantum number. Thus 
the theorem says that the 'grand canonical ensemble' at X" (v > vo) with special 
values of temperature and chemical potential evolves to a 'grand canonical ensemble' 
at T + with the same temperature and the same chemical potential. Note that the 
special values (3g\j and SIbh are determined by the background geometry: /3g H is 
the Hawking temperature and fl BB is the angular velocity of the black hole formed. 



This result is used in subsection 2.3.2 to prove the generalized second law for the 
quasi-stationary black hole. 



2.3.2 A proof of the generalized second law 

The generalized second law of black hole thermodynamics is 

ASbh + AS matter > 0, (2.72) 

where A denotes a change of quantities under an evolution of the system, Sbh and 
Smatter are black hole entropy of the black hole and thermodynamic entropy of the 
matter fields, respectively. For a quasi-stationary black hole, using the first law of 
the black hole thermodynamics 

ASbh = (3 BH (AM BH - n BH AJ BH ), 

the conservation of total energy 

AM BB + AE mat ter = 
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and the conservation of total angular momentum 

A JbH + ALmatter = 0, 

it is easily shown that the generalized second law is equivalent to the following 
inequality: 

matter — ^BHAL matter ) > 0, (2.73) 

where (3b h, Qbh, Mbh and Jbh are the inverse temperature, the angular velocity, 
the mass and the angular momentum of the black hole; E mat ter and L matte r are the 
energy and the azimuthal component of angular momentum of the matter fields. 
Thus this is of the form 

U[po;(3 B hMbh) > U[p ; Pbh&bh), (2.74) 

where U is a functional of a density matrix of the matter fields defined by 

U[ P -Pbh^bh) = Tr \p]np] - f3 BH (Tr[Ep] - n BH Tr[L z p\) , (2.75) 

Po and po are an initial density matrix and the corresponding final density matrix 
respectively. In the expression E and L z are operators corresponding to the energy 
and the azimuthal component of the angular momentum. Note that ( [2.74 ) is an 



inequality between functionals of a density matrix of matter fields ^. We will prove 
the generalized second law by showing that this inequality holds. Actually we do 
it for a quasi-stationary black hole which is formed by gravitational collapse, using 



the results of subsection 2.3.1 and a theorem given in the following. 



Non-decreasing functional 

In this subsection a theorem which makes it possible to construct a functional 
which does not decrease by a physical evolution. It is a generalization of a result 



of |66j. After that, we derive ( 2.74 ) for a quasi-stationary black hole which arises 



from gravitational collapse, applying the theorem to the scalar field investigated in 



subsection 2.3.1 



Let us consider Hilbert spaces T and T . First we give some definitions needed 
for the theorem. 

Definition 4 A linear bounded operator p on T is called a density matrix, if it is 
self-adjoint, positive semi-definite and satisfies 

Trp = 1. 

In the rest of this section we denote a space of all density matrices on F as V{T). 
Evidently V{T} is a linear convex set rather than a linear set. 

Definition 5 A map T ofV(J-) into V{T} is called linear, if 

T [api + (1 - a)p 2 ) = aT{ Pl ) + (1 - a)T(p 2 ) 

/or0< v a<l and y Pl , y p2 

By this definition it is easily proved by induction that 

N 

= J2*iT(pi), (2-76) 

i=l 

if a% > 0, Sili o» = 1 and pi e V{T). 

Now we prove the following lemma which concerns the N — > oo limit of the left 
hand side of (f2j|). We use this lemma in the proof of theorem [?]. 




4 Information about the background geometry appears in the inequality as the variables $bh 
and Qbh which parameterize the functional. 
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Lemma 6 Consider a linear map T of V(T) into V{T) and an element p of 
V{T). For a diagonal decomposition 

oo 

Po = ^2pi\i)(i\, 

i=l 

define a series of density matrices of the form 

n 

p n = Y,Pi/an\i)(i\ (n = N,N+l,---), (2.77) 

i=l 

where 

n 

a n = ^2pi 

»=i 

and N is large enough that a at > 0. Then 

lim ($|T(^)|*) = ($|T(po)|*) (2.78) 

n^oo 

/or arbitrary elements |<I>) anrf |^) o/ J 7 . 

This lemma says that T{p n ) has a weak-operator-topology-limit T(p ). 
Proof 

By definition, 

po = anpn + (1 - a n )p' n , (2.79) 

where 

J = f E,~n+lP*/( 1 -««)|j)(i| K < 1) 
^" 1 K = 1) ' 

Then the linearity of T shows 

(#|T(a,)|*) = a„($|T(p„)|*> + (1 - a„)($|T(^)|v&). 
Thus, if (<&|T(p^)|\I>) is finite in n — > oo limit, then the lemma is established since 

lim a n = 1. 

71— >00 

For the purpose of proving the finiteness of ($|T(p^)|'I'), it is sufficient to show that 
|(4>|/5|\I')| is bounded from above by ||$|| ||^|| for an arbitrary element p of V{T). 
This is easy to prove as follows. 

\mm)\ = i e«wx ? i*>i ^ E i wx^>i < 11*11. ( 2 - 8 °) 

where we have used a diagonal decomposition 

p = ^ j Pi\i){i\. 



□ 

Theorem 7 Assume the following three assumptions: a. T is a linear map of 
V{T) into V{T), b. f is a continuous function convex to below and there are 
non-negative constants C\, C2 and C3 such that |/((1 — e)x) — f(x)\ < |e|(ci|/(a;)| + 
C2\x\ + C3) for y x (> 0) and sufficiently small \e\, c. there are positive definite 
density matrices p^ (e and p^ (g V{T)) such that T{p 00 ) = p^. 
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If [poo, po] = [Poc,T(po)} = and Tr[p oo \f{p p o ^)\] < oo, then 

U [T(p )] > U [p ] , (2.81) 

where 

U[p] = Tr [ Poo f (pp^ 1 )] , 

U[p] = -Tr [ Poo f (pp^ 1 )] . (2.82) 

As stated in the first paragraph of this subsection, theorem |?] is used in subsection 
[| to prove the generalized second law for a quasi-stationary black hole which arises 
from gravitational collapse. 
Proof 

First let us decompose the density matrices diagonally as follows: 

oo oo 

Pa = ^2p%\i)(i\, Poo = ^2qi\i)(i\, 

oo oo 

T(p ) = X>|?><i|, T (P~) = £gi|i><i|. (2.83) 



Then by lemma |5| and ( 2.76 ) , 

n 

pi = (i\T(p )\i) = lim Viyft/fl,,, (2.84) 

n— >oo * — ' 

3=1 

where a n = X)£=iPi an d A^ = (i\T(\j)(j\)\i). has the following properties: 

oo 

j2 A ij = i , o<Ay<i. 



1=1 

Similarly it is shown that 



qi = lim V" Aijqj/br, 



3=1 



where b n = Y17=i 1i- By (2.84) and the continuity of /, it is shown that 

/&/&) = lim / (f>^) . (2.85) 
Next define Cf and <5f by 

n 

Cf=^^/%: C«^Cr/a„, (2.86) 
j'=i 



then the convex property of / means 



3=1 

since 



E ^'jgj _ i ^3<?3 > q 
j=l i qi i * 
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Hence 



■U[T{p )] = £ &/(&/») < £Hr°o£ ^rf(.C?Pi/Qj)- (2-87) 



i=l j = l 

Since C™ and C™ satisfy 

lim CJ 1 = lim Cf = 1, 
it is implied by the assumption about / that 



/ 



cm 



f 



<|1-Cf | ( Cl l/fo/fc)! +c m /q j+ c 3 ) 



for sufficiently large n. Therefore 
Aijqj 



< 



< 



3=1 

|i-CTI 



Cf 



ci X 4y&l/(Pj/&)l + c 2 X - 1 + c 3 £ 

J = l J'=l 3=1 

n n n 

C l £ 1o\f(Pj/Uj)\ + C 2 £ Pj + C 3 £ 9j 

3=1 3=1 3=1 



where we have used < Ay < 1 to obtain the last inequality. Since the first term 
in the brace in the last expression is finite in n — > oo limit by the assumption of the 
absolute convergence of U [po] and all the other terms in the brace are finite, 



3=1 1 

Moreover, by the absolute convergence of U [po], it is easily shown that 

1 \ ™ 



lim 

n — >oo 



Thus 



OO CXD 



■W[%)]<EE^«i^i/«i)' 
i=i j=i 



(2.88) 



We can interchange sum over i and sum over j in the right hand side of (2. 
it converges absolutely by the absolute convergence of U [po] ■ Hence 

oo 

-U[T(p Q )} <J2l3f(P3/l3) = ~U\P0]. 
3=1 



□ 
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Proof of the generalized second law 

Let us combine theorem [| with theorem ^ to prove the generalized second law. In 
theorem set the linear map T, the convex function f(x) and the density matrices 
Poo and poo as follows. 

T = T, 
f(x) = xlnx, 

Poo = Pthi^BHi^BH): 

Poo = Pth(pBH,n BH ). (2.89) 

Note that it is theorem || that makes such a setting possible. Hence, if an initial 
state po and the corresponding final state T(po) satisfy 



\pO,Pth{pBH,&BH)] = [T(po),Pth{/3BH,^BH)} = 



(2.90) 



and U[po] converges absolutely, theorem |?] can be applied to the system of the 
quasi- stationary black hole and the scalar field around it. Now 

U[p ] = Tr [p lnp ] - f3 BH (Tr [Ep ] - n BH Tr [L z p Q ]) -InZ 
= U [p^Pbh^bh) -InZ, 



U[T(p 



Tr [polnpo] - Pbh (Tr 



Epo 



n BH Tr 



L z p~o 



= U[T(po);p B H,ilBH)-lnZ, 



\nZ 
(2.91) 



where 



and 



E 



{n jT } V i / 

Y [Y^™* ) \{ T kf})({ T ht}l 
{n iT } V i / 

Y Y 71 ^ \{ n ip})(i n ip}i 



{ n iP } 



Thus the inequality ( |2.8l| ) in this case is itself, which in turn is equivalent to 

the generalized second law. Finally, theorem ^ proves the generalized second law 
for a quasi-stationary black hole which is formed by gravitational collapse, provided 
that an initial density matrix pg of the scalar field satisfies the above assumptions. 
For example, it is guaranteed by lemma || that if po is diagonal in the basis 
then T(po) is also diagonal in the basis {|{?\p})} and ( 2.90] ) is satisfied. The as- 
sumption of the absolute convergence of U[po; (3bh,^bh) holds whenever initial 
state po at ^~ contains at most finite number of excitations. Therefore the assump- 
tions are satisfied when po is diagonal in the basis {K?"^})} and contains at most 
finite number of excitations. 



Here we have to admit that In Z in Eqs. (2.91 ) diverges formally due to the infi- 
nite volume of the system. However, it should be possible to avoid this divergence by 
a suitable regularization. It will be valuable to analyze what kind of regularization 
does well. 
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2.3.3 Concluding remark 

Now we make a comments on Frolov and Page's statement that their proof of the 
generalized second law may be applied to the case of the black hole formed by 
gravitational collapse ]52| . Their proof for a quasi-stationary eternal black hole is 
based on the following two assumptions: (1) a state of matter fields on the past 
horizon is thermal one; (2) a set of radiation modes on the past horizon and a set of 
modes on the past null infinity are quantum mechanically uncorrelated. These two 
assumptions are reasonable for the eternal case since a black hole emit a thermal 
radiation. In the case of a black hole formed by gravitational collapse, we might 
expect that things would go well by simply replacing the past horizon with a null 



surface at a moment of a formation of a horizon (v = vq surface in Figure 2.1). 



However, a state of the matter fields on the past horizon is completely determined 



by a state of the fields before the horizon formation (v < vq in Figure 2.1), in 
which there is no causal effect of the existence of future horizon. Since the essential 
origin of the thermal radiation from a black hole is the existence of the horizon, 
the state of the fields on the null surface has not to be a thermal one. Hence the 
assumption (1) does not hold in this case. Although the above replacement may be 
the most extreme one, a replacement of the past horizon by an intermediate null 
hypersurface causes an intermediate violation of the assumption (1) and (2) due to 
the correlation between modes on the future null infinity and modes on the past 
null infinity located after the horizon formation. The correlation can be seen in 



( 2.65 ) explicitly in the case of the replacement by the future null infinity. Thus we 
conclude that their proof can not be applied to the case of the black hole formed 
by a gravitational collapse. 

Finally we discuss a generalization of our proof to a dynamical background. 
For the case of a dynamical background, Pbh and Qbh are changed from time 
to time by a possible backreaction. Thus, to prove the generalized second law for 
the dynamical background, we have to generalize theorem |^ to the dynamical case 
consistently with the backreaction. Once this can be achieved, theorem [?], combined 



with the quasi- local first law of black hole dynamics derived in subsection 2.2.2, 
seems useful to prove the generalized second law for the dynamical background. 

Here we have to admit that, in the dynamical situation, entropy S ma tter of 
matter fields might lose its physical meaning, while dynamical black hole entropy 



can be defined as a quarter of trapping horizon area (see subsection 2.2.1). However, 
if the dynamical version of the generalized second law would be proved for some 
definition of S ma uer then, by integrating it from an initial stationary state to a final 
stationary state, we would be able to obtain the generalized second law for finite 
changes of black hole parameters. (Note that in the proof given in this section, it has 
been assumed implicitly by the quasi-stationarity that ASbh "C Sbh, &Mbh -C 
Mbh, etc.) For example, if we take a black hole as the initial state and a flat 
spacetime as the final one, then the finite version of the generalized second law 
insists that matter entropy is increased by evaporation of the initial black hole and 
that the produced entropy is greater than the initial black hole entropy. Thus, 
the generalization is necessary for a detailed investigation of the information loss 
problem. 



Chapter 3 

Black hole entropy 



3.1 D-brane statistical- mechanics 

3.1.1 Black brane solution in the type IIB superstring 

The low-energy effective theory of type IIB superstring contains, as its bosonic 
part, a metric field, a dilaton field, a R-R field and NS-NS fields. By setting the 
NS-NS fields to be zero, we obtain the following low-energy effective action in the 
10-dimcnsional Einstein frame. 



1 



16ttGi 



7l0„ 



(3.1) 



where <f> and H are the dilaton field and the R-R three-form field strength, respec- 
tively. (The 10-dimensional Newton's constant Gio is defined so that the dilaton <f> 
vanishes asymptotically. ) 

For this effective theory we consider toroidal compactification to five dimensions 
with an S 1 of length 2ttR, a T 4 of four-volume (2tt) 4 V, and momentum along the 
S 1 : we assume the 10-dimensional metric of the form 



e ^x+4>)/3 g (5) dx n dx v + e W( d x 5 + A^dx») 2 + e 2 * ^{dX, 



(3.2) 



where = 0, 1, • • • , 4, i = 6, • • • , 9, and all fields depend only on x M . Here we also 
assume that X$ is periodically identified with period 2itR and that each of Xj is 
identified with period Note that the conformal factor e _2 ( 4x+! ^/ 3 makes 

(5) 

the metric gjludx^dx 11 be in the 5-dimensional Einstein frame, and that the field 
becomes a U(l) gauge field in 5-dimensions. 

In Ref. [p8l a six parameter family of black brane solutions of the equation of 



motion following from the action (3.1) was given. The metric is of the form (3.2) 
with 



g$dx»dx" = -f-w(l- r -±y t 2 +f/ 3 



p 2 X 



„2-0 _ 



f -i r l sinh2cr ^ 

A 1/4 / 2 - 1/4 , 

,-3/4 ,-1/4 , 
Jl J 2 J3, 



r 2 
— 



dr 2 + r 2 dQf 3) 



(3.3) 
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where 

Tq sinh 2 a 
h = H ^ ' 

r 2 sinh 2 7 
J 2 = 1 + ^ > 

rg sinh 2 cr 

/3 = 1 + ^ > 

/ = /1/2/3. (3.4) 

The solution is parameterized by the six independent quantities a, 7, a, r^, R and 
V. 

As seen as a black hole in 5-dimensions, this black branc solution has two charges 
Qi and Q 2 associated with the R-R field H and the charge N associated with 
the Kaluza-Klein gauge field A^. These charges are written in terms of the six 
parameters as 

Vr 2 

Qi = — 2- sinh 2a, 
2ff 

r 2 

Q 5 = sinh 2 7 , 
2.9 

R 2 Vr 2 

N = fLllo wah2a, (3.5) 

where g is the 10-dimensional string coupling, by which the 10-dimensional Newton's 
constant is written as do = &n 6 g 2 . Note that, from the general argument of Kaluza- 
Klein compactification, the charge N is related to the momentum P around the S 1 
as P = N/R. The 5-dimensional black hole also has charges Py, and P x related to 
the asymptotic fall-off of ip and x, which are given by 

RVr 2 1 
P$ = 2 ° [cosh 2a — -(cosh 2a + cosh 27)], 
2g 2 



RVr 2 

P x = (cosh 2a -cosh 2 7 ). (3.6) 

25 



In 10-dimcnsions P^ and P x are pressures which describe how the energy changes 
for iscntropic variations in R and V. The ADM energy associated with this 5- 
dimensional black hole is 

RVr 2 

E = —f- (cosh 2a + cosh 2 7 + cosh 2a). (3.7) 
25 

Hence the black brane solution can be parameterized by the five charges (Qi, Q5, 
N, P^, P x ) and the ADM energy E. 

For this solution as a 5-dimensional black hole the Bekenstein-Hawking entropy 
Sbh and the Hawking temperature Tbh are given by 

Sbh = = MV^l~+V^r)(Vn~i + V / ^l)(^V5 + ^V§), 



1 nR ( 1 1 

?BH 2 



f^ + ^Vv / ^+V / ^l)(v / ^5 + y^), (3-8) 



where A is the area of the horizon and G5 is the 5-dimensional Newton's constant 
given by G10 = Gc,(2ir) 5 RV . Here N\, Aj, A 5 , Ag, nj, and n_R are defined by 



Q 1 = N\ — Nf, 
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Q 5 = N 5 -N- 5 , 
N = n L - n R , 

P^P = ~{N 1 + Ni)-^-(N 6 + N E ) + ^(n Il + n R ), 

P x = X(N 1 + N l )-—(N 5 + N- 5 ), 
9 9 

E = -(N 1 + N l ) + —(N 5 + N- 5 ) + hn L +n R ). (3.9) 
9 9 R 



The scales of compactification are written in terms of these quantities as 

R = 



2 \ 1/4 



N x Ni 



V = — — . (3.10) 

The black brane solution characterized by (Ni, Ni, iVg, ul, n R ) can be 
interpreted as a configuration of strings and solitons in the type IIB superstring 
theory f6q| , provided that interactions among string and soliton can be neglected. 
The configuration is composed of N\ "constituent" D-strings wrapped on the S 1 , 
Ni "constituent" anti-D-strings wrapped on the S 1 , N$ "constituent" D-fivebranes 
wrapped on the T 5 = T 4 x S 1 , and N% "constituent" anti- D-fivebranes wrapped 
on the T 5 . Open strings on the D-branes have momentum along the 5* 1 so that 
Pl = til I R is the momentum along the S 1 which is a sum over left-moving massless 
modes of open strings and that P R = n R /R is the momentum which is a sum over 
right-moving massless modes. This interpretation is based on the following three 
facts: 

(i) A single-wound D-string or anti-D-string has 



Q 1 = ±1,Q 5 = Q,N = 0, 

— P =- E= - 

9 9 9 



where plus sign is for the D-string and minus sign is for the anti-D-string. 
(ii) A single-wound D-fivebrane or anti-D-fivebrane has 



Oi = 0,Q 5 = ±l,iV = 0, 

™L p - RV f-M 

29' X ~ 9 ' 9 



P 4 , = - — ,P X = —,E=—, (3.12) 



where plus sign is for the D-fivebrane and minus sign is for the anti-D- 
fivebrane. 

(iii) A left- or right- moving string with momentum P = ztn/R along the S 1 has 



= 0, Q 5 = 0, N = ±n, 



Ti TL 

P^ = -,P x = 0,E=-, (3.13) 
where plus sign is for the left-mover and minus sign is for the right-mover. 



In the following arguments, we consider the case that there are no anti-D-branes: 
Nf = TVg = 0, and thus Qi = Ni, Q5 = N&. In this case the Bekenstein-Hawking 
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entropy and the Hawking temperature are given by 




(3.14) 



The purpose of the remaining part of this section is to explain these expressions of 
the entropy and the temperature by using the D-brane technology. 



3.1.2 Number of microscopic states 

Consider a set of Qi single-wound D-strings wrapped on the S 1 and a set of Q$ 
single- wound D-fivebranes wrapped on T 5 — T 4 x S 1 . The D-strings may be con- 
nected up to form a set of multiply-wound D-strings, which is composed of Nql^ 
D-strings of length 2irRqi (qi — 1, 2, • • •). The D-fivebranes may also be connected 

(5) 

up to form a set of multiply-wound D-fivebranes, which is composed of Nq 5 D- 
fivebranes of length 2irRqs (q$ — 1, 2, • • •) along the S 1 . This general configuration 
of D-strings and D-fivebranes was first considered in Ref. ^3|. Note that, by the 
conservation of charges Qi and Q 5 , the integers and N^' must satisfy the 

following constraints: 

Qi = X>w£\ 

91 

95 

since a D-string which winds q\ times has charge (Qi, Q5)—(qi, 0) and a D-fivebrane 
which winds q§ times has charge (Qi, Q,5) = (0, 95). 

For example, the configuration considered by Callan and Maldacena |7(J is the 
case that A^ = Qi, = Q 5 and all other and N^' are zero: no D-branes 



are connected up. (Figure |3.l|(a) is a schematic picture of this situation for Q\ = 2, 



?5 = 3.) On the contrary Maldacena and Susskind M considered the configuration 



such that Nq^ — Nq^ — 1 and all oth er and Nqy are zero: one long D-string 



and one long D-fivebrane. (Figure 3.1(b) is a schematic picture of this situation for 
Qi = 2, Qg = 3.) The general class of configurations introduced above includes 
more abundant situations. For example, the situation in Figure pT|(c) is an example 



for Qi = 2, Q 5 = 3: 7V 2 (1) = N[ 5) = 7V 2 (5) = 1 and others are zero. 

As explained in the previous subsection, the charges til and tir are implemented 
in the D-brane picture so that Pj, = tir/R is the momentum along the S* 1 which 
is a sum over left-moving massless modes of open strings on the D-branes and that 
Pr = ul/R is the momentum which is a sum over right-moving massless modes. 

There are several types of open strings on the D-brane configuration: those 
with both boundaries being on a common D-string; those with one boundary being 
on a D-string and another boundary being on a different D-string; those with one 
boundary being on a D-string and another boundary being on a D-fivebrane; etc. 
Among them we only consider a set of strings each of which connects a D-string and 
a D-fivebrane, since if many such strings are excited then all other strings become 



massive 71 . As in the situation considered in Ref. |70||, each string connecting a 



D-string and a D-fivebrane has 4 bosonic and 4 fermionic degrees of freedom. Hence 
we have AN^ N^} bosonic and ANqy Nqy fermionic degrees of freedom for strings 
with one boundary on one of Nq^ qi-wound D-strings and another boundary on 

(5) 

one of Nq 5 q§ -wound D-fivebranes. 
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X 6 



(b) 

















D 


} /m 



,x 5 



(c) 



X 6 



Figure 3.1: (a) The situation considered by Callan and Maldacena is shown for 
Qi = 2, Q$ = 3. There are three D-fivebranes wrapped on T 5 (X$-g are coordinates 
for the torus.) and two D-strings wrapped on S 1 of the torus. (X 5 is a coordinate 
for the S 1 .) (b) The situation considered by Maldacena and Susskind is shown 
for Qi — 2, Q5 = 3. There are a triple-wound D-fivebrane and a double-wound 
D-string. (c) Another example of the generalized configurations of D-branes for 
Qi = 2, Qh = 3. 
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For each of fermionic or bosonic degrees of freedom of the string connecting a 
gi-wound D-string and a 95-wound D-hvebrane, the spectrum of the momentum 
along the S 1 is 



('+' for left moving 



R(qi,qs)LCM 
' —' for right moving), 



(n = l,2,---) 



(3.16) 
(3.17) 



where (3i,3s)iCM is the least common multiple of q\ and q^ 1 since the boundary 
condition of the strings is 



X 5 ~X 5 + 2nR( qi ,q 5 ) LC M. 



(3.18) 



We regard the momentum Pl = n^jR (Pr = —ur/R) along the S 1 as a sum 
over all left (right) moving massless modes of the strings. It is evident from the 
assumption of neglecting interactions between strings that the number d(n,L, tir) of 
string states consistent with given values of ul and tir is a product of the numbers 
d(rii) and d(nn) of left- and right-moving string states: 



d{n L ,n R ) = d(n L )d(n R ). 



(3.19) 



At this point it is easy to show that d^^) satisfies the following relation for an 
arbitrary value of w satisfying < w < 1. 



z«^=nnn 

1=0 qi 95 n=l 



1 



w 



n4A r(i)Ar(5) 



I — ((j»/(?l*)t(!M 



(3.20) 



where = l/Qi\Q§\. It is evident that d(0) = 1. Note that the denominator in 
the r.h.s. represents a partition function for the bosonic modes of open strings and 
the numerator represents a partition function for the fermionic modes. 

Note that d(ni,) can be considered as a residue of the r.h.s. of Eq. ( 3.20| ) divided 
by x l+1 , where x — loVd'Qs 1 , Hence, d(ni,) can be written as a complex integral: 



d(n L ) 



^-- (m) nnn 



91 95 n=l 



1 + w"/(9i>95)lca/ 

2 _ w ^/(qi,q5)LCM 



— !— (p dw exp[h(w)], 
2-ki J 



where 



91,95 n 



1 _|_ u ,™/(9l,95)iCM 
1 _ w n/(qi,q 5 )LCM 



(til + 1) lnui. 



(3.21) 



(3.22) 



The contour in the integration w.r.t. a; is a closed curve surrounding x = in 
complex x-plane and, the contour in the integration w.r.t. it; is a closed curve 
surrounding w = in complex w-plane. 

Now we give an asymptotic formula for d(riL ) in the limit 



< min (qi 7 q 5 )LCM, 
(91,95)6/ 



(3.23) 



where I = {(qi, qs^Nq^ Nq^ ^ 0}. For this purpose we use an asymptotic expres- 
sion of h(w). Let us consider a function g(x) defined by 



?(*) = n 



n=l L " 



1 +X n 
1 - X™ 



(3.24) 
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where < x < 1 is assumed. For x *~ 1, we can estimate hig{x) as 



hxg{x) = ^[ln(l + z«)-ln(l- 

n=l 
1 



In a; 



dy[\n(l + e~ y ) - ln(l - e~ v )] 



o 



1 



7T 

In a; 4 



(3.25) 



where we have estimated the summation w.r.t. n by the integration w.r.t. y — 
— n\nx. Hence, we obtain the following asymptotic expression of h(w) for w 1 /™ 1 ~ 1, 
where m = mm (quq5)eI (qi,q 5 ) L CM- 



h(w) 



91,95 



(til + 1) lnu>. 



(3.26) 



We shall estimate the integration w.r.t. w in (|3.21| ) by using this expression and the 
suddle-point method. When the condition ( |3.23 ) is satisfied, it is shown by using 
this asymptotic expression of h(w) that exp[h(w)] has a saddle point given by 



lnw 1/r 



n L + 1 



< 



7T , /E 01 , n N^N^ ( qi , q 5 ) LCM ( ^ 7T / Q X Q S 



my n L ~\~ 1 



(3.27) 



Note that w 1 /™ ~ 1 is guaranteed for this saddle point by the condition ( 3.2 3j ) . 
Thus, we obtain the following asymptotic formula of din^) by using the saddle 
point method. 



d(riL) ~ exp 



LCM 



91 95 



From this expression it is easily shown that 



d(n L ) < exp(27iVrt L (2i(25), 



(3.28) 



(3.29) 



where the bound is saturated if and only if all {q\, q§) (G /) are relatively prime. 
Therefore, for a configuration satisfying 



min(n_L,7ifl) (9i,9s)e/ 



< min (qi,q 5 )LCM, 



(3.30) 



\nd(nL,nji) is bounded from above by the Bekenstein-Hawking entropy given by 
( P4l ): 

\nd(n L ,n R ) < S B h, (3.31) 

where the bound is saturated if and only if all ((71,(75) (G /) are relatively prime. 
This is one of the main results in Ref. (4^] and is a generalization of the results in 
Refs. |7^, |6^] to the more abundant configurations of D-strings and D-fivebranes. 
When ur = 0, which corresponds to an extremal black hole, the condition 



(3.3C) does not hold since the l.h.s diverges. Nonetheless, even in this case, the 
bound ( 3.31| ) does hold if the condition (3.23) is satisfied, since d(0) = 1. The 
bound is saturated if and only if all (qi, q§) (s /) are relatively prime. 



CHAPTER 3. BLACK HOLE ENTROPY 



39 



3.1.3 Canonical ensemble of open strings 



The Hawking process can be described in the D-brane picture as a decay process 
of non-BPS excitations of D-branes |70| [72[ [73) . A collision of a right-moving open 
string with a left-moving one on the D-branes results in an emission of a closed 
string leaving away from the D-branes, which is interpreted as Hawking radiation. 
The spectrum of the emission of the closed string can, in principle, be obtained 
from decay rates of the string, provided that the initial state of the open strings on 
the D-branes is given. On the contrary, when we do not know anything about the 
initial state of the open strings, it seems that the best way to evaluate the spectrum 
is summing up decay rates over all consistent initial states of open strings on the 
D-branes. Formally, the summation should be done by using a microcanonical 
ensemble with fixed til and hr. Since this ideal summation is not easy, we adopt 
an approximation. Our strategy of the approximation to perform the summation 
is to replace the microcanonical ensemble by the canonical ensemble, in which each 
expectation value of til and riR coincides with the given fixed value in the original 
microcanonical ensemble. It is easy to see that the corresponding spectrum of closed 
strings is approximately the thermal one with temperature being the same as the 
temperature of this canonical ensemble. 

Thus, in this subsection we consider a canonical ensemble of ope n str ings on 
the D-brane configuration introduce d in the be ginni ng of subsection 3.1.2 . From 
arguments in the derivation of Eqs. ( 3. IS ) and ( 3.20 ), we can obtain the partition 
function of the canonical ensemble as 



Z((3,a) 
Zl,r{P,o) 



Z L (f3,a)Z R ((3,a), 



nnn 

qi qs «=i 



1 



-Pen(qi,q5)-apr l (.qi,qs) 



1 — e -/3en(9l,«5)-aPn(9l,95) 



4AT< 1 > AT< 5 > 



(3.32) 



where (3 is inverse temperature and a is chemical potential w.r.t. the total momen- 
tum along the S 1 . Here p n is given by ( 3.16| ) (the plus sign is for Zl and the minus 
sign is for Zr) and e„ = \p n \- The quantities (3 and a should be determined by 
requiring 



I, , d\nZ((3,a) 
= dp ' 



d\nZ{P,a) 
da 



(3.33) 



Note that these requirements originate from the last and the third equalities in 
(|3.9|), respectively. 

By using the formula ( 3.25| ), we can obtain an asymptotic expression of Zl^r((3, a) 
in the limit ( 3.30 ). The result is 



h\Z L . R (l3, a) 



n 2 R 



0±a- 

qi 95 

('+' for ' L' ■ ' -' for 'R'). 
From the condition ( |3.33| ), /3 and a are determined as 



LCMi 



(3.34) 



P + a 
/3-a 



ttR\ 



ttR\ 



<LCM 



riR 



(3.35) 
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Thus, for a configuration satisfying (3.30), the inverse temperature (3 and the en- 



tropy S of the effective canonical ensemble are 




LCM, 



91 95 



^lnZ(/3,q) d\nZ(0,a) 
S = ^Z{M-f>—^ a - 



27r ( v / nl+ V^r) 



LCAf 

w lnd(n L , n R ), (3.36) 

where d(riL, riij) is the number of micro-states evaluated by using the microcanonical 
ensemble in the previous subsection. In fact, the last equality is trivial from the 
general arguments in statistical mechanics: entropy in microcanonical ensemble 
can be calculated approximately by the corresponding canonical ensemble. The 
temperature is bounded from below by the Hawking temperature given by 



(13.141 ): 

/T 1 > T BH , (3.37) 

where the the bound is saturated if and only if all (91,95) (€ /) are relatively 
prime. Note that this condition for the equality is the same as that for the equality 



in (3.31). Here remember that the temperature (3 1 coincides with the temperature 



of the thermal spectrum of the closed string emission from the D-branes. 
3.1.4 Summary and speculations 

In this section the D-brane statistical-mechanics has been reviewed by using a con- 
figuration of D-strings and D-fivebranes introduced in Ref. p3[ . We have considered 
a set of multiply-wound D-strings, which is composed of Nql^ D-strings of length 
2nRqi (qi = 1, 2, • • ■) and a set of multiply-wound D-fivebranes, which is composed 
of Nq^ D-fivebr anes of length 2irRq 5 (q 5 — 1,2,- • •) along the S 1 . For configura- 



tions satisfying ( 3.30 ), the number of microscopic states d(n£,rt^) of open strings 
on the D-branes is bounded from above by exponential of the Bekenstein-Hawking 
entropy Sbh of the corresponding black hole, and the temperature /3 _1 of a decay 
of D-brane excitations to closed strings is bounded fr om b elow by t he Hawking 



temperature Tbh of the corresponding black hole: Eqs.( 3.31 ) and ( 3.37 ). Note that 
d(nL,nR) and are evaluated microscopically while Sbh and Tbh are defined 
in terms of macroscopic quantities (area and surface gravity of the horizon of the 



corresponding black hole). The bounds ( 3.31 ) and ( 3.37 ) are saturated if and only 
if all (gi, 95) (s /) are relatively prime. 

Thus several speculations may be possible. 

• Inside a black hole characterized by the four parameters (Qi, Q5, til, tir), 
some dynamical processes may occur. The processes may be described in the 
D-brane picture: D-branes repeat fission and fusion to settle down to one of 
the states for which all (qi, 55) (€ /) are relatively prime. 

• During the process the microscopic entropy increases to reach the Bckcnstcin- 
Hawking entropy of the corresponding black hole. Moreover the temperature 
of closed string radiation from the D-branes decreases to reach the Hawking 
temperature of the black hole. 



These speculations may be significant to investigate the microstates of dynamical 
black holes. For example let us consider a merger of two black holes B\ and B 2 
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and suppose that B\ corresponds to a D-brane configuration {-/Vg* , N^} and 



B2 corresponds to a configuration {N^ 1 ', Nq£ }■ Just after merging, the large 



black hole B formed by the merger corresponds to a configuration {Nq^ = Nq] 1 ^ + 
Nqf 1 ^ , Nq^ = Nql 5 ' + Nq$ } , provided that directions of the D-strings are the same 
for B\ and £?2- In general the last configuration does not saturate the bounds (3.31) 
and (3.37) even if the configurations for B\ and Bi saturate the bounds. Thus, in 
general just after the merger the microscopic entropy of B does not agree with the 
corresponding Bekenstein-Hawking entropy and the temperature of the closed string 
radiation does not agree with the corresponding Hawking temperature. However, 
after a sufficiently long time the D-branes' fission and fusion settle the entropy and 
the temperature to the Bekenstein-Hawking entropy and the Hawking temperature. 



3.2 Brick wall model 



The entanglement interpretation, which will be investigated in detail in section 3.3 , 
seems to be implicit in, and is certainly closely related to a pioneering calculation 
done by Gerard 'tHooft [[55) in 1985. He considered the statistical thermodynamics 
of quantum fields in the Hartle-Hawking state (i.e. having the Hawking temperature 
Tbh at large radii) propagating on a fixed Schwarzschild background of mass M . To 
control divergences, 'tHooft introduced a "brick wall" — actually a static spherical 
mirror at which the fields are required to satisfy Dirichlet or Neumann boundary 
conditions — with radius a little larger than the gravitational radius 2M. He found, 
in addition to the expected volume-dependent thermodynamical quantities describ- 
ing hot fields in a nearly flat space, additional wall contributions proportional to 
the area. These contributions are, however, also proportional to a~ 2 , where a is the 
proper altitude of the wall above the gravitational radius, and thus diverge in the 
limit a — > 0. For a specific choice of a (which depends on the number of fields, etc., 
but is generally of order l p i), 'tHooft was able to recover the Bekenstein-Hawking 
formula with the correct coefficient. 

However, this calculation raises a number of questions which have caused many, 
including 'tHooft himself, to have reservations about its validity and consistency. 

(a) Sbh is here obtained as a one-loop effect, originating from thermal excitations 
of the quantum fields. Does this material contribution to Sbh have to be 
added to the zero-loop Gibbons-Hawking contribution which arises from the 
gravitational part of the action and already by itself accounts for the full value 
of S BH ? [@ 

(b) The ambient quantum fields were assumed to be in the Hartle-Hawking state. 
Their stress-energy should therefore be bounded (of order M~ 4 in Planck 
units) near the gravitational radius, and negligibly small for large masses. 
However, 'tHooft's calculation assigns to them enormous (Planck-level) energy 
densities near the wall. 

(c) The integrated field energy gives a wall contribution to the mass 

AM = (3.38) 

when a is adjusted to give the correct value of Sbh- This suggests a sub- 
stantial gravitational back-reaction p5| and that the assumption of a fixed 
geometrical background may be inconsistent Q Jr6|. 

Our main purpose in this section is to point out that these difficulties are only 
apparent and easy to resolve Ejj. The basic remark is that the brick-wall model 
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strictly interpreted does not represent a black hole. It represents the exterior of a 
starlike object with a reflecting surface, compressed to nearly (but not quite) its 
gravitational radius. The ground state for quantum fields propagating around this 
star is not the Hartle-Hawking state [jr7| but the Boulware state [|78| , corresponding 
to zero temperature, which has a quite different behavior near the gravitational 
radius. 



In subsection 3.2.1 we summarize essential properties of the Boulware and Hartle- 



Hawking states that play a role in our arguments. In subsection 3.2.2 we sketch the 



physical essence of the brick- wall model by using a particle description of quantum 



fields. A systematic treatment of the model is deferred to subsection 3.2.3, in which 



the results in subsection 3.2.2 are rigorously derived from the quantum field theory 
in curved spacetimes. In subsection |3. 2. 4 we propose a complementary principle be- 



tween the brick wall model and the Gibbons-Hawking instanton. In Appendix A. 3, 
for completeness, we apply the so-called on-shell method to the brick wall model and 
show that in the on-shell method we might miss some physical degrees of freedom. 
Hence, we do not adopt the on-shell method in the main body of this thesis. 

3.2.1 The Boulware and Hartle-Hawking states 

It is useful to begin by summarizing briefly the essential properties of the quantum 
states that will play a role in our discussion. 

In a curved spacetime there is no unique choice of time coordinate. Different 
choices lead to different definitions of positive-frequency modes and different ground 
states. 

In any static spacetime with static (Killing) time parameter t, the Boulware state 
\B) is the one annuled by the annihilation operators a km associated with "Killing 
modes" (positive- frequency in t). In an asymptotically flat space, \B) approaches 
the Minkowski vacuum at infinity. 

In the spacetime of a stationary eternal black hole, the Hartle-Hawking state 
| H H) is the one annuled by a^rus, the annihilation operators associated with 
"Kruskal modes" (positive- frequency in the Kruskal lightlike coordinates U, V). 
This state appears empty of "particles" to free falling observers at the horizon, and 
its stress-energy is bounded there (not quite zero, because of irremovable vacuum 
polarization effects). 

If, just for illustrative purposes, we consider a (1 + l)-dimensional spacetime, it 
is easy to give concrete form to these remarks. We consider a spacetime with metric 

dr 2 

ds 2 = -f(r)dt 2 + — , (3.39) 
f{ r ) 

and denote by n(r) the redshifted gravitational force, i.e., the upward acceleration 
a(r) of a stationary test-particle reduced by the redshift factor f 1 ' 2 ^), so that 
«(r) = a /'(?")■ A horizon is characterized by r = r , f(ro) = 0, and its surface 
gravity defined by k q = |/'(r ). 

Quantum effects induce an effective quantum stress-energy T a {> (a, 6, ■ • • = r, t) 



in the background geometry (3.39). If we assume no net energy flux (T[ = 0) — 
thus excluding the Unruh state — T a b is completely specified by a quantum energy 
density p = —T\ and pressure P — T£. These are completely determined (up to a 
boundary condition) by the conservation law T%. b — and the trace anomaly, which 
is 

T a = (3-40) 
for a massless scalar field, with R = —f"(r) for the metric (3.3£). Integration gives 

f(r)P(r) = -^—{n 2 {r) + const.). (3.41) 

247T 
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Different choices of the constant of integration correspond to different boundary 
conditions, i.e., to different quantum states. 

For the Hartle-Hawking state, we require P and p to be bounded at the horizon 
r = r , giving 



Pi 



HH — 



h Kq — K 2 (r) 
24^ W) : 



Phh = P HH + ^f"{r). (3.42) 
When r — * oo this reduces to (setting f(r) — > f ) 

7T 2 

Phh ^ Phh = Qft T BH> 

T BH = hK /2ir, (3.43) 

which is appropriate for one-dimensional scalar radiation at the Hawking tempera- 
ture T B h- 

For the Boulware sta te, th e boundary condition is P = p = when r = oo. The 
integration constant in ( 3.41 ) must vanish, and we find 

p ft K 2 (r) 

B 24tt f(r) ' 

PB = P B + ^f'(r). (3.44) 

If a horizon were present, ps and Pb would diverge there to — oo. 
For the difference of these two stress tensors, 

AT Q b = {TDhh ~ (T b a ) B , (3.45) 

( j3~42j ) and JO^ ) 

give the exactly thermal form 

AP = Ap=^T 2 (r), (3.46) 



where T(r) = Tbh I \[JW) is the local temperature in the Hartle-Hawking state. 
We recall that thermal equilibrium in any static gravitational field requires the local 
temperature T to rise with depth in accordance with Tolman's law |79| 



T^—goo = const. (3.47) 

We have found, for this (1 + l)-dimcnsional example, that the Hartle-Hawking 
state is thermally excited above the zero-temperature (Boulware) ground state to a 
local temperature T(r) which grows without bound near the horizon. Nevertheless, 
it is the Hartle-Hawking state which best approximates what a gravitational theorist 
would call a "vacuum" at the horizon. 

These remarks remain at least qualitatively valid in (3 + l)-dimensions, with ob- 
vious changes arising from the dimensionality. In particular, the (3 + l)-dimensional 
analogue of ( 3.46| ) for a massless scalar field, 

2 

3AP~ Ap~ ^U.T 4 (r), (3.48) 
30^ 3 ' ' 

holds to a very good approximation, both far from the black hole and near the 
horizon. In the intermediate zone there are deviations, but they always remain 
bounded IpOl, and will not affect our considerations. 
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3.2.2 A brief sketch of the brick wall model 

We shall briefly sketch the physical essen ce of the brick-wall model. (A systematic 
treatment is deferred to subsection |3.2.3 ) 



We wish to study the thermodynamics of hot quantum fields confined to the 
outside of a spherical star with a perfectly reflecting surface whose radius r x is a 
little larger than its gravitational radius rg. To keep the total field energy bounded, 
we suppose the system enclosed in a spherical container of radius L»ri. 

It will be sufficiently general to assume for the geometry outside the star a 
spherical background metric of the form 

dr 2 

ds 2 = -f{r)dt 2 + — + rW. (3.49) 
f(r) 

This includes as special cases the Schwarzschild, Reissner-Nordstrom and de Sitter 
geometries, or any combination of these. 

Into this space we introduce a collection of quantum fields, raised to some tem- 
perature Too at large distances, and in thermal equilibrium. The local temperature 



T(r) is then given by Tolman's law (3.47) 



T(r) = T^r 1 ' 2 (3.50) 

and becomes very large when r — > r\ = ro + Ar. We shall presently identify T^ 
with the Hawking temperature Tbh of the horizon r = ro of the exterior metric 
( 3.49| ), continued (illegitimately) into the internal domain r < r\. 



Characteristic wavelengths A of this radiation are small compared to other rel- 
evant length-scales (curvature, size of container) in the regions of interest to us. 
Near the star's surface, 

\~h/T= f^K/Too « r . (3.51) 

Elsewhere in the large container, at large distances from the star, 

/ ~ 1, A ~ h/T^ ~ ro < L. (3.52) 

Therefore, a particle description should be a good approximation to the statistical 
thermodynamics of the fields (Equivalently, one can arrive at this conclusion by 
considering the WKB solution to the wave equation, cf. 'tHooft |35| and subsection 

EM: 



For particles of rest-mass m, energy E, 3-momentump and 3- velocity v as viewed 
by a local stationary observer, the energy density p, pressure P and entropy density 
s are given by the standard expressions 



P = N 
P 



E ATip 2 dp 
o eP E -eh 3 ' 
AT f 00 vp 4irp 2 dp 



3 J e? E -e h 3 ' 
0{p + P). (3.53) 



Here, as usual, 



E 2 -p 2 =m 2 , v=p/E, (3 — T _1 ; (3.54) 

e is +1 for bosons and —1 for fermions and the factor M takes care of helicities and 
the number of particle species. The total entropy is given by the integral 

S= [ s(r)47rr 2 dr/v7, (3.55) 
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where we have taken account of the proper volume element as given by the metric 
( 3.49] ). The factor J -1 / 2 does not, however, appear in the integral for the gravi- 
tational mass of the thermal excitations (It is canceled, roughly speaking, by 
negative gravitational potential energy): 



p(r)Airr 2 dr. 



(3.56) 



The integrals ( 3.5 5| ) and ( p. 56] ) are dominated by two contributions for large 
container radius L and for small Ar = r\ — tq\ 

(a) A volume term, proportional to ^L 3 , representing the entropy and mass- 
energy of a homogeneous quantum gas in a flat space (since / ~ 1 almost 
everywhere in the container if L/r — > oo) at a uniform temperature T^. 
This is the result that would have been expected, and we do not need to 
consider it in detail. 

(b) Of more interest is the contribution of gas near the inner wall r = r% , which 
we now proceed to study further. We shall find that it is proportional to the 
wall area, and diverging like (Ar) -1 when Ar — > 0. 

Because of the high local temperatures T near the wall for small Ar, we may 
insert the ultrarelativistic approximations 

E 3> m, p ~ E, d~1 



into the integrals (3.53). This gives 



1 



3 dx 



(3.57) 



in Planck units (h = 2nh = 2ir). The purely numerical integral has the value 3! 
multiplied by 1, 7r 4 /90 and |7r 4 /90 for e = 0, 1 and —1 respectively, and we shall 



adopt 3!, absorbing any small discrepancy into TV". Then, from (3.53) 



3W 



rpi 



4W\ 



(3.58) 



in terms of the local temperature T given by (3.50). 

Substituting (3.58) into (3.55) gives for the wall contribution to the total entropy 



C 4A/ A - 7~ 

Owall = — 5-47rr 1 i ( 



2 rp'3 



r!+8 



dr 

JW)' 



(3.59) 



where S is an arbitrary small length subject to Ar <C 5 <C r±. It is useful to express 
this result in terms of the proper altitude a of the inner wall above the horizon 
r = ?*o of the exterior geometry (3.4£). (Since ( 3.49| ) only applies for r > n, the 
physical space does not, of course, contain any horizon.) We assume that /(r) has 
a (simple) zero for r = tq, so we can write 



/(r) ~ 2n (r - r ), k 
where kq is the surface gravity. Then 

f-^dr 



1 



/'(ro)^0 (r^r ), 



Ar = -K a 2 , 



(3.60) 



(3.61) 
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and ( |3.59 ) can be written 



Swall — 



T 



907ra 2 V k o/2tt / 4 



-A 



(3.62) 



in Planck units, where A = 4nrf is the wall area. 

Similarly, we find from ( |3.56 ) and (3.58) that thermal excitations near the wall 
contribute 



AM, 



them,wall — 



To, 



480?ra 2 \k /2tt 

to the gravitational mass of the system. 
The wall contribution to the free energy 

F = AM - T^S 



AT^ 



p 

if 



TV 



AT~ 



14407ra 2 \k /2tt 

The entropy is recoverable from the free energy by the standard prescription 

S wa u — ~dF wa u/dT a 



(3.63) 

(3.64) 
(3.65) 

(3.66) 



(Observe that this is an "off-shell" prescription [g2|, [32|: the geometrical quantities 
A, a and, in particular, the surface gravity kq are kept fixed when the temperature 
is varied in (3.65).) 

Following 'tHooft | j35| , we now introduce a crude cutoff to allow for quantum- 
gravity fluctuations by fixing the wall altitude a so that 



S W aii = Sbh, when T a 



T, 



BH, 



(3.67) 



where the Bekenstein-Hawking entropy Sbh and Hawking temperature Tbh are 
defined to be the purely geometrical quantities in terms of the wall's area A and 
redshifted acceleration (= surface gravity) Ko- From (3.67) and ( 3.62| ), restoring 
conventional units for a moment, we find 



= Z P zVW/907T, 



(3.68) 



so that a is very near the Planck length if the effective number Af of basic quantum 
fields in nature is on the order of 300. 

It is significant and crucial that the normalization (3.68) is universal, depending 
only on fundamental physics, and independent of the mechanical and geometrical 
characteristics of the system. 

With a fixed by (3.68), the wall's free energy (3.64) becomes 



1 / T 



16 \T BH 



ATry. 



(3.69) 



This "off-shell" formula expresses F wa u in terms of three independent variables: the 
temperature and the geometrical characteristics A and Tbh- From ( |3.69 ) we 
can obtain the wall entropy either from the thermodynamical Gibbs relation (3.66) 
(with Too set equal to Tbh after differentiation), or from the Gibbs-Duhem formula 
(3.64) which is equivalent to the statistical-mechanical definition S = — TV(plnp). 
Thus the distinction (8^, [32| between "thermodynamical" and "statistical" entropies 
disappears in this formulation, because the geometrical and thermal variables are 
kept independent. 
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The wall's thermal mass-energy is given "on-shell" (Too = Tbh) by 

3 

AM t herm,wall = -t^ATbH (3.70) 
16 



according to ( 3.63| ) and (3.68). For a wall skirting a Schwarzschild horizon, so that 



T B h = (87TM)- 1 , this reduces to 'tHooft's result Q3.38|) 



As already noted, thermal energy is not the only source of the wall's mass. 
Quantum fields outside the wall have as their ground state the Boulware state, 
which has a negative energy density growing to Planck levels near the wall. On 



shell, this very nearly cancels the t herm al ene rgy d ensity (3.58); their sum is, in 
fact, the Hartle-Hawking value (cf. ( 3.45| ) and ( 3.48) )): 



{T^)ther m ,T^=T BH + {Tp B - (T£) H H , (3-71) 

which remains bounded near horizons, and integrates virtually to zero for a very 
thin layer near the wall. The total gravitational mass of the wall is thus, from ( 3.63|) 
and (p38l) , 



(AM) wall = (AM) thermtWaU + (AM) B ,waU 

= ^AT B h{(T 00 /Tbh) 4 -1), (3.72) 

which vanishes on shell. For a central mass which is large in Planck units, there is 
no appreciable back-reaction of material near the wall on the background geometry 



We may conclude that many earlier concerns [|35j, |7J, [75j were unnecessary: 
'tHooft's brick wall model does provide a perfectly self-consistent description of a 
configuration which is indistinguishable from a black hole to outside observers, and 
which accounts for the Bekenstein-Hawking entropy purely as thermal entropy of 
quantum fields at the Hawking temperature (i.e. in t he Ha rtle-Hawking state), 
providing one accepts the ad hoc but plausible ansatz ( |3.68 ) for a Planck-length 
cutoff near the horizon. 

The model does, however, present us with a feature which is theoretically pos- 
sible but appears strange and counterintuitive from a gravitational theorist's point 
of view. Although the wall is insubstantial (just like a horizon) — i.e., space there is 
practically a vacuum and the local curvature low — it is nevertheless the repository 
of all of the Bekenstein-Hawking entropy in the model. 

It has been argued J|9| that this is just what might be expected of black hole 
entropy in the entanglement picture. Entanglement will arise from virtual pair- 
creation in which one partner is "invisible" and the other "visible" (although only 
temporarily — nearly all get reflected back off the potential barrier). Such virtual 
pairs are all created very near the horizon. Thus, on this picture, the entangle- 
ment entropy (and its divergence) arises almost entirely from the strong correlation 
between nearby field variables on the two sides of the partition, an effect already 
present in flat space Q. 

An alternative (but not necessarily incompatible) possibility is that the concen- 
tration of entropy at the wall is an artifact of the model or of the choice of Fock 
representation (based on a static observer's definition of positive frequency). The 
boundary condition of perfect reflectivity at the wall has no black hole counter- 
part. Moreover, one may well suspect that localization of entanglement entropy 
is not an entirely well-defined concept |47| ] or invariant under changes of the Fock 
representation. 
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3.2.3 The brick wall model reexamined 



In the previous section, we have investigated the statistical mechanics of quantum 
fields in the region r\ < r < L of the spherical background ( 3.49| ) with the Dirich- 
let boundary condition at the boundaries. By using the particle description with 
the local temperature given by the Tolman's law, we have obtained the inner-wall 
contributions of the fields to entropy and thermal energy. When the former is set 
to be equal to the black hole entropy by fixing the cutoff a as (3.68), the later 
becomes comparable with the mass of the background geometry. After that, it has 
been shown that at the Hawking temperature the wall contribution to the thermal 
energy is exactly canceled by the negative energy of the Boulware state, assuming 
implicitly that the ground state of the model is the Boulware state and that the 
gravitational energy appearing in the Einstein equation is a sum of the renormalized 
energy of the Boulware state and the thermal energy of the fields. 

In this section we shall show that these implicit assumptions do hold. In the 
following arguments it will also become clear how the local description used in the 
previous section is derived from the quantum field theory in curved spacetime, which 
is globally defined. 

For simplicity, we consider a real scalar field described by the action 



f 



[ = -- I d^x^-g [g^d^dA + 



On the background given by ( [3.49 ) , the action is reduced to 



(3.73) 



/ = / dtL, 



(3.74) 



with the Lagrangian L given by 



L = / cf a ; ~ 

2 



(3.75) 



Here x % (i — 1,2) are coordinates on the 2-sphere. In order to make our system 
finite let us suppose that two mirror-like boundaries are placed at r = r\ and r = L 
(n < L) , respectively, and investigate the scalar field in the region between the two 
boundaries. In the following arguments we quantize the scalar field with respect to 
the Killing time t. Hence, the ground state obtained below is the Boulware state. 
After the quantization, we investigate the statistical mechanics of the scalar field 
in the Boulware state. It will be shown that the resulting statistical mechanics is 
equivalent to the brick wall model. 

Now let us proceed to the quantization procedure. First, the momentum conju- 
gate to <f>(r, x l ) is 



7r(r, x" 1 ) 

and the Hamiltonian is given by 



Or 



(3.76) 



H = ~ / ,r.r 



1 7T 2 + r 2 VUf (d r (j)) 2 + Vnn ij di ■ 



Next, promote the field <p to an operator and expand it as 
1 



[r,x = 



E 

nlrn 



a n i m <Pni(r)Y lm (x l )e 



i\ —iuj n it 



a llm ( Pnl(r)Yi m {x i )e i ' 



(3.77) 



(3.78) 
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where Yi m (x l ) are real spherical harmonics defined by 
1 



L d t (Vnn^djY^) + i(i + i)Y lm = o, 



J Y lm (x i )Y l , m ,(x i )^J^x 1 )d 2 x = 5 w 5 mml , 



and {<fini{r)} (n = 1, 2, ■ ■ •) is a set of real functions defined below, which is complete 
with respect to the space of L2-functions on the interval r\ < r < L for each I. The 
positive constant uj n i is defined as the corresponding eigenvalue. 



-^d r (r 2 fd r ipnl) 



1(1 + 1) 2 
s — +mi 



Vnl + -y-<£nl 



0, 



(3.79) 



Wnl(n) = iPnl(L) = 0, 



f(r) 



dr 



The corresponding expansion of the operator w(r, x l ) is then: 

- yog) y 



7r(r, x l ) = — i- 




^^^m^K^' .(3.80) 



Hence, the usual equal-time commutation relation 



becomes 



<p(r, x l ),ir(r',x n ) 



Unlm-, a n'l' m > 
l&nlm > tin' I'm' 
u nim' u n'l'm' 



i5{r-r')5 2 (x l - x' 4 ), 



7r(r, x 4 ), 7r(r', x'*) 



= 8 nn >6u>8 v 
= 0, 
= 0. 



= 



(3.81) 



The normal-ordered Hamiltonian is given by 

: H := ^u n ia} nlm a n i r 



nlm 



Thus, the Boulware state \B), which is defined by 

a n i m \B) = 



(3.82) 



(3.83) 



(3.84) 



for v (n, I, m), is an eigenstate of the normal-ordered Hamiltonian with the eigenvalue 
zero. The Hilbert space of all quantum states of the scalar field is constructed as 
a symmetric Fock space on the Boulware state, and the complete basis {\{N n i m })} 
(Nnim = 0, 1, 2, • • •) is defined by 



|{JW> = II 7^ UiS^ \ B ), (3-85) 
and each member of the basis is an eigenstate of the normal-ordered Hamiltonian: 
: H : \{N nlm }) = f £)w n ,iV nIrn J \{N nlm }). (3.86) 

\nlm / 



CHAPTER 3. BLACK HOLE ENTROPY 



50 



Now we shall investigate the statistical mechanics of the quantized scalar field. 
The free energy F is given by 



IV 



[e - Ml] = TT __I 



(3.87) 



nlm 



where floo — T^ 1 is inverse temperature. For explicit calculation of the free energy 
we adopt the WKB approximation. First, we rewrite the mode function (p n i(r) as 



and suppose that the prefactor ip n i(r) varies very slowly: 

d r ipni 



4>7< 



Thence, assuming that 



« |*|, 
d r (r 2 f) 



« \k\, 



« Ifcl" 



the field equation (3.79) of the mode function is reduced to 



k {l,LL>ri 



f 



L 



(3.88) 



(3.89) 



(3.90) 



(3.91) 



Here we mention that the slowly varying condition ( [3.89] ) can be derived from the 
condition ( |3.90 ) and viceversa. The number of modes with frequency less than lo is 
given approximately by 



(3.92) 



g{uj) = J v(l,u){2l + l)dl, 
where v(l,u>) is the number of nodes in the mode with (I, to): 

u(l,u) = - [ y/k 2 (l,uj)dr. 



(3.93) 



Here it is understood that the integration with respect to r and I is taken over those 
values which satisfy n < r < L and k > 0. Thus, when 



d r (r 2 f) 



r 2J 



< 



1 



f(3o 



is satisfied, the free energy is given approximately by 



I r°° 

Poo Jo 



— /8ooW 



dgjuj) 



duj = 



F(r)4Trr 2 dr, 



where the 'free energy density' F(r) is defined by 

1 



F(r) 



Here the 'local inverse temperature' /3(r) is defined by the Tolman's law 



(3.94) 



(3.95) 



(3.96) 
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and E is denned by E = ^Jp 2 + m|. Hence the total energy U (equal to AM t herm 
given by ( 3.56| )) and entropy S are calculated as 

r-L 



u 



Tr 



^(F-.H:) . jj 



d 



df3 



p{r)4wr 2 dr, (3.97) 



S = Tr 

r-L 



^(F-.H:) ln e ^(F-.H:) 



Ho 



-F 



s(r)47rr 2 (ir/\/ /(r), 



(3.98) 



where the 'density' p(r) and the 'entropy density' s(r) are defined by 



p(r) 
s(r) 







d/3(r) 
(3 2 (r) 



(/?(r)F(r)) 



15 A-Kp 2 dp 
o e^W^-l (2tt) 3 ' 



0/8(r) 



P» = /3(r) (p(r) + P(r)) , 



where the 'pressure' P(r) is defined by [] 



P(r) = -F(r) = - 



p 2 /E 4irp 2 dp 
3 J e /3MB-l (2tt) 3 ' 



(3. 



(3.100) 



These expressions are exactly same as expressions ( 3.53 ) for the local quantities in 
the statistical mechanics of gas of particles. 

Thus, we have show n that the local description of the statistical mechanics 
used in subsection 3.2.2 is equivalent to that of the quantized field in the curved 
background, which is defined globally, and whose ground state is the Boulware state. 

The stress energy tensor of the minimally coupled scalar field is given by 



In particular, the (tt)-component is 



-Adt<t>Y + f{d r <t>) 2 + ^™ 1 - ' • ••■-■■ 2 



^d^d^ + m 2 ^ 2 



(3.101) 



(3.102) 



Hence, the contribution AM of the scalar field to the mass of the background 
geometry is equal to the Hamiltonian of the field: 



AM 



TUnr 2 dr = H, 



(3.103) 



where H is given by ( |3.77 ). When we consider the statistical mechanics of the 
hot quantized system, contributions of both vacuum polarization and thermal ex- 
citations must be taken into account. Thus, the contribution to the mass is given 
by 



(AM) = Tr 



,/MF-:if:) AM (ren) 



(3.104) 



where AM' r ™' is an operator defined by the expression ( 3.103j ) with T\ replaced 
by the renormalized stress energy tensor 7^ reTl )*_ From ( 3.103] ), it is easy to show 
that 

AM (ren) = . R . +AMs) (3.105) 
1 To obtain the last expression of P(r) we performed an integration by parts. 
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where : H : is the normal-ordered Hamiltonian given by Q3.83D and AM B is the 
zero-point energy of the Boulware state defined by 



B\T {ren) \\B)^r 2 dr. 



(3.106) 



Hence, (AM) can be decomposed into the contribution of the thermal excitations 
and the contribution from the zero-point energy: 



(AM) = U + AM B , 



(3.107) 



where U is given by (3.97) and equal to AM therm defined in (3.56). 

Finally, we have shown that the gravitational mass appearing in the Einstein 
equation is the sum of the energy of the thermal excitation and the mass-energy of 
the Boulware state. Therefore, as shown in subsection 3.2.2, the wall contribution 
to the total gravitational mass is zero on shell (Too = Tbh) and the backreaction 
can be neglected. Here, we mention that the corresponding thermal state on shell is 
called a topped-up Boulware state |3(| , and can be considered as a generalization to 
spacetimes not necessarily containing a black hole of the Hartle-Hawking state J77j . 



3.2.4 Complementarity 

Attempts to provide a microscopic explanation of the Bekenstein-Hawking entropy 
Sbh initially stemmed from two quite different directions. 

Gibbons and Hawking |2^| took the view that Sbh is of topological origin, 
depending crucially on the presence of a horizon. They showed that Sbh emerges 
as a boundary contribution to the geometrical part of the Euclidean action. (A 
non-extremal horizon is represented by a regular point in the Euclidean sector, so 
the presence of a horizon corresponds to the absence of an inner boundary in this 
sector.) 

'tHooft J35| sought the origin of Sbh in the thermal entropy of ambient quantum 
fields raised to the Hawking temperature. He derived an expression which is indeed 
proportional to the area, but with a diverging coefficient which has to be regulated 
by interposing a "brick wall" just above the gravitational radius and adjusting its 
altitude by hand to reproduce Sbh with the correct coefficient. 

In addition, the brick wall model appears to have several problematical features — 
large thermal energy densities near the wall, producing a substantial mass correction 
from thermal excitations — which have raised questions about its self-consistency as 
a model in which gravitational back-reaction is neglected. 

We have shown that such caveats are seen to be unfounded once the ground 
state of the model is identified correctly. Since there are no horizons above the 
brick wall, the ground state is the Boulware state, whose negative energy almost 
exactly neutralizes the positive energy of the thermal excitations. 'tHooft 's model 
is thus a perfectly self-consistent description of a configuration which to outside 
observers appears as a black hole but does not actually contain horizons. 

It is a fairly widely held opinion (e.g. [^3[ Q) that the entropy contributed 
by thermal excitations or entanglement is a one-loop correction to the zero-loop 
(or "classical") Gibbons-Hawking contribution. The viewpoint advocated in this 
section appears (at least superficially) quite different. We view these two entropy 
sources — (a) brick wall, no horizon, strong thermal excitations near the wall, Boul- 
ware ground state; and (b) black hole, horizon, weak (Hartle-Hawking) stress - 
energy near the horizon, Hartle-Hawking ground state — as ultimately equivalent 
but mutually exclusive (complementary in the sense of Bohr) descriptions of what 
is externally virtually the same physical situation. The near-vacuum experienced 
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by free-falling observers near the horizon is eccentrically but defensibly explain- 
able, in terms of the description (a), as a delicate cancellation between a large 
thermal energy and an equally large and negative ground-state energy — just as the 
Minkowski vacuum is explainable to a uniformly accelerated observer as a thermal 
excitation above his negative-energy (Rindler) ground state. (This corresponds to 
setting /(r) = r in the (1 + l)-dimensional example treated in subsection 3.2.1.) 

That the entropy of thermal excitations can single-handedly account for Sbh 
without cutoffs or other ad hoc adjustments can be shown by a thermodynamical 
argument [j36f . One considers the reversible quasi-static contraction of a massive 
thin spherical shell toward its gravitational radius. The exterior ground state is the 
Boulware state, whose stress-energy diverges to large negative values in the limit. To 
neutralize the resulting backreaction, the exterior is filled with thermal radiation 
to produce a "topped-up" Boulware state (TUB) whose temperature equals the 
acceleration temperature at the shell's radius. To maintain thermal equilibrium 
(and hence applicability of the first law) the shell itself must be raised to the same 
temperature. The first law of thermodynamics then shows that the shell's entropy 
approaches Sbh (in the non-extremal case) for essentially arbitrary equations of 
state. Thus, the (shell + TUB) configuration passes smoothly to a black hole + 
Hartlc-Hawking state in the limit. 

It thus appears that one has two complementary descriptions, (a) and (b), of 
physics near an event horizon, corresponding to different Fock representations, i.e., 
different definitions of positive frequency and ground state. The Bogoliubov trans- 
formation that links these representations is known |8j|. However, because of the 
infinite number of field modes, the two ground states are unitarily inequivalent |86|. 
This signals some kind of phase transition (formation of a condensate) in the passage 
between description (a), which explains Sbh as a thermal effect, and description 
(b), which explains it as geometry. We know that a condensation actually does 
occur at this point; it is more usually called gravitational collapse. 

It will be interesting to explore the deeper implications of these connections. 



3.3 Entanglement entropy and thermodynamics 

As explained several times, entanglement entropy is often speculated as a strong 
candidate for the origin of the black-hole entropy. To judge whether this speculation 
is true or not, it is effective to investigate the whole structure of thermodynamics 
obtained from the entanglement entropy, rather than just to examine the apparent 
structure of the entropy alone or to compare it with that of the black hole entropy. 
It is because entropy acquires a physical significance only when it is related to the 
energy and the temperature of a system. From this point of view, we construct a 'en- 
tanglement thermodynamics' by introducing an entanglement energy and compare 
it with the black- hole thermodynamics. 

Our strategy of the 'construction of entanglement thermodynamics' is as follows. 
(See chapter |l| for the 'construction of black hole thermodynamics' for Schwarzschild 
black holes.) 

1. First we give concepts and definitions of 'entanglement energy'. (We follow 
the usual definition for entanglement entropy.) 

2. Next we calculate entanglement entropy and entanglement energy for tractable 
models. 

3. Finally, we obtain 'entanglement temperature' T ent by assuming the following 
relation analogous to the first law of thermodynamics. 



SE en t — T en tSS en t, 



(3.108) 
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where S en t and E ent denote entanglement entropy and energy, respectively. 
We call this relation the first law of entanglement thermodynamics. 



This section is organized as follows. In subsection 3.3.1 we review the concept of 



the entanglement entropy. In subsection 3.3.2 we propose four definitions of entan- 



glement energy and present general formulas for calculating the energy. In subsec- 
tion 3.3.3 explicit evaluations of entanglement entropy and energy are performed for 



some tract able m odels with the help of the formulas prepared in subsection 3.3.2 . In 
subsection 3.3.4 we construct entanglement thermodynamics and compare it with 
the black hole thermodynamics. 



3.3.1 Entanglement entropy 

In this subsection we review the definition and basic properties of the entanglement 
entropy. 



Definition of entanglement entropy 

Let T be a Hilbert space constructed from two Hilbert spaces T\ and T as 

T = T X ®T 2 , (3.109) 

where (§) denotes a tensor product followed by a suitable completion. We call an 
element u € T "prime if u can be written as u = v ® w with v € T\ and w G T 2 . For 
example, u = vi®w\-\-2vi®W2-\-V2®wi-\-2v2®W2 is prime since u can be represented 
as w = (vi+v 2 )(3{wi+ 2u>2 ) ■ On the other hand u — vi<£)Wi+V2®w 2 is not prime if 
neither v\ and v 2 nor wi and u>2 are linearly dependent. The entanglement entropy 
Sent ■ F — ► R+ = {non-negative real numbers} defined below can be regarded as a 
measure of the non-prime nature of an element of T = J-\®!F2- 

First of all, from an element u of T with unit norm we construct an operator p 
('density operator') by 



pv = (u, v)u v € J 7 , 

where (u,v) is the inner product which is antilinear with respect to u. 
context p represents a 'pure state'. 

From p we define another operator ('reduced density operator') p2 by 



P2V = ^2 ( e * ® fa ^ e i®y) y y e F 2 , 



(3.110) 
In this 



(3.111) 



where {e^} and {fj} are orthonormal bases of T\ and T2 respectively. Note that 



Tr 2 (p 2 A) - Tr [p(l ® A)} 



for an arbitrary bounded operator A on ^2- 

Finally we define the entanglement entropy with respect to 



p as 



Sent [p] = -fcsTl-2 [p 2 In p 2 ] 



(3.112) 



(3.113) 



We can totally exchange the roles played by T\ and T2 in Eq.( 3.111 ) and Eq.( 3.113] ). 
The entanglement entropy is so defined as to be invariant under the exchange of 
T\ and T 2 when p corresponds to a pure state, i.e., when p is given by Eq.( 3.lfo| ). 
(See Appendix A. 4 for the proof of this property.) 
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A simple example 

As a simple example, let us consider spin states for a system consisting of an electron 
and a proton. We take T\ = {\ f) e , | j) e } for an electron and = {| |)p> | 1} P } for 
a proton, where 'f is for 'up', while 'J.' is for 'down'. Then T — T\ ® Ti is spanned 
by 

{| T)e ® I T)p, I T)e ® I Dp, I 4>e ® I Dp, |i)e® I Dp}- 

Now let us consider a state 

|$ = (a|t>e + /3|i> e )®(7lT>p + *U>p), 
l«| 2 + |/3| 2 = |7| 2 + H 2 = l, 



which is clearly a prime state. According to Eq.(3.111), we then get 

_ / \a\ 2 a(3* 
Pe -{ a*P |/?| 2 

Here 'e' is for 'electron'. By a suitable diagonalization of this matrix, it is easy to 
see that S en t = 0. We can exchange the roles between 'electron' and 'proton': we 
then get 

_ / | 7 | 2 7 5* " 
Pp ~ \ 7 *<5 \6\ 2 

('p' is for 'proton') which again leads to S en t = 0. 
On the contrary, an s-state 

\4>') = a|t)e®|l) P + /3|4)e®|T) P , 

\a\ 2 + |/3| 2 = l,a/3^ 
is not a prime state. For this state the reduced density operators are given by 

O ( W? \ D ( 1/3,2 

pe ~ \ |/3| 2 ) > pp ~ { 
Therefore we get S ent = -fc B (|a| 2 In |a| 2 + |/3| 2 In |/3| 2 ) > 0. 

Formula of entanglement entropy 

Let us consider a system of coupled harmonic oscillators {q A } (A — 1, ••• ,ntot) 
described by the Lagrangian, 

L = ^6 AB q A q B - l -V AB q A q B . (3.114) 

Here 5 A b is Kronecker's delta symbol^; V is a real-symmetric, positive-definite ma- 
trix which does not depend on {g" 4 }- We have introduced a(> 0) as a fundamental 
length characterizing the system.^] The corresponding Hamiltonian becomes 

H tot = ^5 AB PAPB + \v AB q A q B , (3.115) 
la I 

where p A = o,S AB q B is the canonical momentum conjugate to q A . 

2 Prom now on, we choose the units % = c = 1 and apply Einstein's summation convention 
unless otherwise stated. 

3 Thus |<?^} are treated as dimension- free quantities in the present units. 
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Firstly we calculate the wave function ({g' 4 } |0) of the ground state |0). Note 
that Eq.( 3.115) can be written as 



H tot = ^S AB ( PA + tW Ac q C ) ( PB " iW BD q D ) + ^TrW 



(3.116) 



by using the commutation relation [V^ps] = i5g. Here W is a symmetric matrix 
satisfying (W 2 )ab — clVab- The ambiguity in sign is fixed by requiring W to be 
positive definite. Thus, 

W = VaV. (3.117) 
Now ({<? A } |0) is given as a solution to 



d 



since pa is expressed as The solution is 



jt-a+WabQ* <{/}|0>=0, 



(3.118) 



({^}|0)=(det^ 



1/4 



(3.119) 



which is normalized with respect to the standard Lebesgue measure dq 1 ■ ■ ■ dq n 
The density matrix po corresponding to this ground state is represented as 



(U A } \po\ W B }) = ({q A }\0)(0\{q' B }) 



dct I — 



1/2 



x exp 



\W AB (q A q B + q' A q' B ) 



(3.120) 



Now we split {q A } into two subsystems, {q a } (a = 1, •• ■ ,ub) and {q a } (a — 
n B + 1, ■ • ■ ,n to t)- (We assign the labels '1' and '2' to the former and the latter 
subsystems, respectively.) Then we obtain the reduced density matrix associated 
with the subsystem 2 (the subsystem 1), by taking the partial trace of po w.r.t. the 
subsystem 1 (the subsystem 2): 



({q a }\P2\{q ,(3 }) 



Y[dq c ({q^q a }\ Po \{q\q'P}} 



det — 

7T 



x exp 



1/2 



exp 



(B T A~ l B) {q-q'Y (q-q') 



'\0 



(3.121) 



and 



({q a }\pi\W b }) = 



N 



II dq^{q\q a }\p,\{q'\qP}) 
7— n+1 



det — 

7T 



exp 



-\A< ab {q a q b + q< a q< b ) 



x exp 



-±(BD- 1 fl r ) a6 ( ? -«0-(?-9') 6 



(3.122) 
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where A, B, D, A' and D' are defined by 



T , . , / A a {, B a a 

{W AB ) = ( ( V T 



(B ) a b D a p 

A 1 = A — BD~ 1 B T , 

D' = D — B T A~ 1 B . (3.123) 

(The superscript T denotes transposition.) Note that A T = A and D T = D. Here 
we mention that, if we define A, B and D by 

W-^l f £ T , ab % a3 ), (3.124) 

then 

A = A'' 1 , 
D = D'~\ 

B = -A'^BD^ 1 = -A~ L BD'~\ (3.125) 



(B T ) ab D a P J ' 



Now the entanglement entropy S ent := — Trp 2 hip2 is given as follows ]37[ |38| . 
Let {Aj} (i = 1, • • • , N — Ub) be the eigenvalues of a positive definite symmetric 
matrix n A, 

A := D l / 2 B T A- 1 BD 1 ' 2 . (3.126) 
Then it is easily shown that entanglement entropy is given by 

N—riB 



S t = In/Xi-In(l-W). ( 3 - 127 ) 

1 - Mi 

where /i 4 := A^ 1 (-y/1 + A» - l) 2 . (Note that < ^ < 1.) 
Relevance to black hole entropy 

In the case of black-hole physics, the presence of the event horizon causes a natural 
decomposition of a Hilbert space T of all states of matter fields to a tensor product 
of the state spaces inside and outside a black hole as Eq. ( 3.109| ). For example, 



let us take a scalar field. We can suppose that its one-particle Hilbert space Ji is 
decomposed as 

H = Hi®H 2 , (3.128) 

where H.i is a space of mode functions with supports inside the horizon and 7i 2 is a 
space of mode functions with supports outside the horizon. Then we can construct 
new Hilbert spaces ('Fock spaces') J 7 , T\ and from Tt, Hi and H 2 , respectively, 
as 



syra 



Ti = C©Wi©(Wi®Wi) aBTO ©.--, 

T 2 = C®H2®(H 2 m 2 ) sym ®---, (3.129) 



4 The corresponding expression in rcf.Q ("A a (, := (M~ 1 ) ac N cb ") reads A = DB T A l B in 
the present notation. This definition does not giv e a sy mmetric matrix and should be replaced by 



«A_ ab : = (A/- 1 /2)<ic A r cd ( M -i/2)di,.. j namely Eq.(3.126) 
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where (• • -) S y m denotes the symmetrization. Now these three Hilbert spaces satisfy 



the relation (3.109). Hence the entanglement entropy S en t is defined by the proce- 



dure given at the beginning of this subsection (Eqs.(3.110)-(3.113)) for each state 
in T . 

The entanglement entropy S en t originates from a tensor product structure of 



the Hilbert space as Eq.(3.109), which is caused by the existence of the boundary 



between two regions (the event horizon) through Eq.(3.128). Furthermore the sym- 
metric property of S en t between T\ and Ti mentioned before also suggests that 
S en t is related with a boundary between two regions. In fact S en t turns out to 
be proportional to the area of such a boundary (a model for the event horizo n) i n 



simple models discussed below. In view of the Bekenstein-Hawking formula (1.1) 
thus, the entanglement entropy has a nature similar to the black hole entropy. 

The relevance of the entanglement entropy to the black hole entropy is also 
suggested by the following observation. Let us consider a free scalar field on a back- 
ground geometry describing a gravitational collapse to a black hole. We compare 
the black hole entropy and the entanglement entropy for this system. We begin with 
the black hole entropy. In the initial region of the spacetime, there is no horizon 
and the entropy around this region can be regarded as zero. In the final region, on 
the other hand, there is an event horizon so that the black-hole possesses non-zero 
entropy. As for the entanglement entropy, the existence of the event horizon natu- 
rally divides the Hilbert space T of all states of the scalar field into T\®T-x. Thus, 
the scalar field in some pure state possesses non-zero entanglement entropy. In this 
manner, we observe that the black-hole entropy and the entanglement entropy come 
from the same origin, i.e. the existence of the event horizon. This is the reason 
why the entanglement entropy is regarded as one of the potential candidates for the 
origin of the black-hole entropy. 

Simple models in Minkowski spacetime 

The relation between the entanglement entropy and the black hole entropy was ana- 
lyzed in terms of simple tractable models by Bombelli, et. al. |37| and Srednicki pS| . 
They considered a free scalar field on a fiat spacelike hypersurfacc embedded in a 
4-dimcnsional Minkowski spacetime, and calculated the entanglement entropy for 
a division of the hypersurface into two regions with a common boundary B. Here 
the two regions and B are, respectively, the models of the interior, the exterior 
of the black holes and the horizon. Ref. [[37j chooses B to be a 2-dimensional flat 
surface R 2 , and the matter state to be the ground state, showing that the resulting 
entanglement entropy becomes proportional to the area of B. Ref.|38| chooses B 
to be a 2-sphere S* 2 in R 3 , and chooses the two regions to be the interior and the 
exterior of the sphere. The matter state is chosen to be the ground state. Then it 
is shown that the resulting entanglement entropy is again proportional to the area 
of B. 

Both of the results can be expressed as 

S ent [p ] ~k B Ms-r -J, (3.130) 

where po is the ground-state density matrix, A is area of the boundary, a is a cutoff 
length, and Afs is a dimensionless numerical constant of order unity. In particular, 
Afs = 0.30 for B = S 2 j3Sj . This coincides with the Bekenstein-Hawking formula if 
the cut-off length a is chosen as 

(3.131) 




where l p i is the Planck length. Here note that a depends only on the Planck length. 



CHAPTER 3. BLACK HOLE ENTROPY 



■59 



Note that the cutoff length ( 3.131 ) is almost same as the cutoff length ( 3.68| ) for 
the brick wall model. 

3.3.2 Entanglement energy 

In this subsection we define entanglement energy to construct entanglement thermo- 
dynamics. We give four possible definitions of entanglement energy. The difference 
between them comes from the difference in the way to formulate the reduction of 
a system (caused by, for instance, the formation of an event horizon) . In the first 
to third definitions we assume that some operators drop out from the set of all 
observables, while the state of a total system is regarded as unchanged. In the 
fourth definition, on the contrary, we assume that the state undergoes a change in 
the course of the reduction of the system (so that the density matrix of the system 
changes actually), while operators are regarded as unchanged. Since at present we 
cannot judge whether and which one of these treatments reflects the true process of 
reduction, the best way is to investigate all options. As we shall see in subsection 



3.3.4, the universal behavior of the entanglement thermodynamics does not depend 
on the choice of the entanglement energy. 

Let us consider a system described by a Fock space T constructed from a onc- 
particlc Hubert space TL in the previous section. Let Htot be a total Hamiltonian 
acting on T . We assume that the Hamiltonian Htot is naturally decomposed as 

Htot = H i + H 2 + Hi n t j (3.132) 

where Hi and H 2 are parts acting on T\ and T%, respectively, and H int is a part 
representing the interaction of two regions. 

The first to third definitions of the entanglement energy 

The first to third definitions of entanglement energy follow when we regard that 
the operators connecting the two subsystems drop out from the set of observables 
(when, for instance, an event horizon is formed), while the state of the system is 
regarded as unchanged. To be more precise, we assume that H\ and H 2 remain to 
be observables but that Hi nt is no longer an observable. In this case it is natural 
to define the entanglement energy by one of the following three. 

(a) E end = (:H! :> = Tr[: H x : p], 

(b) E end = (: H 2 :) = Tr[: H 2 : p], 

(c) E ent = (:H 1 :) + (:H 2 :), 

where p is a density matrix of the total system and the two normal orderings mean 
to subtract the minimum eigenvalues of H\ and H 2 respectively. 

The forth definition of the entanglement energy 

Next let us consider the case in which the total density operator p actually changes 
to the product of reduced density operators of each subsystems, p\ and p 2 , (when, 
for instance, an event horizon is formed), while the observables remain unchanged. 
In this case p reduces to p' given by 

p' = Pi®p2 ■ (3.133) 

It is easy to see that the entropy associated with this density matrix becomes 

- fc s Tr [p 1 \np'] = S ent [p] + S' ent [p], (3-134) 



CHAPTER 3. BLACK HOLE ENTROPY 



60 



where S en t[p] and S' ent [p] are entanglement entropy obtained through p\ and p2, 
respectively. S en t [p] and S' ent [p] are identical if p is a pure state (see the argument 
below Eq. Q3.113D ). 

Since we are assuming that the observables do not change, we are led to the 
following forth definition of entanglement energy: 

(d) E ent = (: Htot :) P ' = Tr [: H tot : p'], 

where : — : denotes the usual normal ordering (a subtraction of the ground state 
energy). 



Formula of (: Hi :) for the ground state 

What we should do next is to give formulas of entanglement energy explicitly by 
choosing p as the ground state po of H to t ■ In the next subsection we consider a free 
scalar field and discretize it with some spatial separation for regularization. Since 
the system thus obtained is equivalent to a set of harmonic oscillators, here, we 
give a formulas of entanglement energy for the ground state of coupled harmonic 
oscillators described by the Hamiltonian ( |3.115 ), splitting the total system as in the 
previous subsection. 

First, we derive formulas for the entanglement energies corresponding to the 
definitions (a), (b) and (c). 



Firstly we divide the Hamiltonian (3.11 5|) into three terms as Eq.(3.132): 



Hi = ^S ab PaPb + lv^ ab q a q b 



2a 



TT b [Pa 

la \ 



■ (1) d 



Ho 



Hu 



1 

2a. 
1 

2a 
Htot 



5^ (p a - 
Hi - 



2 

vw. 
Hi 



fsq <r 



■ (2) 6 



— Trw (1) , 
2a 



= Vint a0 a q P , 

where V {1 \ and V mt are blocks in the matrix V given by 



(Vab) = 



V 



(i) 



ab 



(Vint) a t 



(Vint )ab 



v. 



(2) 



(3.135) 



(3.136) 



a/3 



and u/ 1 ) and m/ 2 ' are, respectively, the positive square-roots of aV^ and aV^ 2 \ 
Although there exists freedom in the way of the division, the above division seems 
to be the most natural one. Here and throughout this section we adopt it. 
By rescaling the variables {q A } as 



W 1 ' 2 



q 



c 



BC 



the expression of the density matrix for the vacuum state Eq.( 3.12C ) gets simplified 
as 

r 1 



({^}|po|{g-' s }>=n^ 1/2 exp 



{(q°) 2 + (q' C ) 2 } 



and the normal ordered Hamiltonian : Hi : is represented as 



: H, 



= -—S ab 
2a 



2a 



U 



d 
dq* 

ab ( J?_ 

dq A 



-(1) -C 
W A C 1 







d 



-(1) -D 
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Here the matrices U and u/ 1 ) are denned as 

V AB ._ 5 AC ( w l/2\ gab ( W l/2) g 
V J Co, V / bD 



0« - 
AB ■ — 



DB 
bD 



/ , \ Ca /-, x / , \ bD 

Sac(w- 1 / 2 ) w^(w-^) 5 db . 



Hence the matrix elements of : Hi : p with respect to the basis \q ) are expressed 
as 



({q A }\:H 1 :p \{q' B }) 

= ^{[(^ + i)f/(^)-i; 



AB 



q q 



X 



N 

Hvr-^expHcf) 2 ] 

c=i 

From this we obtain R 



Hi 



/ (ti A }\:H 1 :p \{<f B }) 
-L[aVy b ) (Ar b + A ab 5 ab -2w$6 a ' 



Similarly (: Hi :) is expressed as 
where w^ 2 ^ is the positive square-root of aV^ 2 



(3.137) 



(3.138) 



Formula of (: H to t -)p' for the ground state 

By using the formulas ( 3.122 ) and ( |3.121 ), p' defined by Eq.( 3.133 ) is represented 



as 



(U A }\P'\W B )) = (det^) 1/2 exp 



M AB (q A ~ B • ~> A ~i^ 



[q q +q q 



x exp 



--N AB (q-q') A (q-q') 



l\B 



(3.139) 



where 



(Mas) 
(N ab ) 



Kb o 

D' af3 

(bd~ 1 b t ^ 







' ab 







(3.140) 



We can diagonalize M and N simultaneously by the following non-orthogonal 
transformation: 

B , (3.141) 



/ - q A 



(uM 1 ' 2 ) 

where U is a real orthogonal matrix satisfying 

M- 1/2 NM~ 1/2 = U T \U, 

Ai 

A = 



A 2 



(3.142) 



Sec Eqs. fc. 12.^ )" and (p. 124|) for the definitions of the matrices A, A, D and D. 
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Now in terms of {(T 4 }, H to t is represented as 
Htot = -7T (UMU T ) —r - W Ac f — b + W BD q D )+—TiW, (3.143) 



2a V 



dq 



dq 1 



thus, 



D 



where 



W = um- 1/2 wm- 1,2 u t . 

Hence the density matrix p' is expressed in terms of | { q A } ) a: 

JN 



(3.144) 
(3.145) 



<{^}ip'i{9 ,b }}= n^ i/2ex p 



c=i 



(3.146) 

This density matrix is normalized with respect to the measure dq 1 ■ ■ ■ dq N . 

Now it is easy to calculate the entanglement energy. First the matrix components 
of : H to t '■ p' with respect to {q A } are given by 



({q A }\:H tot :p'\{q B }) 



UMU T - IUM- 1/2 VM- 1/2 U T 

N 



AB 



q A q B 



2a I 



-Tr [W - N/2 - M] X J| tt" 1 / 2 exp [-(q c ) 2 ] . (3.147) 



Hence the entanglement energy (: H to t ') P > is expressed as 



tot •/ p' 



N 

(Y[df){{f}\:H tot :p'\{q B }) 
J c=i 



4o 



Tr [aVM- 1 +M + N- 2W] 



(3.148) 



Here we have used the formula J dx x • Ax exp[— x ■ x\ = ^tt n ' 2 TtA, where N is the 
dimension of x. With the help of the identity Tr [M + N] = TrA + TrD = TrW, we 
finally arrive at the following formula for (: H to t -)p' 



Ht ot :) P ' = — Tr [aVM' 1 - W] 
= -Tr [ViA-r 1 - W- 1 )} 



= -2 Tr 



(3.149) 



Alternative formula of (: Hi 



H 2 :) for the ground state 



We can calculate (: Hi :} + (: H 2 :) for the ground state as a sum of Eqs. (3.137) 
and (3.138). Nonetheless, for a check of numerical calculations, it is useful to give 
an alternative formula of (: Hi :) + (: H 2 :). 



Einstein's summation convention is not applied to Eq.(3.146). 
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In terms of {o' 4 }, the operator : Hi : + : H^ ■ is written 



as 



where w) is defined by 

Sac (rtr 1 / 2 )"'" 

« o6 

o w ( 2 ) Q/3 



(3.150) 



From the expression we obtain 

({f A }|(:IT 1 : + :ir a :)po|{ff B }) 

= {[(«; + - 1)] AB q^q B - Tr [W(to - 1)]} 

x JjTT-Vaexpf-Cg ) 2 ]. (3.152) 

c=i 

Hence we arrive at the following expression (: H\ :) + (: H% :) for 

(:iM + <:# 2 :) = / ( JJ ^)<{? A } 10 #i : + 0po| {? B }) 

= _L Tr T™ 2 ^ 1 - Wl -— Ti[to-W]. (3.153) 
4a 2a 

With the help of the relation Tr[w 2 W^ _1 ] = Tr[al/M _1 ] which follows from the 
definitions of w and M, this formula is simplified as 

(: Hi :) + (: H 2 :) = ^Tr [aWkT 1 - W] - -^Tr [w - W] 

= (:Htot--) P '-^Tr[w-W], (3.154) 



where Eq.( 3.149 ) has been used to obtain the last line. 

3.3.3 Explicit evaluation of the entanglement entropy and 
energy for a tractable model in some stationary space- 
times 

With the help of the formulas derived in the previous two subsections, we now 
calculate entanglement entropy and energy explicitly to construct entanglement 
thermodynamics for a tractable model in Minkowski, Schwarzschild and Reissner- 
Nordstrom spacetimes. 

Model description 

The basic idea of entanglement thermodynamics is to express the thermodynamic 
quantities for a black hole in terms of expectation values of quantum operators 
dependent on the spacetime division as in the statistical mechanics modeling of the 
thermodynamics for ordinary systems. Therefore we must specify how to divide 
spacetime into two regions and with respect to what kind of state the expectation 
values are taken. 
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According to the original idea of entanglement, it is clearly most natural to 
consider a dynamical spacetime describing black hole formation from a nearly flat 
spacetime in the past infinity and divide the spacetime into the regions inside and 
outside the horizon. In this situation, if we start from the asymptotic Minkowski 
vacuum in the past, the entanglement entropy associated with the division of space- 
time by the horizon acquires a clear physical meaning. However, this ideal modeling 
seems difficult. 

Hence, we are obliged to consider a stationary spacetime, and take as the quan- 
tum state a stationary one. To be specific, we consider the vacuum with respect to 
a Killing time in a spherically symmetric static black hole described by a metric of 
the form 

ds 2 = -N(p) 2 dt 2 + dp 2 + r{pf (dB + sin 2 0di/; 2 ) . (3.155) 



The vacuum is well known as the Boulware state. As stated in subsections 3.2.1 



and 3.2.2, the Boulware state has negative energy density and its contribution to 



gravitational mass diverges for a stationary black hole background. It has been 



shown in section 3.2 that the ground state of the brick wall model is the Boulware 
state and that the negative divergence in gravitational mass due to the Boulware 
state is canceled by a positive divergence due to thermal excitations above the Boul- 
ware state. Hence, we can safely state that the brick wall model is (the Boulware 
state) + (thermal excitations). Thermal characters of a black hole system can be 
explained by the latter: entropy and temperature of a black hole are modeled by 
entropy and temperature of the thermal excitations. 

Our present purpose in this section is to explain the thermal characters of the 
black hole system in terms of entanglement. From the above arguments on the 
brick wall model, we expect that entanglement in the Boulware state does well for 
this purpose: entanglement entropy and entanglement temperature might explain 
entropy and temperature of a black hole; the negative divergence in gravitational 
mass due to the Boulware state might be canceled by entanglement energy, pro- 
vided that the sum of the vacuum energy and entanglement energy contributes to 
gravitational mass. 

Here a subtlety occurs: in, for example, a Schwarzschild spacetime, if we require 
that the quantum field is in the Boulware state on the whole extended Kruskal 
spacetime and take the bifurcation surface as the boundary surface, then the en- 
tanglement entropy vanishes because the state is expressed as the tensor product 



of the Boulware states in the regions I and II of Figure 3.2. 

In order to avoid this, we restrict the spacetime into the region I, and replace 
the boundary by a timelike surface S at a proper distance of the order of the cut-off 



length of the theory to the horizon (see Figure 3.2). This prescription corresponds 
to taking fluctuations of the horizon into consideration: quantum fluctuations of 
geometry near a horizon will prevent events closer to the horizon than about the 
Planck length from being seen on the outside. However, we should keep in mind 
that we have no definite criterion regarding the exact position of the boundary. To 
minimize this ambiguity we will also investigate the influence of the variation of 
the boundary position. 

As a matter content we consider a real scalar field described by 

1 



- , v ^dx 4 [d^d^ + ml). (3.156) 

The mass does not play an essential role since a typical length scale controlling 
the entanglement thermodynamics is much smaller than the Compton length of an 
usual field. Therefore we just set = in the following arguments. 
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Figure 3.2: The Kruskal extension of the Schwarzschild spacetime. We consider 
only the region / (the shaded region) . As the boundary £ we take the hypersurface 
r = tb- 



Discretized theory of a scalar field 

Consider a massless real scalar field described by the action 

1 ~ ~2 



dx 4 y ^jd fi (j,d tl (j), 



(3.157) 



where the background geome try is fixed to be a spherically symmetric static space- 
time with the metric ( 3.155 ). For this system we calculate entanglement entropy 
and entanglement energy to construct entanglement thermodynamics using the 
methods developed in the previous two subsections. Those methods are both 
based on a discre te sys tem {q A } (A = 1, 2, • • • , n to t) described by a Hamilto- 
nian of the form ( 3.115 ). For t his di screte system it is easy to divide the whole 
Hilbcrt space T into the form ( |3.10S| ): T is defined as a Fock space constructed 
from {q A } (A — 1, 2, • • • , n to t)', F\ is defined as a Fock space constructed from 
{q a } (a = 1,2, •••,ng); Ti is defined as a Fock space constructed from {q a } 
(a = tib + 1, • • • ,ntot)- In order to apply this scheme to our problem we have 
to co nstruct a discretized theory of the scalar field whose Hamiltonian is of the 
form (3.115). 

First we expand the field in terms of the spherical harmonics as 



<P(j>, = 2^ —^—4>lrn{p)Zl m {0-, 1p), 

Lm 



(3.158) 



where Z lm = \/2^tYi m for to > 0, \/2"<sYi m for to < 0, and Z l0 — Y t0 . Then the 



Hamiltonian corresponding to the Killing time for the action (3.157) is decomposed 



into a direct sum of contributions from each harmonics component H[ m as 



H — 22 



(3.159) 



in i 
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Here Hi rn is given by 



Hlr, 



- 2 ldp 



Nr< dp 



8 /TV 1 / 2 



+ 1(1 + 1) 



, (3.160) 



where iri m (p) is a momentum conjugate to 4>i m (p)- 

Note that for any (l,m) the Hamiltonian ( 3.160 ) of the subsystem has the form 



I I <WmO>) + \ [ dpdp^ lm (p)V [Un \p,p)4>i m {p), (3.161) 



Hl, n — 



where the following algebra of Poisson brackets is understood: 

{<t>lm(p),n>m>(p')} = 6{p-p'), 
{<t>lm(p),<i>l>m>(p')} = 0, 

{n m (p),n'm'(p')} = 0. 



(3.162) 



Each subsystem described by the Hamiltonian ( 3.161 ) can be discretized by the 
following procedure: 



P 

S(p P') 



(A-l/2)a , 
Sab I a , 



(3.163) 



where A, B = 1, 2, • • • and a is a cut-off length. The corresponding Hamiltonian of 
the discretized system is of the form ( 3.115 ) with 



0i m (p) ->■ q , 
n m {p) -* pa/ a , 

W' m W) -> V AB /a 2 



(3.164) 



In this way we obtain a discretized system with the total Hamiltonian ( |3.115| ) with 
the matrix V given by the direct sum 



V = © V {l ' m \ 



(3.165) 



where V^ l,m ^ is independent of m and is explicitly expressed as 



V 



{l,m),A iB 



AB VlmVlm 



Here 



A=l 



( XA+1/2 
^A+l/2 I 



N 



1/2 



A A A 
Irn 



X A 



(3.166) 



x A 

XA+l/2 

N a 



N 



A+l/2 



r(p=(A-l/2)a)/r , 
r(p = Aa)/r Q , 
N(p=(A-l/2)a) , 
N(p = Aa) , 
4>lm(p = (A— l/2)a) . 



(3.167) 
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In the matrix representation V^ l,m ^ is given by the n to t x ntot matrix 

\ 



V 



(l,m) 



AB 



-.(0 



Ax £« 



2a 
72" 



A, 



Aa-i E<? A a 



/ 



-(r /a) 2 A^ 2 



iV 



A-l/2»A-l/2 + Na+1/2%A+1/ 



+-l(l + l)Nix-/ , 



^(r /a) 2 < /2 iV A+1/2 7Vy + 2 ia ; A 1 4 + i/2^i - (3-168) 



where we have imposed the boundary condition 0^*° t+1 =0. In these expressions Tq 
is an arbitrary constant, which we set to be area radius of a horizon^] for convenience 
in the following arguments. 



Spatial division 

We divide the total system by a stationary hypersurface £ defined by r = tb'- 
we split the system {(j)^} (A = 1, • • • , n to t) into the two subsystems, {<f>f m } (a = 
1, • • • , ub) and {0^} (a = Ub + 1, ■ • • , THot)> the area radius re of the boundary is 
given by tb = r{p = ubo). 

Since this division preserves spherical symmetry of the system, we can still 
apply the expansion by harmonics. Thus, entanglement quantities of the system 
are calculated as sums of all contributions from those subsystems, each of which is 
specified by (l,m) and described by the matrix V^ ,m '. 



Convergence of the summation 



Since V^ l ' m ' is independent of m 
and energy are given by 



(— /, —1 + 1, ■ ■ ■ ,1 — 1, 1), the entanglement entropy 



1=0 

£(2/ + 1)2^. 



(3.169) 



1=0 



where S^ t and E^t are entanglement entropy and energy of the subsystem specified 
by (I, m) and independent of m. 

From Eq.(3.16£), one can easily show that 



O((la/r ) 4 ln(7a/r )) as la/r 
O((la/ro)~ 3 ) as la/ro — > 00. 



00, 



(3.170) 



Thus the infinite sums Eq.( 3.169) actually converge so that we can safely truncate 
them at some appropriate I, depending on the accuracy we require and the ratio 
ro/a we set. 

7 For the model in Minkowski spacetime we take ro = a. 
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Figure 3.3: The numerical evaluations of entanglement energy in Minkowski space- 
time. E ent /(r B /a 3 ) is shown as functions of r B /a, where E ent denotes (: H\ :), 
(: Hi :) and (: H to t '■) p' , respectively, and r B = n B a. We have taken ntot — 60 for 
(: H x :) and (: H 2 :), and n tot = 200 for (: H tot :) p ,. 



Numerical calculations in Minkowski spacetime 

Before analyzing our model in black hole geometries, it is instructive to calculate 
entanglement quantities in the simple models in Minkowski spacetime introduced 
by Bombelli et.al. J3?J and Srednicki |3q| . They calculated entanglement entropy for 
divisions of a flat hypersurface in Minkowski spacetime by R 2 and S 2 , respectively. 



Their results are summarized as Eq. ( |3. 130|) 



Here, let us calculate entanglement energy in the model of Srednicki (a division 
by S 2 ) []. The discretized theory of a scalar field corresponding to this model can be 
easily obtained by the above procedure by setting ro = a, r(p) = p and N(p) = 1. 
The discretized system {^^} {A = 1, • • • , n tot ) is described by Hamiltonian of the 



form ( 3.115 ) with the potential given by ( 3.168| ) 



Figure |3.3j shows the result of numerical calculations of (: Hi :), (: H 2 :) and 
Htot '}p> in this model. 
From this figure we see that E ent is almost proportional to r B /a 3 : 



r 2 



Hn) ~ (:# 2 :)~0.1 ., 

7° 



B 



, r 



H tot :) p , ~ 0.4-| . (3.171) 

a 6 



Numerical calculations in Schwarzschild spacetime 

In Schwarzschild spacetime the metric is given by 



ds 2 = - (1 - ^) dt 2 + (l - ^) 1 dr 2 + r 2 (d9 + sin 2 6d^ 2 ) , (3.172) 



r 



8 Entanglement energy in the model of Bombelli et.al (a division by R 2 ) is analyzed in Ap- 
pendix |A.5| . 
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Figure 3.4: The numerical evaluations for S ent of the discretized theory of the 
scalar field in Schwarzschild spacetime. S ent /(r B /a) 2 for n B = 1,2,5 is shown 
as functions of r B /a, where r B = r(p = n B a). We have taken n to t — 100 and 
performed the summation over / up to 10ro/a. 



where r$ is the area radius of the horizon. As the radial coordinate p we take the 
proper distance from the horizon: 



Vy 2 -i + In (y + Vv 2 - i) 



(3.173) 



where the variable y is defined by y = 2r/r — 1. 

Using formulas given in subsection 3.3.1 , we have evaluated S en t numerically. In 
this calculation the out er num erical boundary is set at n to t = 100. The summation 
with respect to / in Eq.( 3.169 ) is taken up to I — [lOro/a] ([ ] is the Gauss symbol). 
From the above asymptotic beha vior of S^ t , this guarantees the accuracy of 10%. 

The result is shown in Figure 3.4. From this figure we see that S en t is propor- 
tional to (r B /a) 2 if we change ro with n B fixed, and its coefficient has only a weak 
dependence on n B . Thus, we get 



( r JL 

V a 



0.3 — 



(3.174) 



This result is essentially the same as the previous result (3.130) for models in 
Minkowski spacetime including the numerical coefficient. This can be understood 
in the following way. 

Let us make a coordinate change from r to x defined by 



r 

ro 



(x + lf 
Ax 



or 
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Figure 3.5: The numerical evaluations of (: Hi :) for the discretized theory of the 
scalar field in Schwarzschild spacetime. (: Hi :)/(r#/a 2 ) for rig — 1,2,5 is shown 
as functions of Tb/cl, where tb = r(p = nsa). We have taken n to t — 100 and 
performed the summation over I up to 10ro/a. 



Then the metric (3.172) is rewritten as 

2 



ds 2 = 



x-1 
x + l 



dr 



2x 



(dx 2 +x 2 dfl 2 ). 



Note that r = 0, rg and oo correspond to x = 0, 1 and oo, respectively. It is easy 
to see that the Hamiltonian is given in this coordinate system as 



H 



^ Him, 



Hi, 



64x 4 (x- 1) 

C (x + iy 

If (x + l) 
2 



-P, 



4x 



I in 



{dmm) + In <Pl„ 



(3.175) 



where £ := r x, and P im and <p lm ^.^^ ^ , lrn .- 6ix . i ^_i- ) 
tylm : = To4>im in terms of <fti m . Here note that the vacuum state is only weakly 
dependent on the prefactor 64x 



are expressed as P; m := 



rp (x+l) 
64 x 4 (x 



and 



(x+iy 



in Eq.( 3.175 ). If we neglect this prefactor, 



the vacuum state is determined by the Hamiltonian which coincides with that for 
the flat spacetime at x = 1. On the other hand, S en t depends on the modes in a 
thin layer around the boundary E, whose typical thickness is a few times of a ~ Zpj. 
Therefore, when S is near the horizon, the value of Sent should be well approximated 
by the flat spacetime value. 



Next, with the helps of formulas developed in subsection 3.3.2 , we have numer- 
ically evaluated (: Hi :), (: H2 :) and (: H to t -)p'- Now we have taken the numerical 
outer boundary at n tot — 100. The truncation in the ^-summation is the same as 
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Figure 3.6: The numerical evaluations of (: H 2 :) for the discretized theory of the 
scalar field in Schwarzschild spacetime. (: H 2 :)/(r B /a 2 ) for n B = 1,2,5 is shown 
as functions of r B /a, where r B = r(p = n B a). We have taken n to t = 100 and 
performed the summation over I up to 10r /a. 
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Figure 3.7: The numerical evaluations of (: H to t :)/>' f° r the discretized theory of 
the scalar field in Schwarzschild spacetime. (: H tot :) p >/(r B /a 2 ) for n B = 1,2,5 is 
shown as functions of r B /a, where r B = r(p = n B a). We have taken ntot = 100 
and performed the summation over I up to 10ro/a. 
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for S en t (up to I = [10ro/a]), which implies that the accuracy is about 10% from 
the above asymptotic estimate for E t 



(0 



The results of numerical calculations are shown in Figure 3.5, Figure 3.6 and 
Figure 3.7. In these figures, E ent /(r B /a 2 ) is plotted as a function of r B /a for 
ns = 1, 2, 5. All of these figures show that E ent is proportional to r B /a 2 : 



H 



Hi:) 
H 2 :) 

tot •) p' 



0.05(n B - l/2)r B /a 2 , 
0.05(n B + l/2)r B /a 2 , 
0.2n B r B /a 2 . 



(3.176) 



From these equations we immediately see that the values of (: H\ :) and (: 
H 2 :) coincide except for a tiny difference independent of n B . This difference is 
understood by the gravitational red-shift: (: Hi :) comes from the modes just inside 
£ while (: H 2 :) originates from the modes just outside S. In the present numerical 
calculations, it means that (: Hi :) and (: H 2 '■} are determined by the modes at 
P = (n B — l/2)a and p = (n B + l/2)a, respectively(see Eq. ( 3.168 )). Hence, taking 
account of the fact that the contribution of each mode to the entanglement energy 
is proportional to the red-shift factor at its location, the ratio of (: Hi :) and (: H 2 :) 
should be approximately given by 



Hi :) : (: H 2 :) ~ N l ! 2 {p = (n B - l/2)a) : N x / 2 {p 
~ (n B - 1/2) : (n B + 1/2). 



(n B + l/2)o) 

(3.177) 



This is consistent with the above numerical result. 

This argument is also supported by the numerical result for the flat spacetime 
model shown in Figure 3.3. In this figure the values of (: Hi :) and (: H 2 :) for a 
massless scalar field in the Minkowski spacetime with S = i?xR = 5 2 xR arc 
plotted. In this case there is no gravitational red-shift effect, so we expect that 
(: Hi :) ~ (: H 2 :) as confirmed by the numerical calculation. 



Numerical calculations in Reissner-Nordstrom spacetime 

In Reissner-Nordstrom spacetime the metric is given by 

IM Q 2 \ , a / 2M Q 2 Y 1 , . 

— + -V )dt 2 + 1 + ] dr z 

r r* I \ r r z 



1 



(de 2 



] 9dip 2 ) 



(3.178) 

where M and Q are the mass and the charge of the black-hole. The area radius of 

the outer horizon ro is 

r = M + yjM 2 - Q 2 . (3.179) 

As the radial coordinate p we take the proper distance from the outer horizon 

p = ^(M 2 - Q 2 )(y 2 - 1) + M In (y + vV - l) , (3.180) 
where the variable y is defined by 

r — M 

= ■ (3-181) 



Using formulas given in subsection 3.3.1, we have evaluated S en t numerically. In 



this calculation the outer numerical boundary is set at n to t — 100, and the boundary 



of the spatial division is fixed at n B = 1. The summation in I in Eq.(3.169) is taken 
up to I — [10ro/a] ([ ] is the Gauss symbol). From the above asymptotic behavior 
of <Sg2{> tru s guarantees the accuracy of 10%. 
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Figure 3.8: The numerical evaluations of S en t for the discretized theory of the 
scalar field in Reissner-Nordstrom spacetime. S en t / (r B / a) 2 for q = 0,0.9,0.999 is 
shown as functions of r B /a, where r B = r(p — n B a). We have taken n tot = 100 
and performed the summation over / up to lOrp/a. 



The result is shown in Figure 3.8. From this figure we see that S en t is propor- 
tional to (r B /a) 2 if we change ro with q = Q/M fixed, and its coefficient has only 
a weak dependence on q. Thus, we get 



0.3 ( — 

a 



(3.182) 



This result is essentially the same as the results in Minkowski and Schwarzschild 
spacetime. 



Next, by using formulas in subsection 3.3.2, we have numerically evaluated (: 
Hi :), (: H2 :) and (: H tot :) p <. We have taken the numerical outer boundary at 
n tot = 100 and the boundary of the spatial division at n B — 1. The truncation in 
the Z-summation is the same as for S ent (up to / = [10r /a]), which implies that the 



accuracy is about 10% from the above asymptotic estimate for E 



(i) 



The results of numerical calculations are shown in Figure 3.9, Figure 3.10 and 
Figure 3.11. In these figures, E ent /(c(q)r B /a 2 ) is plotted as a function of r B /a for 
n B = 1 and q = 0,0.9,0.999, where c(q) is defined by 



c(q) = 



1 + VT 



(3.183) 



All of these figures show that E ent is proportional to c(q)r B /a 2 . The result is 
summarized as 



(: ff a :) 
(: H 2 :) 
Htot '■} p' 



0.025c(q)r B /a 2 , 
0.075c{q)r B /a 2 , 
{).2c{q)r B /a 2 . 



(3.184) 
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Figure 3.9: The numerical evaluations of (: Hi :) for the discretized theory of 
the scalar field in Reissner-Nordstrom spacetime. (: Hi :) /(c(q)r B /a 2 ) for q = 
0, 0.9, 0.999 is shown as functions of r B /a, where r B = r(p = n^a). We have taken 
ntot = 100 and performed the summation over I up to 10r /a. 



(: H 2 :)/(c(q)r B /a 2 
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Figure 3.10: The numerical evaluations of (: H 2 :) for the discretized theory of 
the scalar field in Reissner-Nordstrom spacetime. (: H 2 :)/(c(q)r B /a 2 ) for q = 
0, 0.9, 0.999 is shown as functions of r B /a, where r B = r(p = n B a). We have taken 
ntot = 100 and performed the summation over / up to 10ro/a. 
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Figure 3.11: The numerical evaluations of (: H tot :)„/ for the discretized theory 
of the scalar field in Reissner-Nordstrom spacetime. (: H tot :) p > / (c(q)r B /a 2 ) for 
q = 0,0.9,0.999 is shown as functions of r B /a, where r B = r(p = n B a). We have 
taken n to t — 100 and performed the summation over I up to lOrg/a. 



3.3.4 Comparison: entanglement thermodynamics and black- 
hole thermodynamics 

Entanglement thermodynamics in Minkowski spacetime 

We have introduced four possible definitions of entanglement energy. By combining 
each of them with the entanglement entropy S en t, we obtain entanglement thermo- 
dynamics. 

For this purpose, we consider an infinitesimal process in which the way of the 
division of the Hilbert space Ti. into Hi and "Hi i s changed smoothly, with the 
'initial' state po being fixed. (See subsection 3.3.1 .) Let SS en t and 8E ent be the 



resultant infinitesimal changes in the entanglement entropy and in the entanglement 
energy, respectively. We are dealing with a 1-parameter family of the infinitesimal 
changes for the entanglement. The parameter is chosen to be the area radius of 
the boundary sphere. Thus, the construction o f the t hermodynamics means to use 
the first law of entanglement thermodynami cs (3.108 ) t o deter mine T ent , whic h we 
call entanglement temperature. Combining ( 3.130 ) and ( 3.171 ) with Eq.( 3.108 ), we 
thus get f\ 

k B T ent = ^--, (3.185) 

Ms a 



where Me is a numerical factor in Eq. ( |3 .171 ). Note that entanglement temperatures 
Tent obtained from the four definitions of the entanglement energy coincides up to 

numerical factors of order unity. 

Let us interpret entanglement thermodynamics given by ( |3.130 ), ( |3.17l| ) and 



9 In this subsection, we recover h and c. 
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( 3.185| ). It is helpful to introduce the quantities 



_ (r B \ 2 
n en t = \~) i 

e ent = ^. (3.186) 

Here n en t is regarded as an effective number of degrees of freedom of matter on the 
boundary 2-surface £?, and e en t is a typical energy scale of each degree of freedom 

on B. 

From Eqs. ( |3.130D , ( |3.17l| ) and fl3.185| ), we find that 



Sent ~ kgTl en t, 
Eent ^ ^-ent^enti 
ksTent ~ eent- (3.187) 

Therefore our results can be interpreted as follows]^ : The entanglement entropy is 
a measure for the number of microscopic degrees of freedom on the boundary B; the 
entanglement energy is a measure for the total energy carried by all of the degrees 
of freedom on B\ entanglement temperature is measure for the energy carried by 
each degree of freedom on B. 

Discrepancy between entanglement thermodynamics in Minkowski space- 
time and black-hole thermodynamics 

Now we compare these results with the case of black holes. For that purpose 
we express the black-hole thermodynamics in the same form as in the previous 
subsection. 

Let us introduce the quantities 



rtBH 



Ipl 

e BH = (3.188) 

where re, is area radius of an event horizon of a black hole. We can interpret that 
riBH corresponds to the effective number of degrees of freedom on the event horizon 
and eBH is a typical energy scale for each degree of freedom of matter on the horizon. 
The black-hole thermodynamics can be recast in terms of these quantities as 

SbH ~ ksUBH, 
Ebh ~ "fBHCBHriBH, 
ksTBH ~ Jbhsbh, (3.189) 

where jbh = Ipi/ro- The factor jbh can be understood as a magnification of energy 
due to an addition of gravitational energy or a red-shift factor of temperature since 
\/—gtt ~ Ipl/ ro at r ~ vq + I^Jtq, which corresponds to a stationary observer at 
the proper distance l p i away from the horizon. Thus the following interpretation is 
possible^: The black-hole entropy is a measure for the number of the microscopic 
degrees of freedom on the event horizon; the black-hole energy is a measure at 
infinity for the total energy carried by all of the degrees of freedom on the event 

10 It is safer, however, to regard such an interpretation just as a convenient way of representing 
our results. (See the next section for another interpretation of entangleme nt ent ropy. This note 



in particular applies to the rase of tl 
11 See the footnote after Eq.( 3.187) 



e black-hole thermodynamics (see Eq.( 3.189)) 
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horizon; the black-hole temperature is a measure at infinity for the energy carried 
by each degree of freedom. 



Now we compare the two types of thermodynamics characterized by Eq.(3.187) 
and Eq.( 3.189] ), respectively. Both of them allow the interpretation that they de- 
scribe the behavior of the effective microscopic degrees of freedom on the boundary 
B, or on the horizon. Because of the factor jbh, however, they are hardly under- 
stood in a unified picture. This strongly suggests that an inclusion of gravitational 
effects is necessary for agreement between them. 

The discrepancy is highlighted in the context of the third law of thermodynamics. 
Both types of thermodynamics fail to follow Planck's version of the third law 0, 
but in quite different manners. 

From Eq. fl3.185D , we see that T ent remains constant if =^?& is assumed to be 
constant On the other hand, Eq.(3.f3C) shows that S en t tends to zero as A — ► 0, 
where A = AiTr 2 B is area of the boundary. Therefore we obtain the following A- 
dependence: 



Sent OC A, 
E ent CX A, 
kBT en t OC A°. 



(3.190) 



The system behaves as though it is kept in touch with a thermal bath with temper- 
ature T ent . 



In co ntrast, for the black-hole thermodynamics, Eq.(L5) and Eq.(L6 ) along with 
Eq.(1.7) give the behavior (note that A cx M^ H ) 



Sbh oc A, 
E B h °c VA, 
IzbTbh oc 1/vA, 



(3.191) 



where A — Atttq is area of the horizon. Thus we see that Sbh — > cx) as Tbh 0. 



The discrepancy between Eq.(|3.190| ) and Eq. fl3.19l| ) is quite impressive. On one 
hand, a well-known behavior (3.191) comes from the fundamental properties of the 
black-hole physics. On the other hand, the behavior characterized by Eq.(3.190) is 
also an universal one in any model of the entanglement: the zero-point energy of the 
system has been subtracted in the definitions of entanglement energy, thus only the 
degrees of freedom on the boundary B contributes to E en t, yielding the behavior 
E ent cx A. The definitions of E ent proposed here look quite reasonable though 
other definitions may be possible. The result E ent cx A also looks natural, being 
compatible with the concept of 'entanglement'. At the same time, S ent also behaves 
universally as S ent cx A, which has been the original motivation for investigating 
the relation between Sbh and S ent pA pq, p9|. 



Restoration of the agreement by inclusion of gravity 

Let us discuss a possible restoration of the agreement between entanglement ther- 
modynamics and black-hole thermodynamics by considering gravitational effects. 

Although the behavior (3.187) of entanglement thermodynamics was derived by 
considering models in flat spacetime, it seems very reasonable that we regard the 
quantities S en t, E ent and T ent as those in a black-hole background measured by a 

12 Here we mention that the black hole system does satisfy the third law in the sense of Nernst M, 



Here we are regarding the cut-off scale a as the fundamental constant of the theory, not to 
be varied. 
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stationary observer located at the proper distance a away from the horizon q Since 
Sbh, Ebh and Tbh in (3.189) are quantities measured at infinity, it is behavior of 
S en t, E ent and T ent at infinity that we have to compare with (3.189). 

S en t at infinity probably has the same behavior as that measured by the ob- 
server near the horizon since a number of degrees of freedom seems independent of 
an observer's view-point. That is consistent with the fact that the entanglement 
entropy on Schwarzschild background has the same behavior Sent ~ ksUent p9| , 
On the other hand it seems natural to add the gravitational energy to the entangle- 
ment energy by replacing E ent with \/—g ttE e nt ■ Then entanglement temperature is 
determined b y use o f the first law ( 3.108 ). Thus the inclusion of gravity may alter 
the behavior (|3.187| ) to 



Sent ~ kgrient ? 
Eent ^ 'Yent&entHent 3 
kllTent 00 ^Jent^ent 3 



(3.192) 



where ~f en t = cl/tq. The factor ~f en t represents the gravitational magnification of 
the entanglement energy due to the addition of gravitational energy since on the 
corresponding Schwarzschild background ^J—gtt ~ afro at r ~ tq + a 2 /ro, which 
corresponds to a stationary observer at the proper distance a away from the horizon 
(see the argument below (3.18E)). Here ro is the area radius of the horizon. 

The revised behavior (3.192) shows a complete agreement with (3.189), provided 
that rs — tq and a ~ l p i. Note that the last equality in (3.192) is consistent with 
an interpretation that the entanglement temperature is red-shifted by the factor 
7e„t . Thus the inclusion of gravitational effects restores the agreement between the 
entanglement thermodynamics and the black-hole thermodynamics at least quali- 
tatively. 



Entanglement thermodynamics in Schwarzschild spacetime 



From the numerical results in subsection 3.3.2, entanglement entropy and entangle- 
ment energy in Schwarzschild spacetime are expressed as 



Sent 
E en t 



k B M s (— 



hcM E -4 



\ a 
a 2 



(3.193) 



for all definitions of E en t, where tb is the area radius of the boundary defined by 
tb = r(p = fifla). Here Afs and Me ar e numerical factors of order 1: Ms = 0.3; 
Me = 0.05(n B - 1/2), 0.05(n B + 1/2), 0.2n B for (: H x :), (: H 2 :) and (: H tot :} p ,, 

respectively. 

From these expressions and the first law of entanglement thermodynamics ( 3.108| ), 
entanglement temperature T en t is determined as 



Te: 



■ MtTbh, 



(3.194) 



where Tbh is the Hawking temperature of the background geometry and Mr is a 
constant defined by Mt = 2ttMe/Ms- Numerical values of Mt is Mr — {tib — 
l/2)7r/3, (n B + l/2)7r/3, 4ns7r/3 for (: Hi :), (: iJ 2 :) and (: H to t ') P ', respectively. 

These results have several interesting features. First of all we immediately see 
that the entanglement thermodynamics on the Schwarzschild spacetime shows ex- 
actly the same behavior as the black hole thermodynamics. This behavior is just 
what we expected from the above intuitive argument on gravitational effects: the 



The author thanks Professor T. Jacobson for helpful comments on this point. 
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gravitational redshift effect modifies the area dependence of E ent so as to make the 
entanglement thermodynamics behave just like the black hole thermodynamics. 

Second it should be noted that the temperature T ent becomes independent of 
the cut-off scale a for all definitions of entanglement energy. 

It is also suggestive that the average of entanglement temperatures for (: Hi :) 
and (: H2 :) with tib = 1 gives almost the same value as Tbh- This averaging 
corresponds to averaging out the difference in the red-shift factors for the one-mesh 
'inside' and the one-mesh 'outside' of the boundary. Therefore such an averaging 
may have some meaning. 



Entanglement thermodynamics in Reissner-Nordstrom spacetime 



From the numerical results in subsection 3.3.3, entanglement entropy and entangle 



ment energy in Reissner-Nordstrom spacetime are expressed as 



Sent — k B Ns ( ~ 

\ a 



2 



E ent ~ hcAfEc(q)-^, (3.195) 
or 

for all definitions of E en t, where tb is the area radius of the boundary defined by 
tb = r(p = hbcl). Here Ms and Me are numerical factors of order 1, whose values 
are the same as those for Schwarzschild spacetime, and the coefficient c(q) is a 



function of q = Q/M given by Eq. (3.183). The coefficient c(q) approaches to zero 
in the limit q — > 1. 

Combining these expressions with the first law of entanglement thermodynamics 



(3.10S), we obtain entanglement temperature T ent as 

Tent — MtTbh, (3.196) 

where Tbh is the Hawking temperature of the background geometry and Mr is a 
constant defined by Mr = 2ttMe/Ms- 

Note that both Tbh and T en t become zero in the extremal limit (g — ► 1) of the 
background spacetime. Therefore as in the case of the Schwarzschild spacetime, the 
entanglement thermodynamics in the Reissner-Nordstrom spacetime has the same 
structure as that of the black- hole thermodynamics. 



3.3.5 Concluding remark 

In this section we have constructed entanglement thermodynamics for a massless 
scalar field in Minkowski, Schwarzschild and Reissner-Nordstrom spacetimes. The 
entanglement thermodynamics in Minkowski spacetime differs significantly from 
black-hole thermodynamics. On the contrary, the entanglement thermodynamics 
in Schwarzschild and Reissner-Nordstrom spacetimes has the same structure as that 
of black-hole thermodynamics. In particular, it has been shown that entanglement 
temperature in the Reissner-Nordstrom spacetime approaches zero in the extremal 
limit. 

Our model analysis strongly suggests a tight connection between the entangle- 
ment thermodynamics and the black hole thermodynamics. Of course, our model is 
too simple to give any definite conclusion based on it. In particular, the ambiguity 
in the definition of the energy comes from neglecting backreaction of the quantum 
field on gravity. 

Finally we comment on possible extensions of the entanglement thermodynam- 
ics. The first is the inclusion of a charged field as matter. In particular, it will 
be valuable to analyze the entanglement thermodynamics for a charged field in 
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Reissner-Nordstrom spacetime. In this case we may be able to define the entan- 
glement charge as an expectation value of the charge of the field for the coarse- 
grained state. The second is a generalization to non-spherically-symmetric space- 
times. For example it is expected that construction of entanglement thermodynam- 
ics in Kerr spacetime requires the introduction of a concept of an entanglement 
angular-momentum. 



3.4 A new interpretation of entanglement entropy 



In this section a new interpretation of entanglement entropy is proposed: entangle- 
ment entropy of a pure state with respect to a division of a Hilbert space into two 
subspaces 1 and 2 is an amount of information, which can be transmitted through 
1 and 2 from a system interacting with 1 to another system interacting with 2. 
The transmission medium is quantum entanglement between 1 and 2. In order 
to support the interpretation, suggestive arguments are given: variational princi- 
ples in entanglement thermodynamics and quantum teleportation. It is shown that 
a quantum state having maximal entanglement entropy plays an important role in 
quantum teleportation. Hence, the entanglement entropy is, in some sense, an index 
of efficiency of quantum teleportation. 



In subsection 3.4.1, based on a relation between the entanglement entropy and 
so-called conditional entropy, we propose an interpretation of the entanglement en- 
tropy. In subsection 3.4.2 variational principles in entanglement thermodynamics 
are used to determine quantum states. In particular, a state having maximal entan- 
glement entropy is determined and is used in subsection 3.4.31 to transmit informa- 



tion about an unknown quantum state. Subsection 3.4.4 is devoted to a summary 
of this section and to discuss implications for the information loss problem and 
Hawking radiation. 



3.4.1 Conditional entropy and entanglement entropy 

Entropy plays important roles not only in statistical mechanics but also in informa- 
tion theory. In the latter, entropy of a random experiment, each of whose outcomes 
has an attached probability, represents uncertainty about the outcome before per- 
forming the experiment p7| . Besides the well-known Shannon entropy, there exist 
various definitions of entropies in information theory. For example, the so-called 
conditional entropy of an experiment A on another experiment B is defined by 
H{A\B) — — ^ ab p{a 1 b)\np{a\b), where a and b represent outcomes of A and B, 
respectively, p{a 1 b) is a joint probability of a and 6, and p(a\b) = p(a,b)/p(b) is a 
conditional probability of a on b. Here p(b) is a probability of b. The conditional 
entropy corresponds to an uncertainty about the outcome of A after the experiment 
B is done. In other words it can be regarded as the amount of information about 
A which cannot be known from the experiment B. The quantum analogue of the 
conditional entropy was considered in references (8^, |8p and is calle d the v on Neu- 
mann conditional entropy. Consider a Hilbert space T of the form ( |3. 109 ) and let 



p be a density matrix on T . The von Neumann conditional entropy of p about the 
subsystem 1 on the subsystem 2 is defined by 

%2 = Tr [pa 1]2 ] , (3.197) 

where cjjia = 1 1 <g> In p2 — In p. The von Neumann conditional entropy S 2 \ i of p about 
the subsystem 2 on the subsystem 1 is defined in a similar way. It is expected that 
Sx|2 ( or ^211) represents the amount of the information about the subsystem 1 (or 
2) which cannot be known from 2 (or 1, respectively). 
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The von Neumann conditional entropy can be negative. In fact, it is easy to see 
that 

Si\2 = S 2 \i = —S en t, (3.198) 

if p is a pure state. Hence, if p is a pure state then the conditional entropy is zero 
or negative. Our question now is 'what is the meaning of the negative conditional 
entropy of a pure state?' It might be expected that |«Si| 2 | (= S en t) is the amount 
of the information about 1 (or 2) which can be known from 2 (or 1, respectively). 
However, this statement is not precise. A precise statement is that it is an amount 
of information, which can be transmitted through 1 and 2 from a system interacting 
with 1 to another system interacting with 2. The transmission medium is quantum 
entanglement between 1 and 2. 

The purpose of the remaining part of this section is to give suggestive arguments 
for this statement. 



3.4.2 Variational principles in entanglement thermodynam- 
ics 

In statistical mechanics, the von Neumann entropy is used to determine an equilib- 
rium state: an equilibrium state of an isolated system is determined by maximizing 
the entropy. Thus, we expect that the entanglement entropy may be used to deter- 
mine a quantum state. 

As an illustration we consider a simple system of two particles, each with spin 



1/2: we consider a Hilbert space T of the form (3.109) and denote an orthonormal 
basis of Ti by {| t)i, | (i = 1,2). Let |</>) be an element of T with unit norm 
and expand it as 

10) = a| t)i <8 I T>2 + b\ T)i <8 I 4)2 + c\ |>! <g> | t) 2 + d\ |) a <g> | |) 2 , (3.199) 

where |a| 2 + |6| 2 + |c| 2 + \d\ 2 — 1 is understood. The corresponding reduced density 
matrix is given by 

P2 = (|a| 2 + | c | 2 )|T) 2 2(TI + («&* + crf*)|T)2 2 ai 

+ (a*b + c*d)\ |) 22 (T | + (\b\ 2 + \d\ 2 )\ |) 22 (| | (3.200) 

and the entanglement entropy can be easily calculated from it. The resulting ex- 
pression for the entanglement entropy is 

1 + 1, fl + x\ 1-x, fl-x\ 
Sent = ^ln ^ -— In — , (3.201) 



where x = yl — 4|arf — bc\ 2 . By requiring dS en t/dx — we obtain the condition 
\ad — bc\ = 1/2. Thus a state maximizing the entanglement entropy is 

|0) = ^(|T>i®|i)a-U)i®|T)a) (3.202) 

up to a unitary transformation in T\ and the corresponding maximal value of the 
entanglement entropy is In 2. This state is well known as the EPR state. 

It is notable that the corresponding reduced density matrix p 2 represents the 
microcanonical ensemble. This fact is related to the fact that the maximum of 
entropy gives the microcanonical ensemble in statistical mechanics. Thus, in gen- 
eral, if T\ and Ti have the same finite dimension N then a state maximizing the 
entanglement entropy is written as 



4=E(i n >i®i n >2) ( 3 - 203 ) 

V 71=1 



CHAPTER 3. BLACK HOLE ENTROPY 



82 



up to a unitary transformation in T\, where |n)i and \n)2 [n — 1,2, ■ ■ ■ , N) are 



orthonormal basis of T\ and J- 2 , respectively. (See Appendix A. 6 for a systematic 



derivation. ) In the next section we use the state (3.203) to transmit information 
about an unknown quantum state. 

In statistical mechanics, free energy F = E — TS can also be used to determine a 
statistical state: its minimum corresponds to an equilibrium state of a subsystem in 
contact with a heat bath of temperature T, provided that T is fixed. This variational 
principle in statistical mechanics is based on the following three assumptions. 

1. The total system (the subsystem + the heat bath) obeys the principle of 
maximum of entropy. 

2. Total energy (energy of the subsystem + energy of the heat bath) is conserved. 

3. The first law of thermodynamics holds for the heat bath. 

Is there a correspon ding va riational principle in the quant um system in the Hilbert 



space T of the form ( p.l09|) ? The answer is yes. In section |3.3| a concept of entangle- 
ment energy has been introduced and a thermodynamical structure, which we call 
entanglement thermodynamics, has been constructed by using the entanglement en- 
tropy and the entanglement energy. Thus we expect that entanglement free energy 
F ent defined as follows plays an important role in entanglement thermodynamics. 

Fent — E ent — T ent S en t, (3.204) 

where E ent is the entanglement energy and T ent is a constant. Among several 
options, we adopt the second definition (b) of the entanglement energy given in 



subsection 3.3.2 



E ent = {:H 2 :). (3.205) 

As shown in the following arguments, by minimizing the entanglement free en- 
ergy, we can obtain a state in T characterized by the constant T ent . Before doing 
it, here we consider a physical meaning of the principle of minimum of the entangle- 
ment free energy. Let us introduce another Hilbert space Tbath, which plays the role 
of the heat bath in the above statistical-mechanical consideration, and decompose 
it to the direct product Tbath — Tbathi ® Tbathi- In this situation it is expected that 
the principle of the minimum entanglement free energy corresponds to the following 
situation. 

1. The total system J- to t = T ® Tbath obeys the principle of maximum of the 
entanglement entropy with respect to the decomposition Ttot — Ttotx ®Ttot2, 
where T to ti = T\ 8> Tbathi and T to t2 = T 2 ® T ba th2- 

2. Total entanglement energy (entanglement energy for T + entanglement energy 
for Tbath) is conserved. 



3. The first law of entanglement thermodynamics ( 3.108 ) holds for Tbath- In this 
situation we call the constant T ent the entanglement temperature. 

It must be mentioned here that the variational principle of minimum of the en- 
tanglement free energy is not as fundamental as the principle of maximum of the 
entanglement entropy but is an approximation to the latter principle for a large sys- 
tem. However, like the principle of minimum free energy in statistical mechanics, 
the former principle should be a very useful tool to determine a quantum state. 

We now calculate F ent for the system of two spin- 1/2 particles and minimize it. 
For simplicity we adopt the following Hamiltonian for the subsystem 2: 



;ti 



2 



H 



H 2 
H 2 



IT>2 = e, 
I D2 = 0, 

I D2 = 0, (3.206) 
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where e is a positive constant. The entanglement free energy F ent for the state 
( 3.199| ) is given by 



Fent = e(\a\ 2 + \c\ 2 )+T ent 



1 + x , / 1 + x 
— z — In 



1 — x ( 1 — x 
— ^ — In 



(3.207) 



By minimizing it we obtain the following expression for the state \(f>) up to a unitary 
transformation in T\ . 



1 

7z 



e/2T ^IT)i®IT)2 + |l)i®|l) 2 



(3.208) 



where Z = e ~ e/T ' nt + 1. 

The corresponding reduced density matrix p2 on represents a canonical en- 
semble with temperature T ent . This fact is related to the fact that the principle 
of minimum of free energy results in a canonical ensemble in statistical mechan- 
ics. Thus, in general, if T\ and Ti have the same finite dimension N then a state 
minimizing the entanglement free energy is written as 



1 N 



-B„/2T e „ 



'In)! ® |ra) s 



(3.209) 



up to a unitary transformation in T\, where Z = Xl^Li e - ^™^""*, and E n and \n/2 
(n = 1,2, ■• ■ ,N) are eigenvalues and orthonormalized eigenstates of the normal- 
ordered Hamiltonian of the subsystem 2. (See Appendix A. 6 for a systematic 

derivation.) 

The state ( |3.20S| ) can be obtained also from another version of the principle of 
maximum of the entanglement entropy: if we maximize S en t with E ent fixed then 
the state ( 3.209D is obtained. In this case, the constant T ent is determined so that 
the entanglement energy coincides with the fixed value. 

Note that in Eq. ( 3.209| ) the infinite-dimensional l imit N — > oo can be taken, 
provided that T ent is bounded. In this limit, the state ( 3.209 ) has the same form as 
those appearing in the thermo field dynamics of black holes [p5| and the quantum 
field theory on a collapsing star background Q. In fact, if we can set the value 
of the entanglement temperature of Tbath to be the black hole temperature then 
the state 



section 



3.209 ) in the limit completely coincides with those in Refs. [|S5|, 



In 



3.2 it has been shown numerically that the entanglement temperature for 
a real massless scalar field in Schwarzschild and Reissner-Nordstrom spacetimes is 
finite and equal to the black hole temperature of the background geometry up to a 
numerical constant of order 1. The finiteness of the entanglement temperature in 
the black hole spacetimes is a result of cancellation of divergences in entanglement 
entropy and entanglement energy pq] . Thus, the finiteness is preserved even in the 
limit of zero cutoff length (a — > 0) . 



3.4.3 Quantum teleportation 

In Ref. ]9(| Bennet et al. proposed a method of teleportation of an unknown 
quantum state from one place to another. It is called quantum teleportation. In 
their method the information about the quantum state is separated into a 'quantum 
channel' and a 'classical channel', and each channel is sent separately from a sender 
"Alice" to a receiver "Bob" . What is important is that the quantum channel is sent 
in a superluminal way by using a quantum correlation or entanglement, while the 
classical channel is transmitted at most in the speed of light. Here we mention that 
causality is not violated in an informational sense since Bob cannot obtain any useful 
information about the unknown state before the arrival of the classical channel. 
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Hence Alice has to deliver the classical channel to Bob without fail. On the contrary 
she does not need to worry about whether the information in the quantum channel 
arrives at Bob's hand since the arrival is guaranteed by the quantum mechanics. It is 
notable that recently quantum teleportation was confirmed by experiments |92| . 

In this section we generalize the arguments in Ref. ]9C| to more abundant situ- 
ations and try to reformulate it in terms of the entanglement entropy. 



Let us consider a Hilbert space T of the form (3.10E) with T\ constructed from 
Hilbert spaces Ti± as 

Ti = T i+ ®J : i-. (3.210) 

For example, consider matter fields in a black hole spacetime formed by gravitational 
collapse. In this situation, let Hi be a space of all wave packets on the future event 
horizon and 7i 2 be a space of all wave packets on the future null infinity, and 
decompose each Tii into a high frequency part H.i+ and a low frequency part H-i- . 
Typically, we suppose the decomposition at an energy scale of Planck order. If wc 
define Ti± as Fock spaces constructed from Hi±, re specti vely, th en the space T of 
all quantum states of the matter fields is given by ( 3.10E ) with ( 3 . 2 1 Cl| ) . Although 



the following arguments do not depend on the construction of the Hilbert space !F, 
this example should be helpful for us to understand the physical meaning of the 
results obtained. 

For simplicity we consider the case that all T\± have the same finite dimension 
N although in the above example of field theory the dimensions of the Hilbert 
spaces are infinite ^|. The main purpose of this section is to see general properties 
of the entanglement entropy by using a finite system. Anyway, the above example 
of field theory may be helpful in understanding the following arguments. In the 
finite dimensional case we assume the following three physical principles. 

(a) A quantum state \<p) in J 7 is a direct product state given by |<^+)+ <8> \4>-)-, 
where \4>±}± are elements of T± = T\±® 3~2±, respectively. 

(b) 1 </>+)+ is determined by the principle of maximum of the entanglement entropy 
with respect to the decomposition — T\+ ® Ti+ ■ 

(c) A complete measurement of the von Neumann type on the joint system T\ 
is performed by a sender (Alice) in the orthonormal basis {|Vw)i}, each of 
which maximizes the entanglement entropy with respect to the decomposition 
T\ = Ti + <g> T\-. 

In other words the assumption (c) is stated as follows: the state \<fi) is projected by 
one of the basis \ip nm )i- 

In the following arguments, under these assumptions, we show a possibility of 
quantum teleportation of the state \4>-)- in T- to JF 2 : we make a clone of | </>_)- by 
using the quantum entanglement which the state \<p+)+ has. Therefore a receiver 
(Bob), who cannot contact with T\ 1 may be able to get all information about the 
state 1 4* J)- in JF_, provided that he can manage to get the classical channel. 

Now let us show that explicitly. By the assumption (b) and the arguments in 
subsection 3.4.2 (see Eq. ( 3.203[ )), the state | </>+)+ can be written as 



N 

1 



j=J2\n) 1+ ®\n) 2+ , (3.211) 

n— 1 



15 In applying the results for finite dimensions to field theory, we have to introduce a regulariza- 
tion scheme to make the system finite. For example, we can discretize the system by introducing 
a cutoff length. After that, we can consider a finite dimensional subspace of the total Hilbert 
space of the discrctized theory, for example, by restricting total energy to be less than the mass of 
the background geometry. After performing all calculations, we have to confirm that the infinite- 
dimensional limit can be taken. See, for example, the final paragraph of the previous subsection. 
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where (n — 1,2, - ■ ■ ,N) are an orthonormal basis of Ti+. Next, expand 

| </>-)_ as 

|0_)_ = ^C nm |n) 1 _O|m) 2 _, (3.212) 

where {|n)j-} (n = 1, 2, • • • , N) are an orthonormal basis of and J2 nm \Cnm\ 2 = 
1 is understood. To impose the assumption (c), we adopt the following basis 
{\ipnm}i} {n,m — 1, 2, ■ • ■ , iV), each of which maximizes the entanglement entropy. 

1 N 

|Vw)i = -r= e 2 " ijn/N \(j + m)modN) 1+ ® |j>i_. (3.213) 

7 = 1 



In Appendix A7 it is proved that ( 3.213 ) is unique up to a unitary transformation 



in Hence, \4>) — \<fi+)+ ® \<f>-}- is written as 

M = ^T,^™^® U ™\fch, (3-214) 

nm 

where ^2)2 is a state in J- 2 given by 

102)2 = J2 C n'm>\n') 2+ <g> |m')2-. ( 3 - 215 ) 

n'm f 

and Unm {n, m — 1, 2, ■ • • , N) are unitary transformations in J- 2+ defined by 



TV 



Uffl = e- 27rlkn/N \(k + m)modN) 2+2+ (k\. (3.216) 



k=l 



(See Appendix A. 7 for an explicit derivation of ( [3.214 ).) 



Thus, after the measurements in the basis {|VWi)i} by the sender (Alice), the 
original state \<f>) jumps to one of the states \4> nm ) defined by 

\4> nm ) = \i>nm)l <8> U n 2 ^\4> 2 ) 2 . (3.217) 

This state can be seen by the receiver (Bob), who cannot contact with T\, as the 

(2-t-l ~ (2 I \ 

state Unm 1 02) 2 in Fi- Here note that the unitary transformation U nm in J- 2 + is 
completely determined by a pair of integers n and m (outcome of the experiment 
by Alice). Thus, if the two integers are sent to the receiver (Bob) in the classical 
channel, then by operating the inverse transformation of the corresponding unitary 
transformation in T 2 + the receiver (Bob) can obtain the 'clone' state |0 2 )2 (G F%) 
of \4>-) (e •?""-)• It is evident that |0 2 )2 has all information about the original state 
10-)- 

It is remarkable that information to be sent to the receiver (Bob) in the classical 
channel is only two integers n and m, while information included in the unknown 
state |0-)- is a set of complex constants {C nm } (n,m — 1, 2, • • • , N) with a con- 
straint ^2 nm \ C nm \ 2 = 1. Thus a large amount of information is sent in the quantum 

channel. Here we mention that tracing out T 2 + from the state U nm ^2)2 or \(f> 2 } 2 
results in the following density matrix p 2 ^ on T 2 -\ 

Pi- =E E C ^« N2-2-H, (3-218) 



which is equivalent to the density matrix obtained by tracing out T\- from the 
original unknown state |0_)_. Hence, if the receiver (Bob) cannot contact with 
J-2+, he does not obtain any information from the sender (Alice). 
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Finally it must be mentioned that the success of quantum teleportation is due to 
quantum entanglement in the state |0+) + which has maximal entanglement entropy. 
If we took with less entanglement entropy then the teleportation would be 

less successful. Therefore, the entanglement entropy can be regarded as an index of 
efficiency of quantum teleportation. This consideration supports the interpretation 



of the entanglement entropy proposed in subsection 3.4.1 



3.4.4 Concluding remark and physical implications 

In this section a new interpretation of entanglement entropy has been proposed 
based on its relation to the so-called conditional entropy and a well-known meaning 
of the latter. It is conjectured that entanglement entropy of a pure state with 
respect to a division of a Hilbert space into two subspaces 1 and 2 is an amount of 
information, which can be transmitted through 1 and 2 from a system interacting 
with 1 to another system interacting with 2. The medium of the transmission is 
quantum entanglement between 1 and 2. 

To support the interpretation we have given the following two suggestive argu- 
ments: variational principles in entanglement thermodynamics and quantum tele- 
portation. The most important variational principle we considered is the principle 
of maximum of entanglement entropy. This principle determines a state uniquely 
up to a unitary transformation in one of the two Hilbert subspaces (not in the 
whole Hilbert space). From the proposed conjecture it is expected that information 
can be transmitted most effectively through the two subspaces by using the max- 
imal entanglement of the state. In fact, reformulating the quantum teleportation 
in terms of the entanglement entropy, we have shown that the state having maxi- 
mal entanglement entropy plays an important role in quantum teleportation. This 
consideration gives strong support to our interpretation. 

As a by-product we have shown that the variational principle of minimum of 
entanglement free energy is useful to determine a quantum state. The resulting 
quantum state has exactly the same form as those appearing in the thermo field 
dynamics of black holes |g5| and the quantum field theory on a collapsing black 
hole background [B3| , provided that the entanglement temperature T ent is set to be 



the black hole temperature. It is remarkable that, as shown in section B.3, T ent for 
a real massless scalar field in Schwarzschild and Reissner-Nordstrom spacetimes is 
equal to the black hole temperature of the background geometry up to a numerical 
constant of order 1. Thus we can say that the variational principle of minimum of 
entanglement free energy gives a new derivation of the Hawking radiation. Finally, 
we mention that with this variational principle the entanglement thermodynamics 



is equivalent to 'tHooft's brick wall model 1 35 . 

It will be valuable to analyze how to generalize arguments in this section to 
the situation that divergences in entropy and energy are absorbed by renormaliza- 
tion |7(| ^| . If the generalization is achieved, the physical meaning of the entan- 
glement entropy in black hole physics will become clearer. It is noteworthy that in 



the brick wall model, as shown in section 3.2, the divergence in thermal energy is 
exactly canceled by divergence in negative energy |Q . 

Now the final comment is in order. It is worthwhile to clarify in what phys- 
ical situations the variational principles can be applicable. (In thermodynamics 
the second law supports the principle of maximum entropy.) In other words, in 
what situations does the entanglement entropy increase? In what situations does 
the entanglement free energy decrease? To answer these questions, theorem u\ in 



subsection 2.3.2 or its generalization may be useful. 



Chapter 4 

Discussions 



In this thesis we have analyzed properties and the origin of the black hole entropy 
in detail from various points of view. 

First, in chapter |^ laws of black hole thermodynamics have been reviewed. In 
particular, the first and generalized second laws have been investigated in detail. It 
is in these laws that the black hole entropy plays key roles. 

In section [2~T| wc have re-analyzed Wald and Iyer-Wald derivation of the first law 
of black hole mechanics in a general covariant theory of gravity, following Ref . ^(J . 
In particular, two issues listed in the beginning of section 2.1 have been discussed 
in detail: (a) gauge conditions and (b) near-stationary black hole entropy. It has 
been shown that the minimal set of gauge conditions necessary for the derivation of 
the first law for stationary black holes is that t a and (p a are fixed at spatial infinity, 
where t a is the stationary Killing field with unit norm at infinity, and ip a denotes 
axial Killing fields. It has also been shown that for non-stationary perturbations 
about a stationary perturbation the first law does hold to first order in perturbation. 

However, this first law cannot be applied to a purely dynamical situation. In this 
sence we have called it the first law of black hole statics. The purpose of section |2.2| 
has been to consider dynamical definition of black hole entropy and to derive a 
dynamical version of the first law of black hole, which we call the quasi-local first 
law of black hole dynamics. For simplicity, we have considered the general relativity 
only. Extension t o a g eneral covariant theory of gravity will be valuable. 

In subsection 2.2.1 we have considered two non-statistical definitions of entropy 
for dynamic (non-stationary) black holes in spherical symmetry. The first is analo- 
gous to the original Clausius definition of thermodynamic entropy: there is a first 
law containing an energy-supply term which equals surface gravity times a total dif- 
ferential. The second is Wald's Noether-charge method, adapted to dynamic black 
holes by using the Kodama flow. It has been shown that both definitions give the 
same answer for Einstein gravity: one-quarter the area of the trapping horizon fill . 

In subsection 2.2.1 the quasi-local first law of black hole dynamics has been 
derived without assuming any symmetry and any asymptotic condition p2[ . In the 
derivation we have given a new definition of dynamical surface gravity. In spherical 
symmetry it reduces to that defined in Ref. j55|. 

In section 2.3 we have proved the generalized second law for a quasi-stationary 
black hole which is formed by gravitational collapse |u|. After that, in subsec- 
tion 2.3.3, we have discussed a generalization of our proof to a dynamical back- 



ground. It has been suggested that the generalization may be possible by using the 



quasi-local first law derived in subsection 2.2.1. 

Next, in chapter ^ three candidates for the origin of the black hole entropy have 
been analyzed: the D-brane statistical-mechanics, the brick wall model, and the 
entanglement thermodynamics. 
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In section 3.1 the D-brane statistical- mechanics has been reviewed by using a 
configuration of D-strings and D-fivebranes wrapped on T 5 = T 4 x S 1 , which was 
introduced in Ref. [43] . We have consider a set of multiply- wound D-strings, which 



is composed of Nql' D-strings of length 2irRqi (q± = 1,2, ■■■) along the S 1 and 

(5) 

a set of multiply- wound D-fivebranes, which is composed of N q5 D-fivebranes of 
length 2-?rRq 5 (q 5 = 1,2, ■■■) along the S 1 . Here {N$} and {N$} are sets of 
arbitrary non-negative integers, and R is radius of the S . It has been shown that 
the number of microscopic states of open strings on the D-branes is bounded from 
above by exponential of the Bekenstein-Hawking entropy of the corresponding black 
hole, and that the temperature of a decay of D-brane excitations to closed strings 
is bounded from below by the Hawking temperature of the corresponding black 



hole. This result has been summarized as Eqs.( 3.31 ) and ( 3.37 ). The necessary and 
sufficient condition for these bounds to be saturated has been shown explicitly and 
some speculations has been given. 



In section 3.2 we have re-examined the brick wall model to solve problems con- 
cerning this model. In particular, it has been shown that the wall contribution to 
the total gravitational mass is zero if and only if temperature of thermal gas mea- 
sured at infinity is set to be the Hawking temperature, and that the backreaction 
can be neglected Q. 



In section 3.3 we have constructed entanglement thermodynamics for a massless 
scalar field in Minkowski p5| , Schwarzschild |l6| and Reissner-Nordstrom space- 
times. The entanglement thermodynamics in Minkowski spacetimc differs signifi- 
cantly from black-hole thermodynamics. On the contrary, the entanglement ther- 
modynamics in Schwarzschild and Reissner-Nordstrom spacetimes has the same 
structure as that of black-hole thermodynamics. In particular, it has been shown 
that entanglement temperature in the Reissner-Nordstrom spacetime approaches 
zero in the extremal limit. 



In section 1A a new interpretation of entanglement entropy has been proposed 
based on its relation to the so-called conditional entropy and a well-known meaning 
of the latter [[l7). It has been conjectured that entanglement entropy of a pure 
state with respect to a division of a Hilbert space into two subspaces 1 and 2 is an 
amount of information, which can be transmitted through 1 and 2 from a system 
interacting with I to another system interacting with 2. The medium of the trans- 
mission is quantum entanglement between 1 and 2. To support the interpretation 
we have given the following two suggestive arguments: variational principles in en- 
tanglement thermodynamics and quantum teleportation. It has been shown that 
the state having maximal entanglement entropy plays an important role in quantum 
teleportation. This consideration gives strong support to our interpretation. 

Now let us discuss about semiclassical consistencies of the brick wall model and 
the entanglement thermodynamics. 

It has been shown that the brick wall model is a consistent semiclassical descrip- 
tion of black hole entropy: thermal excitations raised to the Hawking temperature 
above the Boulware state explains the black hole entropy; the positive divergence 
in thermal energy is canceled by the negative divergence in the vacuum energy of 
the Boulware state. Namely, the following simultaneous equations have a solution: 

Swall — SbH, 

Too — Tbh, 

(AM) therm , wa u + (AM) B , wa ii = 0. (4.1) 



It has been found that the last equation is equivalent to the second one. Thus, the 
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solution is 



Too = T BH . (4.2) 

On the other hand, entanglement thermodynamics has the following properties 
common to all definitions of entanglement energy. 



S, 



E ent ~ iVA^^, 

T ent ~ AA T r BH , (4.3) 

where A/s and A/e are numerical factors of order unity, and Mr = 2itMe /Ms- 
Here, we multiplied entanglement entropy 5 e nt and energy £? e „t by the number 
./V of fields Q These properties are expected to hold also for states different from 
the Boulware state, provided that the definition of entanglement quantities are 
properly modified. For other states, of course, the numerical factors As and Me 
will be different from those for the Boulware state. 

Our question now is whether the entanglement thermodynamics is a consistent 
semiclassical description of black hole thermodynamics or not. Here let us assume 
that the sum of the vacuum energy and the entanglement energy contributes to 
gravitational energy. Hence, the question is restated as whether the following si- 
multaneous equations hold or not. 

Sent — £>BHj 

T en t = Tbh , 

Eent + [AM )B, wall = 

0. (4.4) 



section 3.3 



Since t he c utoff length a in section |3.2| seems to be related to the cutoff length a in 
as 

a ~ nso, (4-5) 



a solution of Eqs. (4.4) is given by 



Me =s -A 




UL! " ^ ^-210AV (46) 

Unfortunately, since Ms = 0.3, the r.h.s of the last equation is less then unity. Hence 
the last equation does not hold for positive integer value of ns- Even if ng would be 
allowed to be less than unity, the value of ns given by the last equation would not be 
consistent with the second equation for all our definitions of entanglement energy 0. 
This discrepancy suggests that our model of entanglement thermodynamics suffers 
from a strong backreaction near horizon. This might mean that our choice of the 
pure state (the Boulware state) would be wrong. However, we strongly expect that 



there is a state for which the consistency condition (4.6) holds. It will be worthwhile 



to analyze whether Eqs. ([!£]) are satisfied for the Hartle-Hawking state or not. 

Note that N = 1 corresponds to one bosonic field. Hence, M = Nir 4 /90. 
2 Note that Me depends on rig . 
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Next, we shall combine a variational principle in entanglement thermodynamics 
with the third equation in the semiclassical consistency conditions. As explained in 
subsection 3.4.2, the principle of maximum of entanglement entropy introduced in 
Ref. p7[ determines a quantum state to be (3.203) up to a unitary transformation in 
one of two Hilbert subspaces. Moreover, the principle of minimum of entanglement 
free energy determines a quantum state to be ( 3.20E ) up to a unitary transforma- 
tion in one of two subspaces, too. It has also been mentioned in the second-to-last 
paragraph of subsection 3.4.2 that, if we m aximize entanglement entropy with en- 
tanglement energy fixed, the state ( 3.209 ) is obtained. In this case, entanglement 
temperature should be determined so that the entanglement energy coincides with 
the fixed value. On the other hand, the third equation in the semiclassical consis- 
tency conditions (4.4) does fix entanglement energy, provided that we do not change 
(AM) b .wail- We shall call this condition the small backreaction condition. Hence, 
the principle of maximum of entanglement entropy combined with the small back- 
reaction condition determines a state of a quantum field near horizon to be (3.209) 
up to a unitary transformation in one of two Hilbert subspaces, provided that a 
suitable regularization scheme is introduced. This state has exactly the same form 
as those appearing in the thermo field dynamics of black holes |85f| and the quantum 
field theory on a collapsing black hole background provided that the entan- 
glement temperature T ent can be set to be the Hawking temperature. In this case, 
since the entanglement energy is the same as thermal energy in the brick wall model 
because of the small backreaction condition, and since the entanglement entropy is 
expected to have the almost same value as the Bekenstein-Hawking entropy, the 
entanglement temperature also seems to coincide with the Hawking temperature 
up to a numerical factor of order unity. Therefore, it seems that the principle of 
maximum of entanglement entropy combined with the small backreaction condition 
may provide a new derivation of the Hawking radiation. Further investigations on 
this point will be valuable. 

Finally, let us discuss the information loss problem. 

We have proposed a new interpretation of entanglement entropy: entangle- 
ment entropy of a pure state with respect to a division of a Hilbert space into 
two subspaces 1 and 2 is an amount of information, which can be transmitted 
through 1 and 2 from a system interacting with 1 to anoth er s ystem interacting 
with 2. On the other hand, it has been confirmed in section 3J3 and in many ref- 
erences [^7[ H^, [?^, that the entanglement entropy has the same value as the 
black hole entropy up to a numerical constant of order unity, provided that a cutoff 
length of Planck order is introduced in the theory. Hence we have a large amount 
of entanglement entropy to transmit information from inside to outside of a black 
hole by using quantum entanglement. 

Hence, in our interpretation, it seems that the entanglement entropy is a quantity 
which cancels the black hole entropy to restore information loss, provided that the 
black hole entropy represents the amount of the information loss. For example, 
suppose that a black hole is formed from an initial state with zero entropy (5 = 0). 
In this case, non-zero black hole entropy is generated (Sbh > 0) from the zero 
entropy state. At the same time, entanglement entropy and negative conditional 
entropy are also generated and their absolute values are as large as the black hole 
entropy (S en t = l-Scondl — Sbh)- After that, the black hole evolves by emitting 
Hawking radiation, changing the value of Sbh and S en t (= |5cond|) with Sbh — S en t 
kept. Finally, when the black hole evaporates, the entanglement entropy cancels the 
black hole entropy to settle the final entropy to be zero (S — 0). To summarize, 
the black hole entropy is an amount of temporarily missing information and the 
entanglement entropy is a quantity which cancels the black hole entropy. Both 
entropies appear and disappear together from the sea of zero entropy state. 

We conclude this thesis by giving some speculations. 
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(i) The black hole entropy should be related to a number of microscopic states in 
the corresponding black hole. The microscopic description should be possible 
by using a quantum theory of gravity (eg. superstring theory, loop quantum 
gravity, etc.). 

(ii) The existence of the horizon prevents the microscopic states in the black hole 
from being seen on the outside. Hence, information about the microscopic 
states is lost at least temporarily. (See Ref. [^H) for Hawking's objection 
to string-theorist's point of view.) At the same time, the horizon generates 
entanglement entropy of matter fields which has the same value as the black 
hole entropy. 

(iii) The entanglement entropy is an amount of information, which can be trans- 
mitted through the matter fields from a system interacting with the matter 
fields in the black hole to another system interacting with the matter fields 
outside the black hole, provided that the classical channel in the sense of 
the quantum teleportation can be transmitted properly. The former system 
carries the temporarily missing information corresponding to the black hole 
entropy. Hence the quantum channel of all temporarily missing information 
about the microscopic states can be transmitted to the outside of the black 
hole before the black hole evaporates completely. It is the classical chan- 
nel that is necessary for restoring the original information from the quantum 
channel. 

(iv) We propose the conjecture that one of the following two should be realized. 

(a) The black hole evaporates completely, and the classical channel is trans- 
mitted at the final stage of the black hole evaporation. 

(b) There remains a remnant at the end of the Hawking radiation, and all 
or a part of the classical channel is carried by the remnant forever. 

If (a) is correct, then it is in principle possible to restore all information about 
the microscopic states of the black hole after the evaporation. If (b) is correct, 
then all or a part of the information about the microscopic states cannot be 
restored although all the information remain to exist: the information cannot 
be decoded since the classical and quantum channels are located separately. 
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Appendix 



A.l The conditional probability 

In this appendix we reduce ( |2.62| ) to ( |2.64| ) . First the S- matrix obtained by |^37| is 



71=0 



2 n n\ 



where e and N are a bivector and a normalization constant defined by 
e = 2 Xj(j A <g> jT) S ym , N = Y[ \/l - x~i, 



(A.7) 



where 



Xi = exp {—ir(u>i - VlBHmi)/n) . 



In the expression, u>i and m, arc a frequency and an azimuthal angular momentum 
quantum number of a mode specified by integer i, Qbh and n are an angular 
velocity and a surface gravity of the black hole. ci, ip, i<j, i\ and^r are unit vectors 
in Ttz+ © Hh+ defined in 67 , and the former four are related as follows: 



a a = Ti iO a + Ri ip a , 



(A.- 



where ti, Ti are transmission coefficients for the mode specified by the integer i on 
the Schwarzschild metric [b7| and r%, Ri are reflection coefficients. They satisfy 3 



\U\ 2 + \n\ 2 

t- — T 

By using the S-matrix, we obtain 



mi 



- m 

-R*TJT* 



1. 



(A.9) 



S|{n i7 }) 



N 



TT-i=[i? 4 a t ( 4 p)+T l at( l( T)]"'-' 

n=0 



E 



2 n n\ 



(5e) 



n 



3 The last two equations are consequences of the time reflection symmetry of the Schwarzschild 
metric. 
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/ (2n + ^n, 7 )! 
{2n)\ 



, r m li+rrn rii-li+n^-rrn 



^E E Ex 

n<=Om< =0 2^=0 \ 



n 





1 








/ n i7 ! 







r>mirri" 

Hi 1 s 



, r ni h+rm ni-h+n.^-mt 



(A.10) 



where denotes a summation with respect to rij, m, and ij over the following 
range: Yli rii — n, rii > 0, < rrii < n n , < U < rii. In (A. 10), those orthonormal 
basis vectors in |{7i iP }) <8> T(TLh+) that have a non- vanishing inner product with 
5|{n i7 }) appear in the formQ 



/ (E t (2n t + n i7 ))! 



n 



(A.ll) 



sym 



Thus, when we calculate (({« ( p}| <8 S1{n i7 }) , the summation in ( |A.10| ) is re- 
duced to a summation with respect to n% and over the range n$ > max(0,rt iP — 
n i7 ), max(0, n iP — rii) < fn% < min(n iP , n i7 ) with ij = rt iP — m, ^. Here is an 
element of J 7 (Hh+)- Paying attention to this fact, we can obtain the following 
expression of the conditional probability. 



P(M\{n n }) 

= |Af E 

{ni>max(0,n. p -n. 7 )} 

min(n. p ,n. - f ) 



E( 2 ^+^-v)j ! 



n 



E 



«7 

r„ - rrii 



2t 



(E l (2n, + n i7 ))! 



n 



]] Ili [niln p l(n t + n i7 - ti iP )! 

n 

i 

The inner product in the last expression equals to [] 



iP h i<j 



(A.12) 



' IL [rjjHpKnj + n n - n iP )!] 
(Ei(2«i + ^))' 



4 The number of the 'particle' iO in ( A,lo| ) is nj + « i7 -n iP , setting the number of the 'particle' 
iP to n iP . 

5 The range is obtained by inequalities n; > 0, < rrii < ra i7 , < h < rii, + mj = n iP and 
- 1 - ";T — "iP > 0. 

6 (A.ll) is normalized to have unit norm. 
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Finally, by using (A. 9) and exchanging the order of the summation suitably, we can 
obtain 



|2\™i7+";p 



P({n iP }\{n i7 }) 

=n[( i -^K>(i-i^i 2 )' 

i 

n*n^,)u>Ww,) [-1^17(1 -\R i \^)] k ^r kl \n iP \ 



E E 



L=0 



E 



ni=n. p — min(/i ,mi) 

This is what wc had to show 



-± Q kl(n a - k)l(n iP - li)\mi\(n zl - m,i)\(n zP - to;)! 

nj\(ni - n zP + n n )\ 
(rii - n iP + k)\(ni - n iP + m*)! 



2 I d 1 2 \nv — n. , 



A. 2 A proof of Lemma 2 

In this appendix we give a proof of Lemma ^. 
Proof 

Since a set of all,T and i<7 generates TLh+ @]> the definition of T^^^^y leads 



T {nJ,}{nf p } = E ({ n >T,n P ,^ a }\S\{n n })({ri^}\S\{n iT ,rt p ,n ia }), (A.13) 

{"i^} ! {n i x} 



where 



\{n iT ,n. p ,n. a }) = ]J 



(at( iT )) n - (aW* (aV))"" 



|0). 



In the expression, S'| {n. 7 }) is given by ( |A.10|) and S\{n' 7 }} is obtained by replac- 
ing n a with n' in ( |A.1Q ). Now, those orthonormal basis vectors of the form 
|{n. r , n iP , n i(T }) that have a non-zero inner product with S , |{« i7 }) must also be of 



{n. }{n' } 

the form (A.ll). Thus, T^y^J; 1 ^ vanishes unless there exist such a set of integers 



{rii, n,-} (i = 1, 2, • • •) that 



m + n,~ — n A . 



(A. 14) 



for Vi. The existence of {n-i} and {r^} is equivalent to the condition 



n iP — n' p for Vt. 

□ 



A. 3 On-shell brick wall model 

When we performed the diff eren tiation with respect to /3oo to obtain the total energy 
and the entropy in section ^2, the surface gravity Ko of the black hole and the 
inverse temperature [3^ of gas on the black hole background were considered as 
independent quantities. Since in equilibrium these quantities are related by (3^ = 
kq/2tt, we have imposed this relation, which wc call the on-shell condition, after the 
differentiation. In fact, we have shown that the wall contribution to gravitational 
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energy is zero and the backreaction can be neglected, if and only if the on-shell 
condition is satisfied. 

On the other hand, in the so-called on-shell method |^2], [75|, the on-shell 
condition is implemented before the differentiation. Now let us investigate what 
we might call an on-shell brick wall model. With the on-shell condition, the wall 
contribution to the free energy of the scalar field considered in subsection |3.2.3 is 
calculated as 



p(on-shell) _ A j3^J I . 

wal1 ~ 4 36(W 1 °> 

If we define total energy and entropy in the on-shell method by 



j r(on— shell) 9 f ^(on— shell) 

U wall — qo yP^-^wall 

-7(071 — shell) ^2 ^ j^(on — shell) 

^00 



^wall — Poo qB ^wall ' (A.10J 



then these quantities can be calculated as 

jj{on — shell) ~ 

U wall = U ' 

q(on-shell) _ A I 1 _ I 

where S wa ii is the wall contribution ( 3.62| ) to entropy of the scalar field with = 
Tbh- 

It is notable that the total energy U^™ l ^ shell ' > in the on-shell method is zero 
irrespective of the value of the cutoff a. However, S^ l [ shell ' > is always smaller 
than S wa ii- It is because some physical degrees of freedom are frozen by imposing 
the on-shell condition before the differentiation. Thus, we might miss the physical 
degrees of freedom in the on-shell method. 



A. 4 Symmetric property of the entanglement en- 
tropy for a pure state 

In this appendix we first give an abstract expression for the reduced density oper- 
ators p\ and P2 corresponding to a pure state u in T = J-\®T2, which do not use 
the subtrace. Then with the help of them we prove that S en t obtained from p\ and 
P2 coincide with each other. We follow the notations in §3.3.1. 

Proposition 1 For an arbitrary element u of T = T\®J-2, there are antilinear 
bounded operators A u S B{!F\, JF 2 ) and A* u e B(J-2,J-i) such that 

(A u x, y) = (A* u y, x) = {u,x® y) (A.18) 

for v x G T\ and v y € T2 ■ 

Proof. Fix an arbitrary element x of T\. Then {u, x®y) gives a linear bounded 
functional of y(£ J- 2) since 

\(u,x®y)\ < ||u||||x||||t/||. 
Hence by Riesz's theorem there is a unique element z UlX of T2 such that 

(z u ,x, y) = (u, x ® y) (A.19) 
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for v y £ T%. Let us define ^4 U by A u : x 
bounded operator from T\ to J~i since 

\\A-vA\ = ll^xll = IMIIMI- 

Exchanging the roles played by T\ and T% in the above argument, it is shown 
that there is an antilinear bounded operator A* u from T% to T\ such that (A* u y, x) = 
(u,x®y). □ 

Note that A u and A* u defined above are written as 

A u x = y^fj(x®fj,u) , 

3 

A uV = ^ei{ei^y,u) . (A.20) 

i 

Using this expression, it is easily shown that 

A* u A u x = ^ei{ei® fj,{u,x® fj)u) , 

A u A* u y = 5^/j(e i ®/ J -,(«,e < ®y)«) . (A.21) 



These coincide with p\ and pi, respectively, if u has unit norm (see Eq. ( 3.110 ) and 
( 3.111 )). Therefore the following proposition says that p\ and p2 have the same 
spectrum and the same multiplicity and that entropy of them are identical. 

Proposition 2 p^ (g B(T{)) and pi^ (€ B{J 7 2)) defined by 

n W - A* A 

Pl 2) = A U A* U (A.22) 



are non-negative, trace-class self-adjoint operators, where A u and A* u are defined in 
Proposition R for an arbitrary element u of T . The spectrum and the multiplicity 
of p« and are identical for all non-zero eigenvalues. 

Proof. In general 

{x',p[px) = (A u x, A u x') 
for y x,x' (<G T\) by definition. Therefore 

(x,p^x) = \\A u x\\ 2 >0 (A.23) 



and 



i 

= ^^{A u ei, fj){fj, A u ei) 

i,3 
hj 

= INI 2 , (A.24) 



i.e. piP is non-negative and trace-class. In general a non-negative operator is 
self-adjoint and a trace-class operator is compact. Thus the eigenvalue expansion 
theorem for a self-adjoint compact operator says that all eigenvalues of pu are 
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discrete except zero and have finite multiplicity. For a later convenience let us 
denote the non-zero eigenvalues and the corresponding eigenspaces as Xi and T\j, 
= 1,2,-..). 

(2) 

Similarly, it is shown that p„ is non-negative and trace-class and that all eigen- 
values of it are discrete except zero and have finite multiplicity. 

Now kerp\?> = ker A; since p^y = A u A* u y and {y,Pu ] y) = \\A* u y\f for an 



arbitrary element y of Ti by definitions. Moreover, from (A. 18) it is evident that 

y _L RanAu <^ y e ker A*. 

With the help of these two facts T% is decomposed as 

T 2 = ker p^> © Ran A,. (A.25) 

where the overline means to take a closure. 

Moreover, it is easily shown by definitions that 

p\i ^ A u x — \^A u x , 

(A u x,A u x') = \i(x',x) (A.26) 

for v x (e and v x' (<E Fi^')- Hence A u maps the eigenspace to a eigenspace 

of pu with the same eigenvalue, preserving its dimension. Taking account of 



Eq.(A.25), this implies that the spectrum and the multiplicity of pu and p^ are 



identical for all non-zero eigenvalues. □ 

A. 5 Entanglement energy for the case of B = R 2 
in Minkowski spacetime 

First we take B — R 2 . Without loss of generality the resulting two half-spaces are 
represented as {(xi, x%, x%) : x\ > 0} and {(xi, X2, £3) : x\ < 0} p7[ ]. 

Here some comments are in order. Since all the degrees of freedom on and 
across B, which is infinite, contribute to the entanglement energy, a suitable cut- 
off length a(> 0) should be introduced to avoid the ultra-violet divergence. For 
the same reason, the infra-red divergence is also anticipated in advance, since B is 
non-compact in this model. The latter is taken care of by considering the massive 
case since the inverse of the mass characterizes a typical size of the spreading of the 
field. Clearly a should be taken short enough in the unit of the Compton length of 
the field, m^ 1 , to obtain meaningful results. Therefore we shall only pay attention 
to the leading order in the limit m^a — ► in the course of calculation as well as in 
the final results. These remarks are valid in any model of this type, and the same 
remarks apply to the case of the entanglement entropy, too ^8| . 

In order to calculate (: Bitot -)p' for the present case, we first note that the 

term VABQ A q B in Eq.( 3.114 ) corresponds to the expression J (V(/>) 2 + m 2 ^ 2 d 3 x 
read off from Eq.( 3.156| ), which defines an op erator V(x, y) acting on a space W 



({<!>{•)}, d 3 x). In order to use the formula ( 3.149 ), thus, we need the inverse of 
the positive square-root of aV . For this purpose, it is convenient to work in the 
momentum representation of W p7[ given by 



= I <p£exp[ik ■ x], 



d 3 k 
(2tt) 3 

d 3 x cj)(x) exp[-ik ■ x] . (A.27) 
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The results are 

/d 3 k -> 
^-^■(fc 2 + m|)exp[ifc • {x-y)] , 

W~\x,y) = / -^(P+ml)-^ ex p[ik-(S-y)}. (A.28) 

Note that both V(x,y) and y) are symmetric under the exchange of x and 

y. (The cut-off must preserve this property.) Now the formula (3.14E) gives 

(:H tot :) p , = -~ / d 3 y / d 3 z / -^(fc 2 + m 2 ) exp[ifc • (y - f)] 

2J yi< - a J Xl >a JlfciKa-i (, 27r J 



rf 3 fc' 



(fc 2 + m 2 ,)^ 72 exp[ifc' • (x - y)], (A.29) 

im<o- 



where, as discussed above, a cut-off length a was introduced in the integral. 

Since the integrand is invariant under the translation along B, the integral with 
respect to X2 and x$ yields a divergent factor A = J R2 dx2dx$. Clearly this factor 
should be interpreted as the area of B. If this divergent integral A is factored out, 
we obtain the following convergent expression for (: H to t '-)p r - 

, rr \ ilp. f . M 2fc !l dfcl A""' 

(: ±f tot = -— / dyi / ctai 



x (^, 2 + fc! 2 + m 2 )(/c|| 2 + fc'x 2 + m 2 )- 1 / 2 
xexp^fcx-fc'xXz/i-zx)]. (A.30) 

Here fen is a 2-vector lying along f? and k\, k\ are components normal to B (if 
we make an obvious identification of R 3 with its Fourier space). Let us change the 
variables from x\ and y± to z = xi — yi and u = (xi + y\)/2. Then z and w take 
values in the range z < 2a and — (| — a) > u > (| — a), respectively. Hence the 
integration with respect to u yields 

(:H tot :) p , = -- / dz(z-2a) 1 1 ' ' 1 



2 y 2a V (2tt) 2 J_ a -i 2n 2tt 

x (fcfj + fci 2 + m|)(fc 2 + k\ 2 + ml)- 1 ' 2 eatpH(Ai - k\)z] 



— / azyz — 2a) / cos(fciz) / — re(re +Ki ) 

2 J2a J-a- 1 27T J m 27T 



dk'i 
~2i7 



(K 2 + fc'i 2 )- 1/2 cos(fc'iz) (A.31) 



in the leading order, where k is defined by k 2 — k 2 + m|. Here note that in this 
expression, the integration with respect to k\ followed by that with respect to n 
leads to an infra-red divergence if we set m = 0, in accordance with our discussion 
at the beginning of this subsection. 

Now let us recollect some formulas with the modified Bcsscl functions |9ql : 



i cos t 

K Q (x) = / dt- 



L 



dx xKq(x) = xqK\(xq) , 



Xo 

dx x 3 Kq(x) = Xq K 1 (x ) + 2x 2 K 2 {xo). (A.32) 
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With the help of these formulas E { L t is written as 



E ent = — 7T / dz(z-2a) / — cos(fciz) 



1 

X 2^ 



-(/ci 2 + Tr?^)K\ (vnz) + 2-^K 2 (mz) 



Air 3 a 3 

in the leading order. Here we have introduced 



'ai (ma) + 02 (ma) + 03 (m«)] (A. 33) 



ai{x) = -x J ^^Xi«)[2^cose+(e 2 -2)sine], 

f°° £ — 2 
a 2 (x) = -2a; 2 j ^3- A" 2 (asO sin£, 

/•OO to 

a 3 (x) = -x 3 ^ de^-tfifcO sin ^ (A.34) 



A numerical evaluation shows 

[ai(cc) + 0:2(21) + 0:3(2;)] ~ 0.05 as x — > 0. 

Therefore we get [] 



„ > 0.05A 
H *ot :>p' « W . (A.35) 



A. 6 States determined by variational principles 



In this appendix we give derivations of formulas ( 3.203| ) and fl3.209p . 
We consider a Hilbert space T of the form 

T = T\® (A.36) 

where T\ and Ti are Hilbert spaces with the same finite dimension N. An arbitrary 
unit clement \<fi) of T is decomposed as 

N N 

\<t>) = J2 °™\n}i ® |m) 2 , (A.37) 

n=l m=\ 

where |n)i and |n) 2 (?i = 1, 2, • • • , AT) are orthonormal bases of .Fi and ^2, respec- 
tively, and J2 n m \C nm \ 2 = 1 is understood. Here, without loss of generality, we 
can choose the orthonormal basis of T 2 to be eigenstates of the normal-ordered 
Hamiltonian : H 2 : of the sub-system as 

: H 2 : |n} 2 - E n \n) 2 . (A.38) 

Since C^C is a non- negative hermitian matrix, we can define a set of non-negative 
real numbers {p n }, each of which is the eigenvalue of the matrix C^C. Hence, 

CtC = V*PV, (A.39) 

7 If we adopt another regularization scheme with the same cut-off length a, the result may 
change. However, the change is in sub-leading order. 



Appendix 



101 



where V is a unitary matrix and P is a diagonal matrix with diagonal elements 
{p n }- With these definitions, the entanglement entropy S en t and the entanglement 
free energy F ent are calculated as 



JV 



Sent = ~ y^Pnlnp», 

n=l 
N N 

Pent — ^ y ^ ^ -^nPm | ^mn | ~t~ Pent ^ ^ 
n— 1 m— 1 n— 1 

The constraints ^ n m |C«™| 2 = 1 an( l = 1 are equivalent to 



(A.40) 
(A.41) 



E p » = 1 ) 

n=l 

N 



(A.42) 



i=i 



Thus, the variational principles are restated as follows: to maxim ize ( A.4C) under 
the constraints ( A.42 ); to minimize ( A.41 ) unde r the constr aints ( A.42 ). 

Now, we shall show that expressions ( A.40 ) and ( A.41 ) are same as those for 
entropy and free energy in statistical mechanics in the subspace Ti- Let us consider 
a density operator p on Ti'- 



N N 

P = ^ E P nm\n)22(m\ 



(A.43) 



n— 1 m— 1 



where P nm is an N x N non-negative hermitian matrix with unit trace. By diago- 
nalizing the matrix P as 

P = V^PV, (A.44) 
we obtain the following expressions for entropy and free energy. 

N 



S = 



71=1 

N N 



N 



f = E n p m \v mn \ 2 +Ty^p n \np n , 



(A.45) 



where {p n } are diagonal elements of the matrix P and T is temperature. The 
constraints Trp ~ 1 and V'V = 1 are restated as 



JV 



n=l 

JV 



(A.46) 



i=i 



At this point, it is evident that the variational principles of maximum of entropy 
are the same in entanglement thermodynamics and statistical mechanics and that 
the principles of minimum of free energy are also the same in the two schemes. 
Hence, the principle of maximum of the entanglement entropy gives 



(A.47) 
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as the principle of maximum of entropy gives the microcanonical ensemble 



P =-6 

x nm — u nm 



(A.48) 



in statistical mechanics. Similarly, the principle of minimum of the entanglement 
free energy gives 

(&C) nm =Z- l e- B ^ T ^6 nm , (A.49) 
as the principle of minimum of the free energy gives the canonical ensemble 



P - Z~ 1 e~- B » /T J 

L nm — ^ u nm 



(A.50) 



in statistica l mec hanics , whe re Z = e En ' T "'' and Z = e E ™/ T . It is easy 



to see that ( |A.47| ) and ( |A.49| ) are equivalent to ( |3.203| ) and ( |3.209| ), respectively, up 
to a unitary transformation in T\. 

Finally we comment on the generalization of the analysis when the Hilbert space 
is divided into two subspaces with different dimensions (dimT\ > dimTz). In this 
case, by defining S en t from p2, we obtain similar results. 



A. 7 Bell states 



In this appendix we show that in the finite dimensional case the orthonormal basis 
{|Vw)i} defined in the physical principle (c) is given by ( |3.213 ) uniquely up to a 
unitary transformation in !F\+. After that, we derive the equation (3.214). 
We consider the following decomposition of the Hilbert space T\ . 



(A.51) 



where T\+ and T\- are Hilbert spaces with the same finite dimension N. From the 
arguments in Appendix A. 6, each of the basis {|V>nm)i} is obtained by applying a 
unitary transformation in T\+ to the following state in T\. 



1 N 

-717 £ M 

1=1 



1+ 



® b')i-) 



(A.52) 



where and (j = 1,2, •••,N) ar e orth onormal bases of T\ + and .Fi_, 

respectively. Evidently, any states given by ( 3.213 ) are obtained by this p rocedu re. 
Moreover, it is easily confirmed as follows that a set of all states given by ( 3.213 ) is 
a complete orthonormal basis in the N x N dimensional Hilbert space T\ ■ 



1 N 

±_^2 e ^ij(n'-n)/N Sjn 



nn' u mm' 



(A.53) 



Let us suppose another complete orthonormal basis {\ipnm)i} in Fi, each of 
which maximizes the entanglement entropy with respect to the decomposition T\ = 
J-\+ ® T\-. Since both {\ipnm)i} and {\tpnm)i} are complete orthonormal basis in 
T\i they are related by a unitary transformation U in T\. Moreover, U is a unitary 
transformation in !Fx+i since any states maximizing the entanglement entropy are 
related by unitary transformations in !F\+ as shown in Appendix A. 6. Therefore, 
the orthonormal basis {\ipnm)i} defined in the phys ical p rinciple (c) is unique up 
to a unitary transformation in T\- and is given by (3.213). 
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Now let us show the equation (3.214). The right hand side is transformed as 
follows. 



2-Ki(j—k)n/N 



2+(k\n') 2+ C n 



jk \ n / mm'n 1 

x|0' + m)modN) 1+ (g> \ (g> \ {k + m)modN} 2+ <S> \m') 2 
= ^ (^=\(n' + m)rnodN) 1+ <g)\(n' + m)modN) 2+ 

mm'n' 

xC n ' m '\n')i- ® |m') 2 - 
= (^V|m"> 



)i+ <g> |m ) 2 + 



Vn'm' 



<8 |m'>2- 



(A.54) 



The final expression is |</>) itself. 
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